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THE TREND TOWARD AUTOMATION 
(IN ANTICIPATION OF THE 22ND CONGRESS OF THE CPSU) 


(Editorial) 
Translated from Avtomatika i Telemekhanika, Vol.22, No. 10, 
pp. 1265-1268, October, 1961 


In October of this year the 22nd Congress of the Communist Party of the Soviet Union is being convened. The 
congress will study the program for creating the foundations of communist society in our country. 


N. S. Khrushchev vividly compared communist society with a cornucopia of the products of physical and spi- 
ritual labor which will be accessible to every member of society. Everyone shall have enough food, clothing, foot - 
wear, housing, books. N. S. Khrushchev says that the time is not far off when people in our country will work 3 to 
4 hours a day on the basis of the further development of production, science and engineering, and on the basis of 
technical progress and the automatization of production. Men will control complex machines and mechanisms and 
in this short working day will produce appreciably more products than now; this will guarantee a fuller satisfaction 
of the requirements of the country and of all the requirements of each citizen of communist society. 


Automation will help millions of Soviet workers make their work more perfect and intellectual, and will bring 
the level of development and the preparation of the worker closer to the level of the engineer. The creative forces 
and energy of Soviet man in the field of production will become ever more concentrated in scientific and research 
work; it will include the planning and creation of new effective means for achieving material welfare, 


The reduction of working time due to automation is, as we know, the most important condition for true human 
freedom and will permit workers to devote greater attention to the development of their physical and mental capacities, 
to the enrichment of their spiritual surroundings, and to the all-round development of their talents. 


In the USSR automation is becoming the task of the entire people. Because it is the most important part of 
technical progress the automation of production has become the general line of development in industry and is the 
center of attention of the Communist Party. 


The scientific, engineering and organizational sides of automation form a unified whole. This makes it neces- 
sary to develop the proper forms of organization for scientific-research work and automation projects; it is also neces- 
sary to establish spheres of the effective application of automation and the preparation of cadres with the required 
qualifications. 


The practice which has existed heretofore of automatizing individual production operations or parts of tech- 
nological processes, notwithstanding the high engineering and economic effectiveness of individual projects, cannot 
assure the solution of the problem which has been posed. A qualitatively new result can be obtained only by in- 
ttoducing full automation which has not yet been widely developed. The basic obstacle to this development is the 
fact that many forms of existing equipment and technological processes are not suited to automation. In order to 
utilize all the possibilities and advantages of automation more fully and to make it one of the chief factors in the 
economic development of the country it is necessary to achieve a high level of production mechanization; the oper- 
ating machines and the technology itself must be developed in accordance with the requirements and possibilities 
of highly efficient automation. New technological processes and aggregates must be designed in such a way that the 
automatic control systems are not mechanical auxiliaries but organic parts of the process and aggregates. Automation 
must encompass all sections of production with a unified automatic control system that assures integrated operation 
of all aggregates in optimal modes, the complete mechanization of the technological processes, and a highly or- 
ganized flow of production. 


For this purpose it is necessary to employ new forms for coordinating the work of scientists and research and 
design organizations, One such form is the complex development of plans for proposed projects involving new auto- 
matic installations in the basic sections of industry. Such plans must take into account the interrelation of the prob- 
lems involving automatic control, the newest technology and equipment, a high degree of mechanization, and adequate 
Power resources; it must involve the joint efforts of engineers, machine builders and automation specialists. 
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The plans developed by leading specialists in various scientific and engineering fields must provide a clear 
concept of the automated factories of the future and must define those scientific engineering problems which must 
be solved so that the automated factories can be constructed in a brief time. 


At the present time the science of automatic control has achieved considerable successes. The general prin- 
ciples and laws which have been discovered for designing automatic control systems, as well as the engineering 
means for realizing these systems, have made it possible to introduce automatic control into industry and to assure 
not only an appreciable improvement in the engineering -economic production indices but to create a number of new 
processes which were inconceivable under conditions of manual control (atomic power installations, continuous rolling 
mills with a rolling speed exceeding the speed of an express train, installations for the annealing and gasification of 
materials in a boiling layer, etc.). 


However, the functions which a human being performs in controlling technological processes are immeasurably 
more expensive than those which modern automatic control systems perform. The automatic control systems of the 
future must eliminate the effect of human subjective factors on the engineering-economic indices of the process; 
they must process an enormous amount of information, analyze the flow of the process, and plan and program its 
control in advance. 


The design of such automatic systems must be based on the wide utilization of computer engineering techniques 
which will make it possible to remember a greater amount of data, to solve complex logic problems automatically, 
and to perform various types of calculations at a great speed. For comparatively simple deterministic technological 
processes we can build control systems which are program control systems. 


For the more complex processes for which it is impossible to formulate a control program in advance, a program 
will be sought automatically during the course of the process itself on the basis of retrieved information and ap- 
proximate data on the properties and characteristics of the process which is stored in the system. In such automatic 
systems self-adaptive and learning methods will find the widest application. Therefore it is necessary to expand the 
work on the theory and methods for choosing the optimal automatic control strategy under conditions of incomplete 
a priori information on the charactcristics of the objects and under conditions where noise is present. For this purpose 
it is necessary to develop and apply new mathematical techniques (algorithm theory, the theory of operations, games 


theory, information theory, etc.). 


It is necessary to expand work on the development of the theory and principles of designing self-adaptive and 
other automation and remote-control systems which are capable of adapting to changing external conditions, 


In designing complex automatic control systems the study of the characteristics and properties of the auto- 
matized processes acquires exceptional importance. At the present time such a study is being performed using ex- 
tremely primitive means which require a great deal of time expended by the researchers and do not yield reliable 
information. As a result it is sometimes necessary to perform a prolonged (sometimes lasting many years) refine- 
ment of the designed automatic control systems. It is necessary to develop the methods and the automatic equipment 
for studying the dynamic characteristics and energy properties of controlled objects, and to perform a statistical 
analysis of the processes which occur in them and the perturbations which act on them. It is also necessary to develop 
methods for determining the optimal parameters of new automatic machines, production lines, and complexes on the 
basis of fundamental investigations of the kinematics, dynamics and accuracy of the individual mechanisms and 
machines; further, we must develop the theoretical foundations of automatic technological and transportation processes. 


In the near future automatic devices which are used not only for automatizing production processes but also for 
automatizing elements of mental work in operations which involve appraisal, planning, design, and scientific 
research will begin to find wide application. 


In order to solve this problem it is necessary to develop the study of the mental activity of a human being in 
conjunction with physiologists; in particular, it is necessary to investigate the processes involved in the solution of 
problems by the human brain (in such investigations, of course, the examination of all possible variants is impossible). 


It is necessary to develop the principles governing the design and theory of programming control machines and to 
take up the investigation of the physiological features of the human being for the purpose of developing moreefficient 


control systems. 


In order to produce effective control systems we must find new principles governing the design of automation 
elements by basing ourselves on the newest achievements in physics and chemistry, including the use of physico- 
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chemical processes and reactions, intra-molecular and intra -atomic processes, nuclear and other types of radiation, 
etc. 


It is necessary to develop new and highly reliable automatic control devices which are based on the use of 
semiconductor, magnetic and other contactless elements. Pneumatic automatic control systems, which do not con- 
tain moving parts and use the interaction of a jet, must be applied as widely as possible. We must develop micro- 
elements for automatic devices, including logic and memory elements which have ultra-small dimensions and assure 
high speed of response, high reliability, and are easy to manufacture; on the basis of these elements we can develop 
microdevices which simulate nerve tissues, 


The combination of logic and computing operations in complex automatic control systems is being successfully 
solved through the application of digital techniques. Therefore it is necessary to assure the widest application of 
automation devices which employ digital engineering and feature the aggregate principle of design. This, in turn, 
requires the replacement of analog transducers by transducers which represent the values of the measured quantities 
directly in digital form. 


It is necessary to perform investigations which assure transition to a common state-administered system of en- 
gineering means for automation which is based on the wide utilization of digital engineering elements. 


We must develop the theoretical foundations for designing automatic control systems that assure processing of 
the retrieved information in a form that is convenient for the use of the servicing personnel. 


Machines are required for the automatization of the processes involved in synthesizing and analyzing the struc- 
tures of automatic devices, for planning and designing, for realizing complex programs for checking, testing, ad- 
justing, and assembling, It is necessary to develop units which perform the processes involved in analyzing and 
synthesizing the structures of complex control devices by machine. 


Remote control units which make it possible to incorporate machines and technological aggregates distributed 
over wide areas into a single controlled complex will be developed extensively (units for use at oil-well installations, 
irrigation systems, power systems, gas- and water-supply systems, transportation systems, etc.). 


An important condition governing the development of automatization is the assurance of operational reliability 
of automatic and remote control systems, This requires extensive expansion of work in the field of reliability. 


Here we need to develop both the theory of structural reliability and the theory of the reliability of devices 
and their elements as a function of their characteristics, operating conditions and modes of utilization. It is neces- 
sary to develop engineering methods for computing and improving the reliability of automation and remote control 
systems, including methods for constructing highly reliable devices using elements with low reliability. 


The achievement of the enumerated measures will be a great contribution to creation of the material and en- 
gineering base of communist society. 


Soviet specialists on automation are inspired by the new draft program of CPSU which will be discussed at the 
forthcoming 22nd Congress of the CPSU; they will apply all their creative forces, knowledge and energy to achieve 
the problem posed before them. 
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VICTOR SERGEEVICH KULEBAKIN 


(ON HIS SEVENTIETH BIRTHDAY) 


(Editorial) 
Translated from Avtomatika i Telemekhanika, Vol. 22, No. 10, 
pp. 1269-1272, October, 1961 


On October 30, 1961 the outstanding Soviet scientist in the field of electrical engineering and automation, 
Stalin Prize laureate, Academician Victor Sergeevich Kulebakin will reach his seventieth birthday and the anniver- 
sary of 45 years of creative scientific and pedagogical activity. 


Victor Sergeevich was born October 30, 1891 in the family of a rural teacher. In 1909 he entered the Moscow 
Higher Technical School where he specialized in the field of internal combustion engines and electrical engineering. 
In his advanced courses he heard a series of lectures by N, E. Zhukovskii on aviation and made a detailed study of 
the theory of heat engines; he also performed a series of research assignments in the field of mercury rectifiers and 
electrical machines. 


V. S. Kulebakin was an initiator and active participant in the foundation of new scientific establishments (the 
State Experimental Electrical Institute, later renamed the V, I. Lenin All-Union Power Institute; the Institute of 
Automation and Remote Control of the Academy of Sciences, USSR) and of new fields of industry and forms of produc- 
tion (instrument design, the production of special instruments, etc.). He participated actively in the development 
of the Lenin electrification plan for the country ("GOELRO") and performed tasks assigned by the "GOELRO" commission, 


A great contribution to the realization of the "GOELRO” plan was made by a large group of scientists among 
whom V. S. Kulebakin stood in the front ranks, 


It is precisely to participation in the solution of this great task of the economy of the country that one of the 
basic fields of activity of V. S. Kulebakin was devoted. 


Working in the electrical engineering laboratory of the Moscow Higher Technical School and later in the 
laboratories of the All-Union Electrical Engineering Institute and the Moscow Power Institute, he conducted a num- 
ber of research projects in the field of electrical machinery, electrical measurements, electrical equipment, elec- 
trical power, etc. 


In these fields V. S. Kulebakin has written scientific papers devoted to computations and further theoretical 
developments in the fields of stationary and transient processes in dc electrical machinery, synchronous and asyn- 
chronous electrical machines and transformers. He proposed a number of new designs for electrical machines and 
investigated new systems, For his development of an electrical mining locomotive with a capacitor induction motor 
V. S. Kulebakin along with other participants in this project was awarded a Stalin Prize. 


Of great scientific value are the papers written by V. S. Kulebakin on the theory of cross-field rotary am- 
plifiers, on the pulse method of controlling the speed of electric motors, and other problems, Important papers were 
written by V. S. Kulebakin on the kinetics involved in the excitation of electrical machines; these papers form the 
basis for a new branch of the theory of electrical machines which has become especially significant in recent years 
due to the wide use of rotary automation and the development of new large power systems. 


Much space in V. S. Kulebakin's papers is occupied by automatic control of electrical machines. He developed 
the theory and methods of computation for various types of voltage regulators. In his monograph “Electrical Equip- 
ment,” which has sold out in several editions, a systematic presentation of the methods for computing and designing 
starting and control equipment is given. His fundamental work on the computation and design of starting and con- 
trolling rheostats was awarded a prize by the People's Commisariat of Education. 


In a number of his papers V. S. Kulebakin developed the theory and methods of computing devices for the elec- 
trical ignition of internal combustion engines; these papers provided the first scientific bases for this field of equip- 
ment construction. In this field he is one of the greatest specialists, and his papers are still widely used in practice. 
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V. S. Kulebakin devoted a great deal of work, energy and initiative to solving problems in aviation electrical 
equipment; he studied this field from the very first days of his engineering and scientific career. His research en- 
compasses basic problems in this complex and varied branch of engineering; in particular, he dealt with problems of 
lighting (illumination for night flights, reflection of light from terrestrial coverings, the reflection of light from 
rotating propellers, the computation of the illumination of open spaces). 


A very important and essential part of the work performed by Victor Sergeevich Kulebakin is in the field of 
automatic control theory, His most significant work dealt with the perturbed movement of automatic control systems 
in the presence of external signals f(t), In the early period of development of automatic control theory V, S, Kulebakin 
pointed out the importance of the analysis and synthesis of automatic control systems while taking into account per- 
turbations f(t) which act continuously. In 1940 V, S, Kulebakin cited the general conditions for determining the basic 
parameters of automatic controllers in the presence of external perturbations. In treating the characteristics of tran- 
sient responses x(t) in automatic control systems V. S. Kulebakin first proposed the use of integral predictions for 


co) 
processes. In 1940 he established the relationship between the integral f x dt and the parameters of a control system 
0 
(in particular, their static accuracy). As we know, later on the use of the integral prediction method in automatic 
control theory proved to be exceptionally fruitful. 


In general control theory the integral predictions proposed by V. S. Kulebakin were developed for the case 
where the mean-square estimate of system accuracy is to be determined, so that these predictions can be used as 
criteria in the general theory of automatic control system optimalization. Victor Sergeevich Kulebakin was one of 
the few scientists and researchers who came forward 20 years ago with an energetic defense of frequency methods in 
automatic control theory. 


A natural development of the indicated work by Kulebakin was his work in the field of the general theory of 
automatic control system invariance. At the present time the general front of research in this field, both in the 
Soviet Union and abroad, is extremely extensive and V. S, Kulebakin is undoubtedly the leader of this scientific 
school in automatic control theory. V. S, Kulebakin was able to draw in the most outstanding scientists and specialists 
to work on his problem; these included Academician N. N. Luzin who established the mathematical conditions and 
the theorems of absolute invariance and invariance with an accuracy of up to ¢ in the theory of differential equa- 
tions, Later on V. S. Kulebakin used N, N, Luzin's theorem as the basis for proving that in a multiple-loop system 
that is subjected to the simultaneous effect of several perturbations it is possible to create conditions such that any 
coordinate x; becomes independent of one or several perturbations f(t) (i.e., the achievement of so-called selective 
invariance), In the field of invariance theory V. S. Kulebakin has achieved results such as the solution of the prob- 
lem involving the applicability of the principle of absolute invariance in real physical systems, the general nature 
of problems of autonomous control and selective invariance, further development of the theory of combined auto- 
matic control systems the forms and conditions governing invariance in various automatic control systems and its 
utilization in programming and optimalizing mathematical and control machines, In testing the problem involving 
the behavior of continuously perturbed automatized linear systems and depicting the functions by means of the in- 
tegrals of differential equations, V. S, Kulebakin produced the theory of the K(D) transform for functions satisfying 
the condition K(D)f(t) = 0 for which K(D) # 0; f(t) « 0. The indicated condition represents a special form of in- 
variance for automatic control systems that has a great practical significance. 


In summarizing this work V, S, Kulebakin demonstrated that the principle of total or partial invariance is ap~ 
plicable to real physical devices. The established invariance criteria make it possible to devise methods for the 
accurate computation of physical and structural parameters for these devices, and automatic systems designed ac- 
cording to the principles of invariance are refined and rational, In developing this point of view, V. S, Kulebakin 
frequently put forward the idea of designing automatic systems and production objects in direct conjunction, and of 
assigning the properties of controllability to the production object proper while altering these properties in the direc- 
tion of self-adaptation and automatizability. 


V. S. Kulebakin is not only an outstanding scientist but also a talented teacher who is able to organize a faculty 
tightly around him and to inculcate interest in and love for science in young people. 


In his long years of professorial activity at the Moscow Power Institute and the Zhukovskii All-Union Higher 
Aviation Institute Victor Sergeevich developed a number of new courses: a general course in electrical equipment, 
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starting and control equipment, automatic control of electrical machinery, aviation electrical engineering, the 
electrical equipment of airplanes, etc. V, S, Kulebakin also compiled textbooks and various other books on most 

of the subjects mentioned above. As an able and experienced teacher V. S. Kulebakin devotes a great deal of energy 
to the education of young scientific workers and students and supervises their work. Many of them have already be- 
come outstanding figures in the national economy, scientists and engineers, candidates and doctors of the technical 
sciences, professors, and Lenin Prize Laureates. 


Victor Sergeevich is doing a great deal of government and organizational work: for eight years he was a mem- 
ber of the Higher Degree Commission of the All-Union Committee of Higher Education of the USSR; he also headed 
the Methods Commission of the Power Engineering Universities (G.U.V.U.Z, V.S.N.Kh.), worked as a member of the 
Commission for the Mobilization of the UralResources during the Great Patriotic War, participated in the work of the 
Commission on Developing the GOELRO Plan before the war and on the Commission on 100-Cycle Current attached 
to the Council of Labor and Defense, and on the Commission of the Division of Technical Sciences of the Academy 
of Sciences USSR, on the electrification of railway transport. From 1940 on V. S, Kulebakin was a member of the 
Council of Engineering and Economics Experts of the State Plan of the USSR for many years where he participated in 
scientific-engineering planning for the development of the national economy. 


From 1951 on V. S. Kulebakin took an active part in the Academy of Sciences, USSR, in providing scientific 
help and cooperation to the builders of large hydroelectric power stations; he headed the brigade of scientists of the 
Academy of Sciences, USSR, on the Stalingrad hydroelectric complex. 


V. S, Kulebakin is an active public figure. From 1953 through 1959 he was a deputy of the Moscow City 
Council of Workers Deputies. V.S, Kulebakin is one of the leading organizers and active participants of a number 
of scientific conferences and conventions. He is a member of the editorial board of leading scientific journals of 
the USSR and is the editor-in-chief of the journal, "Izvestiya Akad, Nauk SSR, Otdel. Tekh. Nauk." 


The activity of V. S. Kulebakin as a member of the Bureau of the Division of Technical Sciences of the Academy 
of Sciences USSR and as a member of the presidium of the Editorial and Publishing Council of the Academy of Sciences 
USSR, etc. is well known. 


V. S. Kulebakin is the author of more than 300 scientific and engineering papers published in the form of books, 
monographs, textbooks, and articles. Many works have been published under the editorship of V. S. Kulebakin. 


The election of V. S. Kulebakin as a corresponding member of the Academy of Sciences USSR in 1933, and as 
an Academician in 1939 was an acknowledgment of his great services to Soviet science. 


In 1942 Victor Sergeevich was awarded the title of Major General of Aviation Engineering. The activity of 
V. S. Kulebakin has been marked by high government awards: two orders of Lenin, the Order of the Red Banner, the 
Order of the Red Banner of Labor, three Orders of the Red Star, two Orders of the Badge of Merit, and six medals, In 
1961 V, S. Kulebakin was awarded the honorary title of Distinguished Scientist and Engineer of the RSFSR, 


In honoring Academician V. S. Kulebakin on his 70th birthday the editorial board of “Avtomatika i Tele- 
mekhanika” wishes him long life and further success in his scientific work for the good of our country. 
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ANALYTIC DESIGN OF CONTROLLERS IN SYSTEMS 


WITH RANDOM ATTRIBUTES 
Il, EQUATIONS OF OPTIMUM SOLUTIONS. 
APPROXIMATE SOLUTIONS 


N. N. Krasovskif and £. A. Lidskii 


(Sverdlovsk) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 10, 
pp. 1273-1278, October, 1961 

Original article submitted March 18, 1961 


The general approach, as described in [1], is used to derive the equations which determine 
the optimum control rule. The method of obtaining approximate solutions is described. 


1. Equations of the Optimum Lyapunov Function v® and of the Optimum Control 


Rule —° 
1, The principles to be used in determining the optimum control rule were formulated in [1], the latter being 


such that it minimizes the integral criterion of quality“ 
co 
Te = \ M{olz(t), &(01/ 20. tes to = 0) dt = ming (1.1) 


0 





of the stochastic system under control 
dz; 
= = G [2;,---, Zn, H(t), EI, (1.2) 
E = & [a45--+) ns N- (1.3) 
In accordance with these equations, €° is determined by the condition 


[A+ |. = min, (“+ o| = (1.4) 





where v° is the positive-definite optimum Lyapunov function (see (4.1) and (4.2) of [1]. The proof is given in 
Appendix I. We shall now obtain partial differential equations for v’ and €° which follow from (1.4). 


2. The derivative dM{v}/dt at the point (x, n) can in view of the equations (1.2) and (1,3), be writtenas ** 


aM te) -> ad (2, &, 0) + | » (2, 8) de 9 (mB) — 


x (1.5) 
—4(n) »(2, 0) +} : Fea ay 5055 
jm 


* The concepts and notation used are as introduced in [1]. 
** In (1.5) the integral is understood to be a Stieltjes integral [2]. 
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In the particular case of n(t) assuming only a finite number of values n(t)={a,..., %,}, we have (at 
the point x, a ; ) 
dM) — 





a= 2 ae, —% (% Ba) + >? P,; (v (2, a,) — 
—v (z,4)) +5 > eae ky hg9,0 (1.6) 


In the case when the function q(a, 8) has a density p(a,g ), this can be written as 


ae -% dv &. n) 








Be. (2, & 0) + | » (zB) p(n, 8) dB — 


h 3} Mele) dv 5 Dk (1.7) 


—q(n) »(z,9) +5 2, k, io ph,6,6,,. 
j= he 


The derivation of the formulas (1.5)-(1.7) is given in Appendix II. 
3. The equation 





oe Mo, (z, 8 0) + | v(x, 8) dya (ny 8) — 


i=] 





—4(n) (2, + 3% ee Be kyl 8, +o I[z, §, 4] = 0 (1.8) 


(as well as its special forms which can be obtained from (1.6) or (1.7) ) provides the first equation for the deter- 
mination of v(x, n) and €°(x, n). The other equation for v° and €° is derived by differentiating (1.8) with respect 
to € because, according to the condition (1.4), the left-hand side of (1.8) attains its minimum at € = €°: * 


av (x, n) 99; (2, & 0) 8 (pH ;) 
p? E ~. 9 —— + +3 % ek k, 59,9; 7 +3 = 0. “— 





The optimal functions v" and €° satisfy, therefore, the equations (1.8) and (1.9), and moreover this solution 
is adopted where v’ is a positive-definite function. 


2. Arriving at Approximate Solution of the Problem 

The equations (1.8) and (1,9) prove too difficult to solve in the general case. This is also the case when at- 
tempting to find a solution which would ensure the bringing of the trajectory to the point x = 0 (i.e., such that its 
v will be a positive-definite function). The following method, however, which leads to approximate solutions is 
presented. Instead of the equations (1.2) and (1.3) we shall consider the auxiliary system 





= i [z, 7, 6), § = E (z, 7, 9), (2.1) 


where the introduced parameter % is such that when 9= 0 the minimum of the functional 


\ e (x, §, 0) dt 
0 


* Should €° be submitted to an additional constraint (for instance | €°| <1), the minimum in (1.4) should be reached 
with this condition also taken into account, 
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is attained in a straightforward manner and also when 9% varies from zero to unity, the functions » and ¢ approach 
in a continuous manner the functions ; and w respectively. The optimal Lyapunov function v(x, 7, 9) and the 
optimal control rule €%(x, 1, 9) should satisfy the equations (1.8) and (1.9) for every 9; therefore, differentiating 
the latter with respect to 9 one is able to obtain the equations which describe the variation of the solutions 9° and 
E° with 9. Availing oneself of this device may in a number of cases facilitate the solution. In particular, one can 
try to find v° and €° in the series form 


yo = 2 an (6) tx(z,y), &°= Dibs () Cx (2, 0) (2.2) 


in terms of some functions ¢, (x, 1) in such a way that the equations (1.8) and (1.9) are satisfied not exactly by 
the finite aggregate (2.2) but approximately (for example, in the best mean-square approximation sense). The con- 
ditions of such approximation lead to equations describing the variation of the coefficients a(®) and bj() with 
respect to . The method is also notable for the fact that if one starts with the stable solution of a problem for $= 
= 0 and then the problemand also its solution are continuously varied, one is able to obtain the same branch of 
the solution for the equations (1.8) and (1.9) which, for every $, will give solutions ensuring that the trajectories of 
the transient process attain the prescribed motion Z(t) (when x(t) = 0). 


3. Example 
A simple example will illustrate the described methods for obtaining a solution. Let us examine the selection 
of aregulator &° = &°(x, m) for a second-order system whose damping 7 (t) is known only probabilistica lly: 





d*z d: 
at nl) ate=mb+1(@) (m = const), (3.1) 


or in an equivalent form of simultaneous equations 


d. 
4 = 2, o = -—2,— qh (é) ze+ mE + 2 (z). (3.2) 


We shall confine ourselves to the simple case of n(t) assuming only two values: 9 = 7, and 7 = nN, with the 
transition probabilities P[nj--n; within time At]= pyjAt + o(At). The random disturbance y2 shall be characterized 
by the frequency A > 0, by the function p(x, —, 7) = 4% (u# = const) and by the random quantity vz whose variance 
is 02 ([1], §2).Itis required under these conditions to find the control rule 

Eo = E° (1, Za, 0) (3.3) 


such that whatever the initial conditions x49, X%9, 19 = (1; OF Nz), and tg = 0, a probabilistically stable transient 
process should occur in the systems (3,2) and (3.3) which will minimize the mean-square error 


I (zx0, 220, Me) = | M Ce} (2) + 23 (0) + 8 (es (9, 29 (0, 0 (0) dt = ming, (4) 
0 


We shall endeavor to find the optimal Lyapunov function when the latter is a quadratic form 
vo = by; (n) x? + 2by2 (q) xiz2 + b20(y) a. 


In such a case the equations (1.8) and (1.9) take the form 


2 (ian + Of) a2) an + 2 (OL) aa + OG) 29) (— 21 — ty ze+ mb) + 


2 2 
+ Pre dy Of — Of?) 2x +-200 p* 2h08 + >) a? + &* = 0, 


i, jm i=] (3.5) 
(t= 4,2; s = 4,2; Ls), (Of) = by (mp), 
(of) 2. + BM. a2)m + & = 0. (3.6) 
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By eliminating € from the equations, and afterwards by comparing the coefficients of the products X4Xj we ob- 
tain a system of quadratic algebraic equations for the coefficients bf YQ = 1, 2). 


APPENDIX I 
An outline is given for the proof showing that the optimal Lyapunov function v° (1,4) represents the solution, 
The asymptotic stability in probability [3] of the solution x = 0 under the condition that v° is a positive-definite 
function results from the theorems in [3]; these are proved in [3] when the function 9 (t) is a random function of a 
less general type; nonetheless this discussion remains valid for the variable 1 (t). 


We shall show that it follows from the condition (1.4) that (1.1) is minimized. It is assumed that the functions 
v and €° are sufficiently smooth to ensure the existence of solutions to the equations (1.2) and (1.3) and also to en- 
able one to use the derivative dM{v }/dt and to perform the transformations applied subsequently. By averaging 
over the random quantities x4(x9, 19, to, 1, €°), 1 (M9, to, t), we obtain from (1.4) the formula ° 


(* {v (z (t), 1 (t))/z0, No, to= 0} ) 
at ge 





== M (dM {v (x(t), 1 (8)} / dt)g./ 20, No, to = 0} = 


=_—M {o [x (0, go ()) / Zo, No, to}, @) 


which implies that M {v (z(t), 9 (#))/z0, No, te =0},, decreases with the increase of t. Integrating with respect to 
t from t= 0 tot = T we obtain 


M {v(x (T), (T)) / zo, No, to = O}g. — v (z0, No) = 
T 
san a \ M {@ [x (1), €°(0)} / 20, No, te = 0} dt. 
0 


From (II) we draw the conclusion that the integral is convergent when T— ©. But, in view of w being non - 
negative this means that lim M {w},o= 0 when t > ©, that is,that in accordance with the condition b) [1] the 
monotonic function M {p (x(t), 1 (t) / zo, No, to = 0},, approaches zero whent-» », Therefore, 


(it) 


v (xe, me) = | AM fw [x (0), 8°(O1/ 20, no, to = 0) dt, ail) 
0 


that is,the expression I, 0[Xo, 1 9] is finite, and I, 9 = v. 


We shall now use the indirect proof method and assume that there exists a function € * [x, n]#€°[x, n } such 
that when & = €* in (1,3) and (1.4), the solutions { x(t), n(t)}¢ * produce the inequality 


I ge [zo, No] < Je [zo, no] (Iv) 


for some initial conditions x9, 19 (at ty = 0). 
The condition (1.4) implies that 





(M2) > ta 


By averaging the inequality over the random quantities x(x9, 19. to, 2, €° ). 1(M9, te, t) and integrating 
subsequently with respect to t we obtain the. inequality 


M (v(x (T), (T)) / x0, No, to = 0}. — v (zo, No) > 
T 


=2- \ M {@ [z (2), §* (#)) / 20, No, to = 0} dt. 
0 


* The index €° signifies that the quantities are evaluated when the control rule is € = € °[x, n}. 
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As the integral on the right-hand side of the inequality is convergent when T-> « (in view of the assumption 
IV), by an analogous argument we conclude that the following inequality is valid: 


v (zo, No) S \ M (vo [x (t), &* (4)] / x0, No, to = Oh dt = J; [z0, Nol. (V) 
0 


But the inequality (V) is not consistent with the formula (III) and the assumption (IV), Thus the optimality 
of & ° is proved. 


APPENDIX II 
The fundamental idea with regard to the derivation of the formula for dM{v }/ dt is outlined. The case is 
considered when n(t) assumes a finite number of values only. The derivative dM{v }/dt can be obtained by proceed- 
ing to the limit, as described by dM{v }/dt= lim M { Av}/ 4twhen 4t +0, When deriving the formula, the 
terms of the order of smallness higher than At can be ignored. Therefore, if the quantities x(t) = x,(t)= 4 occur 
at the instant t, it is sufficient to take into account during the time interval t= 7 <t+ At the following events: 


A — that the quantity n(t) preserves its value aj, P(A)= 1— (Ss! P13) At; By—that the quantity n(t) changes 


i 
its value n(t) = aj precisely once, and at the instant t + At we have n(t¢+ At)h=a;, P(B)= At (p;;); 
C (v1,--+,¥,) — that within the time interval At an impulse p; vjd(i=1,...,9) occurs at the input, that is, the 


coordinates xj change step-wise by A,.xz, ~p,(z, n, &)¥;, P(C) = AAt. The events A and > B; are complementary, 


and the events C(v,,...,¥n) are independent of A or of Bj. We denote by RQ the common part of the events R and 
Q; by Q” the event complementary to Q; by Avo the increment of the function y along the trajectory of the system 
(1.2), (1.3) provided the event Q has taken place. We have 


k 
M { Av) = M (Avac) P (AC) + >) Avg, P(B) + Avgc+P (AC) = 
v ji 
Ov (x, a;) { = dv (x, a;) 
= u|>, air « 4.2, + zz > “Gz, oz, Az, WwALEY + 


7 8=1 s,l=1 


od " @v (x, a;) 
+ >} [v(@, a) — v(2, a) py; At+ eo @,(z,a,) At vp) 


jei s==1 


with an accuracy to the order of magnitude o(At), 


The formula (VI) is obtained in the following way. The second term is easily obtained and does not require 
any comments; the third term follows from the well-known rules which determine the increment of the Lyapunov 
function for ordinary differential equations [4] in view of the fact that for n = a4 = const the equations (1.2) become 
ordinary differential equations within the interval At. The occurrence of C indicates a step-wise change 4, x; , 
and thus when obtaining the first term of (VI) one has to determine the increment Av corresponding to some value 
of the random quantity vz and to evaluate subsequently the mathematical expectation M over the occurrences v7. 


The increments WV Aci, . ) are expanded in a Taylor's series up to the terms of the second order, neglecting 
seees nn. 


the variations in the coordinates x, , which are due to the fact that the coordinates x; up to the step 4x; and later 
do not remain constant but vary along the trajectory of the system (1.2)-(1.3) when 9 = a4. The terms of the third 
order or higher in the expansion of the increment Av, of the function y due to the step-wise changes 4,,x; in the 
coordinates are simply neglected because when 2 — cc, 6, + 0,Act = const [1] in the limit, these terms contribute 
infinitely small quantities of the order o(At), Because A,z, =p,v, and M {v,)= 0, M (v,v) = k,%% [1], 
therefore by dividing (VI) by At and by proceeding to the limit At - 0 we obtain (1.6), To get (1.6) in a more general 
form (1.5) we proceed in the usual way from a finite sum to the integral. 
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Original article submitted February 7, 1961 


4 linear system of automatic control which is optimal in a certain definite sense is considered. 

The equation for an optimum regulator is derived by applying a method of analytic design 
of regulators; the regulator may be put into effect by means of ordinary detective elements and 
regulating units. 


1. Introduction 

The dynamics of regulation of systems which are submitted to disturbing forces has continually drawn the 
attention of investigators. The trend to endow the self-regulating system with the ability to resist the disturbing 
forces has led to many concepts. There arose, in particular, the idea of controlling by load (and in the general case 
by disturbance) subsequently called the Poncelet-Chikolev principle; it became the subject matter of many in- 
vestigations. 





In order to clarify what is meant by these concepts it should be noted that any system possessing the property 
of asymptotic stability also has the property of resisting the disturbing forces provided the latter remain sufficiently 
small. This, however, does not suffice to be of much practical use. 


Firstly, it is important to note that even though the disturbing forces are usually bounded, their sufficient 
smallness nevertheless cannot always be guaranteed. Secondly, should these disturbing forces remain sufficiently 
small, but still measurable (the usual situation in the case of Poncelet-Chikolev principle), or should they be given 
in the form of known functions of time, it is of considerable importance to reduce, as much as possible, the domain 
of feasible deviations of the system from its undisturbed motion. Poncelet, Chikolev and some more recent in- 
vestigators have, in fact, had this particular goal in mind. It is our aim in the present paper to find the control rule 
in an analytic form and such that a measure which provides a description of the resistance of the system to the ex- 
ternal forces will assume its extremal value. 


The Poncelet-Chikolev principle, as well as the conditions governing its application to various types of con-" 
trol system, was investigated by V. S. Kulebakin, B. N. Petrov, A. Yu. Ishkinski, A. G. Ilvakhnenko, G. M, Ulanov, 
and others. 


In order to be concise the relevant literature is not cited here but can be found in [1, 2]. 


The author would like to express his gratitude to A. M, Letov for the formulation of the problem as well as 
for his aid in completing this work, and to N, N, Krasovskii for his most valuable comments. 


2. Formulation of the Problem 
An object under control is considered such that its disturbed motion is described in the neighborhood of the 
equilibrium position by the equations 





Te = Dd) bata + mE + fe (t) (ke =4,.....n). (2.1) 


a=l1 


Here n,, are the generalized coordinates of the object, & is the coordinate of the regulating device, fj,(t) are the 
continuous and bounded time functions describing the external disturbing forces, and by,, and my ate the constant 
parameters of the object and of the regulating device respectively. 
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Let the equations (2.1) be defined in an open domain N and the coordinates n}, satisfy the obvious boundary 
conditions 


Th (0) = Tho» Ne (0) = Noor +++ 5 Nh (0) on Tho? 


Ta (co) = M2 (0%) =... = 4, (0) = 0. (2.2) 


These signify that should 4 transient process occur in the domain N, it must then terminate whent = 
at the origin of the coordinate system. It is known that the conditions (2.2) can always be satisfied when the system 
(2.1) is closed and stable with fj,(t) = 0 [3]. 


For the sake of simplicity we shall confine ourselves to continuous dwindling* functions f,(t), that is,to func- 
tions fj,(t) = e5t a(t), where ¢(t) is a bounded continuous function, and 6 is a positive quantity, arbitrarily small 
[4-6]. 

Under these conditions we shall endeavor to find the functions £, 0;,...,. 1p of the class Cy related by the 
equations (2.1) and also minimizing the integral 


I(t) = \ Vae (2.3) 
0 
of the positive-definite quadratic form 
V = >) agni + cf. (2.4) 


k=1 


The integral for any t > 0 characterizes the disturbing effect of the external forces, and 1/1(€) therefore 
measures the resistance of the system to these forces. Other characteristics could, of course, be found to serve as 
another resistance measure. The measure presented here is closely related to the results obtained in [3]. 


Thus, the integral (2.3) has been selected to serve as the criterion of optimality of the system; it is a functional 
defined over the C, class of functions. The values of this integral characterize the integral square error weighted 
suitably by the constants ay, and c, the error being accumulated in the system throughout the duration of the transient 
process.” ° 

The problem consists in determining the equation of the regulator which together with the equation (2,1) of 
the object of control represents the stable system with the highest measure of resistance to the given external dis- 
turbing forces fj,(t). 


If such an equation does exist, it must obviously correspond to the principle of combined control as described 
in (1, 2}. 


3. Solution 

The problem under consideration is a variational problem of Lagrange. The required functions are E, 1j,....2n 
of the C, class such that they minimize the functional (2.3), and also satisfy the boundary conditions (2.2) and final- 
ly are connected by the relations 





& =%— Tee me —f,(t)=0 (k=4,...,2). (3.1) 


It is well known how to solve this kind of problem. The function 


* A bounded function is said to be dwindling if it approaches zero when its argument increases without bounds. 
** The choice of the parameters aj, and c is not discussed at this juncture. 
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H=V+)) i, (3.2) 


k=1 
is formed, where A, are arbitrary Lagrange multipliers satisfying the equations 


A, = 2a,n,— >) dA, (em ti---.m), 











a=] 
. (3.3) 
0 = XE — >) mA,. 
@=1 
If we add the equations (2.1) we obtain the complete system of the variational problem, 
The last (3.3) equation gives 
- m 
av a 
E= > 5 Nas (3.4) 
a=1 
Having put this — in (2.1) we obtain the equations of the variational problem in the form 
n n m 
My = Di bpate + My Di Gee + fy (O 
a=] a=l1 (3.5) 
n 
» = 2a,n, — » babe’ 
a=1 
We solve this for n; and Ay. The characteristic determinant of the system is obviously 
m? myn 
bi — B, » On, 3 a 
mm m? (3.6) 
bas , , ban B, 2 : ’ ’ a 
2a1, ’ 0, —r by a ’ a ba be 
0, ’ 2an, aes Din ’ " ae ban —p 
It is easy to show that if #y,..., My, ‘are the simple roots of the equation 


then the quantities —;, . . . , —ip must also be simple roots [3]. We shall assume that #;,..., Mn, are, in fact, 
simple roots, numbered in such a way thai the inequalities 


Rep, <0 (k=1,...,2) (3.8) 


take place and thus the first n roots have negative real parts. 


This can always be done in view of the previously mentioned property of the determinant (3.6), The general 
solution of the homogeneous system obtained from (3.5) will consist of a linear aggregate of exponential functions 
of the kind Cye# k*, Cn+Ke Uk (k= 1,..., m) and will contain 2n arbitrary constants cy, . . . » Cn» Cn +4+-++ Cn? 
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n n 
N= dy Akcee*#! + D) Anpecnpne 


(3.9) 
n n 
he = > A;** c,er et + >) Anticnsse** (a=i,..., n), 
k=l k=1 
where Ay (v,j=.1,..., 2m) are completely determined constants in terms of by, My, ay, C. 
In order to determine the general solution of (3.5) we apply the method of variation of parameters: 
n n 
>) Aicne**! + >) Adanengee ** = f, (0), 
k=] k=1 (3.10) 


n n 
Dd Artec + Ant ienee =O = t,....a). 


k=1 k=1 


With the assumptions made with respect to 1, one is able to determine the coefficients c,,, cp ,, by means 
of the formulas which contain indefinite integrals only, that is, 


= t+ DD, Ses, (Ode, 
i=1 (3.11) 


Cag @ Ni? DD Nes (t) dt (A= 1,...,8). 
i=1 


The Dj, Dp 44(i = 1,... ,m) can be found from the system (3.10) with the aid of Cramer's formulas; they 
will contain the coefficients Ay (v,j=1,..., 2m). Consequently, the general solution of the nonhomogeneous 


system (3,5) is , cong 
n= 2) Aaree’* + >) D AtDe™ {ere f, (t) dt + 
k=1 


+ > Ansatnpxe + 3) D Ange Dayse Ley, (t) at, 


k=1 k=1 i=1 


n n n : 2 
he = > Ayt*y,ere + PS > Ant" D,er" jeve f, (t) dt + o 
k=1 k=1 i=1 
+ D Artirngne™ + DD AMD et Sef, (Q dt (= N-m 
k=l k=] i=1 
In accordance with the condition (2.2) at infinity we obtain 
Tory = 9 (k=1,...,n). (3.13) 


In view of the assumptions made with regard to 1, the constants y}, can be chosen so that they agree with 
the initial conditions, the choice however not being unique. 
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Taking into account (3.13) we obtain from the first n equations of the system (3.12) a relation in the form 


(TY. + > D,\e*" f, (t) dt) et! 


i=1 


and we substitute it into the remaining (3.12) equations. 
It is assumed that it is possible to perform this operation. As a result we obtain 


= > Bin + > Baye >) De \en'n, () dt @=1,-.-.M). (3.14) 


k=1 k=1 i=1 


The coefficients BE. BO .,(K «= 1,...,m) are completely determined, and the D; (i= 1,...,m) are 
finally expressed in terms of the parameters by,, and my of the system, and in terms of the weighting coefficients 
a}, and c of the functional (2,3). 


As a result of this substitution of (3.14) into (3.4) the required equation of the regulator is obtained satisfying 
the Poncelet-Chikolev principle: 


= Dpnet+ D Pare Der! {ere f, (t) dt. (3.15) 


a=1 a=1 i=1 


This equation together with the original equations (2.1) describes a stable system of automatic control when 
f(t) =O(k=1,...,m). aes ttald i. tie of the found optimum system has a linear 
aggregate of exponential functions c,e" k' (k= 1, ..., m) as its general solution where 'y, have the property (3.8). 


The obtained results can easily be generalized for the case when the control is effected by means of n regulating 
devices; in such a case the equations of the variational problem are 


Me = Dy Ogata + Di Mya’ + hy (O 


a=] a=l 


i, = = 2a,n, — s Dada, (3,16) 


a=1 


- 2ctE, — >} Makha (k=1,..+,M). 


a=] 


By solving them in a similar manner we obtain the equations of the regulating devices in the analytic form 


7 > pn, + s pi), >} Dea! |e" f,(t)dt (k=4,...,n). (3.17) 


a=l1 a=] t=1 


The coefficients p"), p() _, By (a, k, i= 1,..., m) are evaluated by the step-by-step method and are final- 
ly expressed in terms of the system parameters by, ites and of the weighting coefficients ay, c) of the minimized 
functional, 


As can be seen from the obtained equations (3.15) and (3.17), the equations of the optimum regulator contain 
the external disturbances fj,(t) under the indefinite integral sign. As previously mentioned the integration constants 
must vanish [see (3.13)], As noted by N. N, Krasovskii the latter can always be attained provided that the functions 
f(t) are known over the whole interval t(0, ), In this connection it is of interest to consider two cases: 
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a) The functions f,(t) are known, The indefinite integrals of the known functions can then be found, and the 
optimum regulator can be designed in accordance with (3.17). Such a regulating device is considered later in the 
sections 5 and 6. 


b) The functions f;,(t) are not known but the values of the functions f(t), f}.(t), fk (t), + +++ can be measured 
at any time instant, 


In this case, after successive integrations by parts, the expression (3.15) becomes * 














[= -> PaNa + 3 Pn+e > D; py ras nh (t) (r =0, 1, 2,...). (3.18) 


a=) a=1 i=1 


The character of the transient process depends on how the index r is chosen, The bigger it is, the nearer is the 
transient process in the system to the optimal one in the sense of (2.3). 


The choice of r should be governed by the specific possibilities inherent in the design of the given automatic 
control system, 
4. Verification of Boundary Conditions 

When the equation of the regulator is (3.15) the coordinates ny, (k= 1, ..., m) of the controlled object vary 
in accordance with the formulas (3.12) and (3.13). 


Of course it can be verified straightaway that the initial conditions are fulfilled when the 1, are simple roots. 
We shall now verify the conditions at infinity. It is sufficient to show that the expression 





e-t | eutf, (t) dt (4.1) 


is a dwindling function of time for any k= 1,..., 0. 


When formulating the problem it was stated that the functions f,,(t) (k= 1, ..., m) are dwindling functions, 
that is, 


fe (t) = e*qpy (t) (k = 1,..., ), 
where | y(t)| = L, L is a positive quantity, and & is positive and sufficiently small, Then 


e-w! lef (t)dt| <L owt | etme | - more (4.2) 


Consequently, the expression (4.1) represents a dwindling function. 


5. Example 
A system of the first order whose variational problem equations are 








R= bn + mE + f(), 4 = 2an — dA, 0 = WE — mda (5.1) 
will be examined as an example. 
The corresponding characteristic equation is 
pt —o— "= 0, pe = F ep me (5.2) 
and the regulator equation is 
ga Pg — tl my (oat. (5.3) 


+ fe) denotes the rth decivative with respect 10 5 of f(t). Of course, the series 3 fet ff (t) is convergent 
prt 


for any t = 0. 
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In the case of f(t) = 0 the equation (5.3) is identical with the one obtained in [3], 


This example can also serve to demonstrate a different approach to the solution: use is made of a fun- 
damental concept in the analytic theory of differential equations; that is,that, with some restrictions on f(t), the 
latter may be considered as a particular solution of another differential equation [7]. 


For the sake of simplicity, it is assumed that f(t) = ¢ satisfies a very simple equation of the form 


FC 9=0 CM=/O, (5.4) 


which can now be treated as an equation of an element connected to the object under control. 


As, by assumption, the function f(t) is a dwindling function, we have ¢ (#)= 0. Then (5.1) is replaced by the 
following equations of the variational problem: 


n= bn ++ mt, FC, F 0 =0, 4, = 2en—dA,, 
(5.5) 


d oF oF 
a(*%)> —A+ ha 5E 0 = 2cE — my. 


The boundary conditions have already been discussed, The first, third and the fifth equations can each be in- 
tegrated independently of one another. Under these conditions the equations for the optimal system can now be 
written in the form 


n=bn+C+mt, Ft, & o=0, 
(5.6) 





gm Sg SS ool este, 


m 


The last formula in the regulator equation may be looked upon as a controlling signal indicating the deviations 
from the equilibrium position of the adjoint element £, which can be used to increase the measure of resistance of 
the controlled object to the disturbing forces. 


6. A Simple Self-Adaptive System 

Many tvpes of self-adaptive systems can be found in literature [8], We shall describe a simple case shown dia- 
gramatically in the figure below [16). 

The object-regulator part of the scheme operating in the ordinary manner, without adaptation, shall first be 
considered, If the characteristics of the object are known and are known to be stable, a regulator can be designed 
which will provide the system with an entirely satisfactory measure of control. 
However, it often occurs that the parameters of the object are not stable and vary 



































; Dx 4 in a haphazard way. In the simple case, for instance, when the impressed motion 
of the object can be described by the following equations: 
é * 
ts x 
: ‘ Ty = Dy byaNy + mE = A..--+m), (6.1) 














1) Controlled abject; _ 2) where 1, are the generalized coordinates of the object, € the coordinate of the 


— = pen ag controlling unit, and the parameters byq of the object and my, of the controlling 
anslyzer ; unit may vaty within a given bounded set D of real numbers [9]. The parameters 


by, and my, vary in a random manner but in such a way that they remain constant 
within a sufficiently long time interval, The probabilistic characteristics of this 
process remain unknown, 

For any fixed combination of by, Mm, €D the matching parameters pj, pz, . .. » Pn Of the regulator can 
be found such that the system as a whole has the required characteristics. The latter, however, shall not be main- 
tained if the parameters b;,,, and my vary in the set D, Under these circumstances there arises the problem of adap- 
ting the system to the variations of by,, and m x, that is,of having a new set of matching parameters of the regulator 
such that the whole system possesses again the required characteristics, 
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The adaptation process depends primarily on how the characteristics which the system should possess are form- 
ulated, This may be done as follows, We assume that a standard [9] of the object is defined in the set D, and that 
it is introduced into the system in the form of an electronic analogue corresponding to the equations 


Nk La 2 been (A = 1,..., n). (6.2) 
We also assume that the equalities 
nix (0) = n, (0) (k = 4,...,n), (6,3) 


are true for the analogue, that is,that themotionsof the object and of its analogue start from the same initial devia - 

tions of the regulated coordinates. It seems feasible to suppose that the analogue of the standard satisfies the require- 

ments of the Lyapunov asymptotic stability, and alsothat the transientprocess 1 satisfies the previously for- 

mulated requirements of the standard. An analyzer is also connected to the analogue which compares the transient 

processes 1’, of the standard and n}, of the real object of control. The comparison can be carried out in several 

ways. Let us consider the case when the analyzer produces such an equation of the regulator that the least value of 

the integral 

co n 

1=([ Dd, (, — apt + et] (6.4) 
5 k=1 

is attained, 


The main difficulty lies in determining such an equation. We now assume that in some way we have gotten to 
know the values of by, and my. | 


The new variables 9}, are introduced,defined by the formulas 
,=n, — (Aa i,..., a”), (6.5) 


With the aid of the formulas (6.1), (6.2) and (6.5) one can find the equations describing 9: 


6 = 310.0, + mE + Di (ee — die) Ne (R= 1... 2). 6.6) 
a=] 


The equation (6,6) differs from the equation (6.1) by the presence of the last component. Consequently, in 
finding a regulator which minimizes the value of the integral (6.4), one has to submit the object of control to con- 
stant disturbances equal to 


fy) = >) Ora — dea) Ne (= 1,... 0). (6.7) 


a=] 


The functions fj,(t) must be dwindling in view of the fact that the functions n*,, , which are the solutions of 
the system (6,2),are stable by our assumption. Consequently, the equations of the required regulator can be written 
in the form (3.15),where fj,(t) are determined by (6.7), 


The parameters by, and m; may be determined experimentally. To this end one can avail oneself of the 
methods developed in [10-15]. 


The system is in such case submitted to tests, and the information obtained yields the values of by, and my. 
After each process a correction should be introduced into the equation of the regulator. 


This process can be made to be automatic by means of a computer working out the values of the parameters 
with the aid of an agreed algorithm for evaluation from experimental data of the parameters of the object. 
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The problem of optimal control is considered for systems with distributed parameters, the 
motion of which is described by nonlinear integral equations. A theorem formulated as a 
maximum principle is introduced which can be used to determine the extremals for any 
general variational problems, where it is desired to minimize an arbitrary functional with 
supplementary boundaries. 


The essential and most general results in the theory of optimal automatic control systems [1, 2] have been ob- 
tained for objects of control with lumped parameters, the motion of which is described by a system of ordinary dif- 
ferential equations 


z= f (x,u), (1) 
where the matrix function 
f(z, u) 
f{(z,u) =|] : , (2) 
}" (z, 4) 


of the variables x and u is defined on the direct product 
(x, u)EX*Q, czEX", weEQ. (3) 


Here X" is an n-dimensional phase space, and Q an arbitrary, topological, Hausdorff space of the permissible 
values of the control parameter u. The optimal control problem is that of choosing from the class of permissible 
controls (for example, in the class of measurable, bounded functions, or in the class of piecewise continuous functions, 
with values from ) a function u(t) (tg = t = t,), for which the control trajectory x(t) for Eq. (1) traverses a path 
from the point xX» = Xg(tg) to some manifold M,. The dimension s of this manifold should not exceed n—1, and the 
integral 


P= { (x, uae (4) 
te 


along the path should be a minimum, where f° (x, u) is a known function of its arguments, defined on the same set 
asthe function f(x, u). 


In practice, however, objects of control with distributed parameters occur much more often than objects with 
lumped parameters or objects that can be reduced to the case of lumped parameters. It is therefore of great im- 
portance to obtain an optimal control theory for systems with such distributions. To this class belong, for example, 
continuous production lines, and also a whole series of objects involved in production. 








In many cases, a design problem consists of the choice of a certain definite optimal distribution of parameters 
of the assembly to be constructed. For such systems, a control system must be constructed which will give the best 
possible performance in a definite sense, when certain conditions are fulfilled, When it is a question of design, then 
the optimal distribution must be determined of those parameters that are at the choice of the designer. 


We will consider an object of control, the state of which is characterized by a one-column matrix function 
Q'(P) 
Q"(P) 


where P is a point in a certain m-dimensional region D of the Euclidian space E(y;, .... Ym). 
The matrix function Q = Q(P) will be called the state of the system. 


Objects having parameters distributed in the space can be affected by certain controlling actions, also dis- 
tributed in space and time. In other words, the controlling action is described by the vector function 


u=u (S), (6) 
where the point S € D. 
Each component of this control vector has the form 
Wag = Ugg (Ys - +++ Yq) = Ug (S,) (@=1,2,....m; B= 1,2,....94), (7) 


where p,, is a positive integer dependent on a. 


Thus the vector u = u(S) consists of 
m 
rs >) P. (8) 
a=] 


components. Here the functions u,g (S,,) are defined only on the set of points S, = (yy,..-» Yq) € N,CcD 

(a= 1,2,...,m), where N, is some manifold of dimension a. We will assume that the functions u,, g are measur- 
able, bounded, and square-integrable functions of S,, and that their values are in some admissible set 2, which can, 
in general, be any topological Hausdorff space. In particular, an important case is when Q is a closed region of some 
r-dimensional space Eo. 


We call each vector function u = u(S) taking values in the space 2 a control. 


The process of controlling the object consists of finding a state Q = Q(P) to correspond to each control u=u(S) 
according to a definite law. The law defining such a correspondence can be specified by a certain operation A: 


Q= Au. (9) 


In other cases, the law of correspondence between u and Q can also be given by more complicated functional 
equations, 


A wide class of objects of control can be described by nonlinear integral equations of the type 
Q (P) = \K (P,S, Q(S), w(S))ds, (10) 


where Q(P) and u(S) are defined by (5)-(7), and K = K (P, S, Q, u) is a one-column matrix function of the four vector 
variables: 
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K' (P, S, Q, u) 


(11) 
K (P,S, Q, u) =| - ° 


K" (P, S, Q, u) 


The functions K!(P, S, Q, u) (i= 1, 2, . .. , m) belong to the class Lz, and have the continuous partial derivatives 


OK’ _ _K'(P, T, Q, u) 
aQi aqi 
We now introduce the square matrix 


(j =1,2,...,”)- (12) 





aK _ |\ax* 3 
70 = [505 | (i,7 =1,2,..., m)- (3) 
The superscript i, here and in the sequel, indicates the row number, while j gives the column of the matrix. 
In this case, the problem of optimal control with distributed parameters can be stated as follows. 


On a set of states Q = Q(P) and controls u = u(S), related by the integral equation (10), let there be defined q 
functionals having continuous gradients: = 





r= '( Q(P)) @=0,1,2,....0, (14) 
P=P(Q(P),a(P)) =O (2) G=I+1,....9, (15) 
where 
- | ris, Q (S), u (S)) ds 
res es =\F(s Q (S), u (S)) as, 16) 
s* {Fs Q(S), w (S))ds 
F°(S, Q, wu) 
F(S, Q, u) = 
(17) 
F*(S, Q, u) 


The functions @(z)(i = 72 +1,..., 4) and Fi(s, Q, u) = 0, 1,..., k) are continuous, and have continuous 
partial derivatives with respect to the components of the matrices z and Q respectively. 


Here the functionals 144 =7 +1, ..., 4) can clearly depend on the control u. 


We note that the functionals H(i = 0, 1, .. +,» @) can also be determined when the point P does not necessarily 
run through the whole region D, but only through a subset (possibly of lower dimension) of this region. 








Among the permissible controls, we wish to find a control u = u(P) (u € Q), for which 
r=0 (i=0,1,...,p—4, p+4,...,9), (18) 


with the functional IP taking its least possible value. Here p is the index of the functional, 0 =p <q, 
A control u = u(P) yielding the solution of this problem we call an optimal control. 
We introduce the following rectangular matrices: 


a® || ao’ 


de Hoa] HOt. HOH... ey ts 











i 
70 — | 50 (i=0,4,...,4; $= 1,2,...0), oe 


grad J =|jgrad;J‘| (@=!+1,...,9;/=4,2,...,n). (21) 


Here grad I denotes the jth component of the vector grad I! in the coordinates Qi. 


The solution of the above-formulated problem can be based on the following theorem, the proof of which is 
given in Appendix 1, 


Theorem. Let u= u(S)(u € Q) be a permissible control, such that, in view of the equation (10), the conditions 


f=uQ =6,1,...,p—%, p44... 9) 


are satisfied, and let the matrix function M(P, R) = IIMgj(P, R)|l Gi, j= 1, 2, ..., m) satisfy the integral equation 
(linear in M(P, R)) 


aK (S,R,Q(R),ul(R 
M (Py Ry + gee ODS (ae PNET Q (A), ¥(R)) og — 











a 5. (5). ¥ (5) 1 (S, R) dS. (22) 


Then, in order that this control u = u(S) be optimal, it is necessary that there exist non-zero, one-column, 
numerical matrices 


a=|o,a,...,c}, 56 =f Croy---s teh (23) 
with c, = ~1, such that for almost all fixed values of the argument S€D the function* 


ae | | F (P, Q(P), u(P)) a) : 











IT (S, u) = a= 
x (50 O20) 8) & (PLS, Q (S)s 2) = (mer, R) K (R,S, Q(S),u)dR)+ 
20 / f F(P,Q(P), u(P)) oP 
+4 3 F (S, Q (5), 4) + 
+ b (grad J (Q(P)), K(P, S, QS), 2) 7 M (P, R)K(R,S, Q(S)u)dR) (24) 


of the variable u ¢€ Q attains its maximum, i.e., for almost all $ €D , the relation 


(25) 
holds, II (S, u) = H (S) 


* Two expressions in brackets, separated by commas, denote a scalar product. 
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The function H = H(S) is defined to be the upper bound of the values of the function 1(S, u) for u€Q and a 
fixed value of S: 


H (S) = sup II (S, u). "2 (26) 


The above theorem can be used to obtain a system of equations which is satisfied by an optimal control. 
This result is more general than the result obtained in [4]. 
The following problem, in fact, is posed in [4]. The object of control is described by the integral relation 
uN 
Qe = ot) r K' ((t, t, ui (t), «5 Uy (4) ae. (27) 


It is desired to find an optimal control u = u(t) = (u;(t), . . . , up(t)), t» St Sty, ué Q, where Q is a closed 
region of r-dimensional Euclidian space, such that the functional 


t 
ge =| F(x, Q(x), w(t) de (28) 


takes its minimum value for the condition that the point Q(t,) €M, where M is some manifold of arbitrary dimension 
§ not exceeding n—1 in the n-dimensional Euclidian space E(Q’, .. . , @®). 


We will show that Theorems 1 and 2 of [4] follow from the above-stated theorem. 
Let the equations of the manifold M have the form 
{ { n { 
@ =@ (%,.-.-, A) =@(Q) =1,2,...,8—8), (29) 
where the functions g!(Q) are continuous, and have continuous partial derivatives. We introduce the matrices 
8-35 
a 


K (t,t, u) =| K' (t,t, u)} (= 1,2,..., a). 








(i=1,2,...,8—s8, f=1,2,...,a), (30) 


(31) 
We also define I and 8(z) = z’: 

I? = @ -| F(t, Q(t), u(t) de = © (z) =| F°(, Q(t), & (%)) dr, (32) 

O' (2) =@'(z) (6 =1,2,...,n—2), (33) 


In this notation, the last problem can be stated in the following way. Let the region D be an interval [to, t,], 


and let the point S = r €[tg, ty]. We wish to find an optimal control u = u(r), t) = 7 Sty, u €Q, such that the func- 
tional 


I* = @® (z) “|r (t, Q(t), w(t), w(t) de 


has a maximum for the condition 


t 
I* = @' (2) = (9 ()= o'(\ x" (t, tu (x) dx (i=i,2,...,), (34) 
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Here the functions F(r, Q, u) and Ki(t,, r, u) play the role of the functions F'(s, Q, w) in the equation (16). 


Before we apply the theorem, we note that, since the matrix function K(t,, 7, u) in this case is independent 
of the variable Q, the only matrix function satisfying the equations (22) will be the function 


(M(P,R)) =0, PED, RED. (35) 


Applying the theorem stated above, we find that for the existence of an optimal control u = u(r), it is neces- 
sary that there exist a non-zero, one-row, numerical matrix 


a =| Co, C1,.- +» Cre (36) 


with cg= -1, such that, for almost all fixed values of the argument T the function I(T, u), having the expanded form 


i on 


II (t, u) = co\ >) OF QO.) et (t, t, wu) dt 4 





1G ag 
ae aat ts), 
+ eo F(t, Q (x), a) + 2 Dee eK! (t;, T, u), (37) 
j=1 i-«1 
gives the maximum for the variable u€Q. 
We now set 
n—s F 
gE, = nae (jit) 2>..., n), (38) 
i=1 


The last equality means that the vector — = (E;, ..., &,) must be orthogonal to the manifold M at the point 
Q(ty) €M, i.e., at the point Q(t,) the transversality condition is satisfied. 


The conditions (36)-(38) obtained above form the results of [4] stated in the theorems 1 and 2, 


We now show how to apply the theorem obtained above to obtain the conditions for a maximum and for trans- 
versality, obtained previously in [1], for the problems (1)-(4) formulated in the present article. 


We note, first of all, that the following lemma, proved in Appendix 2, is valid, 
Lemma. Let the matrix functions pw = p (t, r) and vy = v(t, r), t and r €[t, ty) be related by the equality 


u(t, ) + v(t, v) = |p (4.0) ¥@, 2) ao. (39) 


If the function v = v(t, r) is defined to be 


v (t) for 2 >t>t>b, 
0 for iu >tS>t>b, 


v = v(t, t) = (40) 


where D (1) is square-integrable in [tg, ty), then the function p(t, r) satisfies, forty =t = r = tg, the equation 


t 
p(t, t) + v(t) = \u (t, 6) dO v (vt) (41) 


t 


and for ty =T>t2ty 


t,t) = 0. 
p (t, t) 42) 
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We now turn to the derivation of the conditions for a maximum and for transversality, obtained in [1]. 
Let the equations of the manifold M in the space X" be 


g' =g' (z) (=1, 2,..., s—s). (43) 
The region D is the interval [t, t;], and the points P and S are given by 
P =t€ Ite, tilandS = t € Ito, ts). (44) 


We write the differential equation (1) and the equality (4) in the form of an integral equation 




















z() = \ T(t, t, (x), w(x) de, (45) 
where 
= 
F(t, t, = (x), 4 ()) = 1 hog oes (46) 
0 fri >t>tpt 
(i =0,1,2,..., m5 28 =0). 
If we now set 
2° (t) (47) 
Q (t) -_ —_ z (t), 
x" (t) 
K(t% Qe) = IF (tL = TU 4 Hu) C=O tm = 
= 9 (2) = f° (x (x), u(t) de = © (2) = \P (x, Q(x), u (x) dr, (49) 
t, 
gp (2) = 2°, (50) 
ty i 
f=q (Vy (x(x), u(x) de) = O'(\ F (x, Q(x), 4 (X))dr) 
Te eer Be (51) 
then with this notation the function I(r, u) defined in (24) will become 
t 
G9 (f 1 (=), u(e)) at) 
Il (x, u) as me | stews (7, t, Zz (x), u) Tr 
29 (f f(@ (t), u (¢)) at) ra 
-\u (¢,0) #10, «, © (x), w) dO) dt) + a I(t ©, (x), u) (62) 
Thus, using the lemma, and taking the expression f (z(t), u) +--> out from the bracket on the right, 
we obtain 
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(+, w) = 9 (x)((f (2 (x), t+), (63) 


where 


t; 
dq (ff (x(t), u (ty) at) 











~ (1) = a—— x 
ty 
are ca t 20 (J /(@(), 4) at) 
z ® e 
x(% (1—| M (¢, 0) dd) dt + a5. " (54) 


It only remains to show that the matrix function » (7) satisfies almost everywhere on the interval [tg, t,] the 
differential equation 


o=—yZ. (55) 
We calculate the derivative of the function »(1r), using (54): 


ty 
aq (7 (x(t), u(t) dt) t, 
Oz 


p(t) =a \ She. 6) 4 (¢, x) dt — 





t 


t 
t), t)) dt 
0 (fe u (t)) ) ote, wie) tse) 


=e 5 Oz 





According to the lemma, the matrix function M(t, r) for tj =>t = r2=tg is determined from the condition 


ty 


M (t, t) + of(2 (e) u (t)) a M (t, 0) dé 4f (2 (e. u(x) 








(57) 


Transferring the term gue. “(t)) from the left to the right side of the equation, and then taking it out 
as a common factor, we obtain 





M (t, t) = (\ (t, 0) de — 1) FEC) ui), (58) 








When we substitute this expression in (56), take the term — GLE “ (*) “out from under the sign of in- 
tegration, and use brackets, we find that almost everywhere on [ts, ty] we have 


(2) = - 





4 
dq (f f(z (), u(t) dt) 
lost \ ee “ M (t, ®) d0) dt +. 


t 





ag (j t (z(t), u @®) "| 


8f (2 (x), & (*)) (59) 
+a=- is . 


But from (54) we see that the expression in the square brackets is none other than »(r), Thus (55) has been 
proved. 
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It is evident that the maximum of the function H with respect to u, considered in [1], coincides with the maxi- 
mum of the function [II(r, u) in (53), and the function y (r) in this expression satisfies the differential equation (56), 
Thus the maximum principle for the problem stated at the beginning of this article has been proved. 


Moreover, if in the expression (54) for the function » = »(r), we set r = t;, then we obtain 


4 
dq (ff (x(t), u(t) dt) 


(4) =a—; (60) 





This means that the vector »(t,) is orthogonal, at the point x(t,), to the manifold M defined by the-equations 
(43). 


We have thus proved the transversality condition obtained in [1]. 


As an example of optimal control of an object with distributed parameters, we consider the problem of heating 
a “thick” body with zero initial temperature, through a thickness x(0 = x = 1) for a giventime T. The condition 
must be satisfied that the deviation of the temperature distribution Q = Q(x, t) (0 =t = T) of the body from a cer- 
tain given distribution Q° (x) at time T be minimal, in the sense that the functional 


1 
r= \(Q(@)—Q (2, T)I"de (61) 


isa minimum. Here y is a number greater than zero. 


In this case, the controlling parameter u = u(t) (0 =t = T) is the temperature of the external medium which 
is heating the surface of the body. 


It is natural that the temperature of the heating medium u(t) be a bounded function, i.e,, Ay = u(t) SA, (0 = 
st = T), where A, Az are given, that is,the set Q in this case is the interval [A;, Az]. 
As is known, the temperature distribution Q = Q(x, t) (0 =x =1,0 =t = T) can be given by the relation 


T 
Q(z, t) = \ K (zx, t, t)u (x) dt, (62) 


where K(x, t, +) is a known scalar function. * 


Since in this case, the kernel K(x, t, r) is independent of the controlled function Q(x, t), then, according to 
the equation (22), the function M (t, r) = 0 for any t and r in [0, T], According to the theorem (p, 1166), the 
optimal control u= u(r)(0 =r = T) at any fixed instant of time r must yield the maximum of the function 


1 
Il (x, u) = cy a7 (z, Q(z, T)) K (x, T, t) udz (63) 





or, in the notation of the last problem, 





1 in 
I (x, u) = col SIC) — ON K (x, 7, 2) udz = 
0 


— — yoou \ fe (z) —Q (z, n|""K (x, T, v) dz. (64) 


* By the Green's function method. 

































Here u does not depend on x, and can be taken out from under the integration sign. 
). Since cp = —1, according to the theorem, then —ycg > 0, and the maximum of the function II (r, u) relative 
to u is evidently attained for the condition A; =u =A, when 
1 
u(r) = BEM 4 Assign (| 9 (2) — Qe, 7) | R, 7, eae. (65) 
0 

If, instead of the function Q(x, T) in the equation (65), we substitute the expression given for it in (62), then 
we obtain an integral equation for the optimal control u= u(r). For example, if y = 2, A;=—1, Az = 1 (minimizing 
the square of the deviation), then 

1 T 
u(t) = sign || Q” (z) —{ K (x, T, *) u(t) dx] K(z, T, x) dz. (66) 
8 0 0 
Rearranging and changing the order of integration, we obtain 
T 
u(t) = sign [B (x) — \ L(x, 9) u (0) ds, (67) 
0 
where L (r, 9) is a symmetric kernel 
1 
L (x, 9) = \K (x, T, t) K (x, T, 9) dz, (68) 
0 
. 1 
B(x) = \ Q° (x) K (z, T, 1) dz. (69) 
0 
) This result coincides exactly with the result of solving the problem considered in [4], but it has been obtained 
more simply and rapidly. 

We now show how the result we have obtained can be applied in the solution of a known minimax problem, 
which can be formulated in the following way. Let the controlled coordinate x = x(t) be related to the controlling 
action u = u(t) by the integral relation 

T 
z(t) =@ (t) + \ K (t, *) u (x) dr, (70) 
0 
) where we will assume that 
9) #0 for OSig?, 
K (t, 1) #0 for OQ t<qtcT, (71) 
K (t,t) =0 for OQt*<tr<T. 
We also assume that 
4) lu@i<1 fr Ot. (72) 


We are to find a control u(t) (0 =t = T), satisfying the condition (72), such that the functional 
[° = J° (z, (t) = max |z (t) (73) 
0<t<T 


has its minimum value, whereupon 








x(T)= 0, (74) 


This last condition is equivalent to the condition 


T 
P =P (z(t))= \ 8(T — x) x(t) dt =0. (15) 


0 
The gradient of the functional I° is 
grad I° (x (t)) = 8 (tm — #), (76) 


where ty, is the time when the function |x(t)| (0 <t <= T) attains its maximum. 
The gradient of the functional I is 


grad J* (z (t)) = 8(T — 1). (77) 
Thus the function [I (r, u) takes the form 


T 
Il (x, u) = co | 3 (tm — t) [re (t)+K (t, t)u]dt + 


T 
+ a[ 8(7—2[F9 © + K (t, 0) wl dt = cog (tm) + 


+ co K (tm, t)u +a 9 (7) + cK (7, 1). (78) 


Since cg = —1, the maximum of the function I(r, u) is reached for 
u = u(t) = sign [K (tm, t) — ak (7, 1%)] (<t<7). (79) 


Since for T = Tr > ty we have 
K (tm, t) = 0. 


u = u(t) = — sign aX (7, t) for T >t > tm. (80) 


But the last formula for the control u(r) coincides with the formula for the control in the case when the time 
of the transition process is minimized. Thus, after the instant t = tr, when the function | x(t) | (0 =t < T) has at- 
tained its maximum value, the optimal control relative to rapidity of action does not lead to any worse value of func- 
tional 1° than the value | x(t,,) | in the optimal process minimizing this functional, We should remark that the 
time t,, can be determined by the method of successive approximations. 


APPENDIX I 
Proof of the theorem. We will assume that there exist an optimal control u = u(S) (SED, u€Q) and an optimal 
state Q = Q(S), such that, in view of the equation (10), the conditions (15) are satisfied: 


fn=0 (i=0, 1,..,p—4, pti,...4% 0<PCg 


and the functional IP attains its minimum value. 
We introduce the following functional: 
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imal 





q 
Te = >) Alt = M = Vit Be, (81) 
i=0 


where 
h= Who. Agls Mo =P do. ee Mg, AY DAs a Agde (82) 
I 1° o'+ (2) 
f=] .,h=l|.1,4—-0=]]. ; (83) 
19 r ©? (s) 


Further, we set up a correspondence between each point P€D and points S, €Nq (a = 1, 2,...,m), such 
that when the point P traverses the whole region D, then the point S,, traverses the whole set Ny (a = 1, 2,...,m). 


We fix an arbitrary point P’ €D and the points S'., € N, corresponding to it, which are regular for the controls 
uap(S,,), and surround the point S €Ng (a = 1, 2,...,m) with a region 5. C Ny, the volume of which tends to 
zero when the diameter d,, of the region tends to zero, 


We now determine the control i = u (S), giving the modification of the optimal control u = u(S): 


sey for S,EN, — 5, 


F 84 
>” ire at (a = 1,2,...,m), (84) 


u =u (Ss) = 
where v is an arbitrary point of the control region Q, 


The value of the functional I, is now calculated from the modified control (84), with an accuracy up to small 
quantities of a higher order than €, where € > 0 tends to zero when the diameter d,, of the region 5g © N, (a= 1, 
2, . + « » Mm) tends to zero: 


I, (u (S)) = 47 (Q (S) + 4Q(S)) + 
+10( 75,918) + 401, 719 as — 


='1(Q (8) +.4Q (9) + 1-0 (\ F(S, Q(S), u(s))as + 





+{ pee =D sowpes Heurwy'o7.0= 


— F(Ss’, Q(S"), u (S")))).- (85) 


Here the increment AQ(S), obtained from the modification of the optimal control u(S) in the equation (10), 
satisfies, with an accuracy up to small quantities of higher order than ¢, the nonhomogeneous Fredholm integral 
equation, linear in AQ(S): 


AQ(P) =e(K(P, S’, Q(S’), ) —K(P, S’, Q(S’), u(S')) + 





OK (P,S (86) 
se (35 5, QS 415) no eas 
As is known [5], the solution of this equation for AQ(P) is given by 
AQP) = e[KP, 5, QS, — KP, 8, Q, us) — 
4 \ a (P, 5) (K (S, S*, Q(S%, 9) —K (S, S, Q(S'), u (Sas, (87) 


where the matrix function M(P, R) satisfies the equation (22 
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Substituting the value of AQ(P) from (87) into expression (85), and expanding the functions @ and 1, in 
powers of €, we obtain 
Ie (@(S)) =4'T, (Q (5) + vo Fes, Q(S), u(S)as) + 
+ en (grad 1,(Q(P)), K (P, S’, Q(S"), ») — 
— K (S, S’, Q(S’), u (S’)) —\ M (P, S) (K (S, S*, Q(S’), ») — 


D 
— K (S, S"’, Q(S’), u (S’)dS) + 
am F (S, Q(S), u(S)) dS) 
+ eh" Oz = 





{ |S (P,Q (P),u(P)) 
5Q 


x [K (P, Ss’, Q (S’), v) — K (P, Ss’, Q (S"), u (S’)) _ 
—\ we, S)(K (S, S', Q(S"), ») — K (S, S$’, Q(S"), u(S")) dS] dS+ 


+ F(S', Q(S’), ) —F(S', Q(S’), u sy} (88) 
Since the functional I, takes the least possible value, then the main part of the increase in this functional for 

a variation in the control u = u(S) will always be non-negative, i.e., 
Al, = 1, (u(s)) — 1, (u (8) > 0. (89) 


When we replace I, (u(S)) in (89) by the expression in (88) and use (24) and (81), we find that the main part of 
the increase in the functional can be expressed as 


Al, =e [I (s’, v) —1(S’, u(s’))]>0 (90) 


where we have set a = A", b=A' and Ap = 1. 
Multiplying both sides of the inequality by —1, and assuming that ¢ > 0, we obtain 


II (S’, v) <1 (S’, u (S’)). (91) 


Since the inequality (91) is valid for any point v€Q, then the function I(S', u),defined by (24), attains a maxi 
mum with respect to the argument u for fixed S', i.e., for almost all S' €D 


Ti (S, u) = A (S), (92) 


which was to be proved. 


APPENDIX II 
Proof of the lemma, Let 


t 
» (Et, t) + v(t, t) =Snue 8) v (8, t) dO 


4 


as _fv@ fae t>t>t>h, 
vevag={) for >t >t>be, 


where U(r) is a square~integrable function on the interval [t, ty]. 
We consider the case when ty} >t =r =tg, Then the function p (t, r) must obviously satisfy the equation 
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y(t, t) + ¥() = Jet, 6) d0 ¥ (x). (98) 


t 


We now consider the case when ty = T> t = tg. In this case, we obtain 


t 
u(t, t) = \ (xp (t, 0) dO. (94) 


t 


For every fixed value of t, this is a Volterra integral equation. Because of the assumption in the lemma con- 
cerning the integrability of the square of /(r) in [to, ty], and the uniqueness of the solution of (94) [5], we have 


p(t, t) =O for hot > tp bt. 
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The maximum principle is applied in the investigation of a problem in the synthesis of optimal 
control, The optimization criterion is given by the integral of the sum of the squares of the de- 


viations. 


1, The Problem 
We consider optimal transient processes in a system of automatic control, described by the differential equations 





dz 
ae Seite isan ith a) 
j=1 


where the x; are the phase coordinates of the system and the a;; are constant coefficients. The function u describes 
the controlling action. Controls u are permissible that belong to the set U of piecewise-continuous functions satisfy- 
ing the condition 


|u(t)| <1. (2) 


There is, at the present time, a powerful mathematical method—the maximum principle of L. S. Pontryagin 
[1]—by means of which we can synthesize optimal controls, in particular optimal controls relative to rapidity of 
action, As is known, the problem of synthesizing controls for optimal rapidity of action can be solved, generally 
speaking, only by the use of high-speed computers. Even in cases where it is possible to find the analytical form 
of the function u(x;, . . . , Xp), the form of such functions is very complicated, It is therefore interesting to consider 
the problem of synthesizing systems that are rather simple, but at the same time optimal in a definite sense for a 
transient process. We consider the following problem in synthesizing an optimal control. For the system (1) it is 
desired to find, among all the permissible controls, an optimal control function u(x;, . . . , X,), such that the corre- 
sponding trajectory [x;(t), . . . , Xp(t)] of system (1), starting at any initial position [x,(0), ..., X,(0)], tends for 
t-» « to the origin, and such that the functional 


= \( So + cu* ) dt (3) 


converges and takes its smallest possible value. 


The functional (3) — the integral of the sum of the squares of the deviations — is an indirect criterion for op- 
timizing the quality of the transitional process. The coefficients aj and c are weighting constants. 


This type of problem in optimal control synthesis was investigated by A. M. Letov in [2], using the classical 
calculus of variations method, The difference is in the problem, Here the function u, describing the control action, 
is permitted to have discontinuities of the first kind, the presence of which will not permit the use of the classical 
calculus of variations method. The problem of minimizing the functional (3) for the system (1) was also investigated 
by A. re Hee ae eeeeeee OF Cynnantc poeguemming. In the present article, the synthesis problem will 
be solved by the maximum principle. 
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2. Solution of the Problem 








It turns out that the problem posed above can easily be solved by the maximum principle. We first of all 
determine the optimal control in the form of a function of the time, The method of solution, using the maximum 
principle [1], is as follows. 


We introduce the notation 


fo (24, . . +) Zn, U) = 3} aja + eu 


j=1 


fa (Spy 2 oy Zn Uw) = Dd} 2425 + du (i=1,...,m). 
j=1 


In addition to the phase coordinates x, . . . , X, for the system (1), we introduce one extra coordinate Xp, 
satisfying the equation 


ar er 


In the (n + 1)-dimensional phase space X, there is a vector X= { Xp, X,...» Xp, } with rate of change given by 


op = 1 (tas + + ++ Tne Y)s (4) 


where f(xy, . . . , Xp, u) is a vector with coordinates f(xy, ...,. Xp, U),.. +» fylX%q,. + - » » Xp» U)- 


We now form the system of differential equations adjoint to the system (4), with the new variables » = {p»o, 
¥y--- Ont 


dy) a 
— 2 in ye = 0, 


T--3e Ya = — 2agory — >} ays (é=1,...,m). (5) 


j=1 


We also form the Hamiltonian, which is the scalar product of the vector » and the velocity vector dX/dt. We 
denote the Hamiltonian by 


H (z,, u) = Deg t= Dh (24, » + «9 Ruy B) = 


_ ¥(3 uxt + cu*) + >? ¥, (> Qyj25 + bu). (6) 


j=1 


For fixed xy, ..., Xp and Po, Oa, - - - + » Wp» the function H becomes a function of u; we denote the upper 
bound of the values of this function by 


= (z,) = | H (z, », u). 


According to the theorem giving the maximum principle [1], the condition necessary for the optimization of 
the control u(t) is as follows. Let u(t) be the optimal control and let x(t) = { xg(t), x4(t), . . . » X(t) } be the corre- 
sponding trajectory of system (1). Then there exists a non-zero, continuous, vector function »(t)= { polt), o2(t), 
++» ®p(t)}, such that the functions x(t), p(t), and u(t) satisfy the systems (4) and (5), while at any instant of 
time the Hamiltonian H is maximum relative to u, i.e., 
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H (z(t), p (), u ()) = M Iz, (t), 9). (7) 
If x(t), p(t), and u(t) satisfy these conditions, then 


tpo (t) = const <0, M [zx (t), p(d)) = 0. (8) 


It is easy to show, from (7) and (8), that the optimal control u(t) is given by 


f n n 
DM soe [7 DH] <1, 
=1 =1 
u(t)={ 1 for z Dom > 1, (9) 
—1 fe 5 Yoy,<—1. 





i=1 


We note that in the right-hand side of (6), the terms containing u in the expression for H are 


crpo (t)u® (t) + u (t) >) b,», (2). (10) 


i=1 


The expression (10) evidently has a maximum at the same time as the function H for fixed po, .... ¥n: 
Xq, . +» Xp- Therefore the necessary condition (7) becomes 


cae (t) wt (0) + (®) 5 6,9, ) = sup [ey (fw? + wD B,¥, (0). any 
i==1 


i=1 


We rewrite the expression (11) in the form 


cro (t)u® + u >) b,, (t) = cm, (2) [u + 3% 2 bp, (| at ie yy “a of. 


t=] 


The function H (x, », u) in (6) is homogeneous relative to p, (a = 0,1,...,m). Therefore, according to the 
condition (8), we may assume that p, (t)= —1. It is then easy to see that if u(t) is determined according to (9), the 
expression (11) attains its maximum, 


Now, having the expression for u(t) in (9), we can establish the relation petween the functions » ;(t),....¥n(t) 
and the functions x;(t), . . . , Xp(t), and this directly determines the optimal control u in terms of the coordinates 
Xyp +++» Xpy- We first of all assume that after the time of the transition process, the control u has not reached its 
limiting value. Then, substituting the first expression in (9) in the right-hand sides of the system of equations (5), 
we obtain the linear system of 2n differential equations 


Mm bb < ; 
Fm Dowty Dey Gm 2a — Bay Cte an 


d=1 


in the 2n unknown functions x;(t), . . . . Xp (t), a(t), . - +» y(t). The system (12) can be solved for given values 
of the boundary conditions. For given, fixed initial conditions x,(0), ..., X,(0), we must choose the corresponding 
initial values 9 3(0), . . « » ®n(0), so that the solution x;(t), . . . , Xp(t) of the system (12) tends to zero for t-> «. 
The characteristic equation of the system (12) is 
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11) 


(12) 


lues 
nding 








st 2am? b 6 2 £8. 643 eee 











Using elementary methods of transformation of the determinant (13), we can show that the relation A (A) = 
= A(-A)= 0 holds, i,e., the equation has roots in pairs, equal in magnitude but with opposite signs. We will assume 
that these roots are simple, and will denote them by Ay, .... Ani ~Age. +++ Ape Let the condition RedAy < 0 
(i=1,...,M) be satisfied, which is evidently always possible. Then the general solution of the system (12) can be 
written in the form 


n 


Aat 
Xj (t) = > cA; (Aj) e ’ |- 
j=1 


1 —A;l 
CnpjA(— Abe  , (14) 
1 


M: 


~~. 
if 


YY A jt 
Pi (t) = Dj cjAngi (Aj) + 


j=1 j 


—A;t 
Cn4jAn+i (— Aj) e ’ 


Ms: 


Ly 


(15) 
(i= {,...,a) 


where the A, (Aj) (k=1,2,...,2n;j=1,...,m) are (n~1)th order minors of the rows of the determinant 
4 (Aj) for which they are not all zero [3]. 


The arbitrary constants cy, ... , Cgp are obtained from the initial values x,(0), . . . . Xn(0); p3(0), . . - » Bp (0). 
In order that for t +» « the functions x;(t), . . . , Xn(t) tend to zero, it is necessary that cn 4j=O(j=1,...,n). It 
can be seen from (15) that the cp 4; are zero if 


%, (0) ass p> Cj)An+i (A;) (i =1i,..., n) (16) 


j=1 


(we assume that | A,+; (— A;)| +0). 
On the other hand, the arbitrary constants cj (j= 1,... , m) must satisfy the conditions 


x; (0) = >) cjAy (Aj) G=4.---. a). (17) 


j=1 


When we eliminate cj from (16) and (17), we obtain the relation between x; (0) and (0): 
°p, (0) = >) 7,25) (0) = 4,.--.m). (18) 
j=t 


It should be remarked that this relation must hold at some instant of time, and we can always take this instant 
for the initial time. Thus the optimal control u inside the boundary has been obtained in the form 
n 
{ : 
u=5- 2 km, (19) 
i=1 


u 
k= > bi¥ 5: 


j=1 











It is obvious that the relation (19) only holds in the case when the absolute value of the sum = (ky +. 
+ + « + KpXp) ds not greater than one. We can easily see that when |——(kyxy +... + kan) | > 1, then the optimal 


control must be beyond the boundary(i.e., lul = 1)if the right-hand boundary conditions are fulfilled, that is, the solu- 
tion of system (1), which is nonlinear due to the nonlinear relation between_u and x), . . ., Xp, tends to zero, In other 
words, solution of this system must be asymptotically stable. Thus the optimal control is given by 


A Spas fo l= Dam <1, 
i=1 i= 


u={ 1 fr £>kn>1, (20) 


i=i 


—1 for Yan <—1. 


. i={ 





The question of the stability of the system (1) with the control (20) has to be investigated separately. We see 
that, because of the continuity, stability is guaranteed in the region Lote +. . « + KpXp)] = 1+ ©, where ¢ is 
sufficiently small, The question of the stability of the system (1) with the control (20) thus reduces to the problem 
of finding a region of attraction, the boundaries of which lie beyond the limits Lars +ooet kXn)I =1. This 
problem will be the object of a separate investigation. 


3. Comments 

It is interesting to compare the result we have obtained with the result found in [2] by the classical method 
of calculus of variations, for a more restricted class of control functions than we have used in the present article. 
We find that the results are the same. The system (12) coincides with the Euler system of equations (3.3) in [2]; 
in the former, the functions y(t) play the part of Lagrange multipliers, and the results are obtained by solving these 
systems of differential equations. Thus no new results are obtained for the problem by generalizing the class of 


permissible functions. 
The author wishes to thank A, M. Letov for several valuable suggestions which were of use in the present article, 








4. Example 
We consider the system of the second order 
oF + 2a + ws =u (Jul <3). @1) 
The function to be minimized is 
co 
{ (az* + curat. 
0 
When we transform (21) to the normal form, we obtain the system in the form 
a >. Zaz, + wa, 
(22) 


1 
Fm «o's, — 2asy + 5s > _ —), + 2ahs. 


The roots of the characteristic equation for the system (22) are 


= —oa+ iB, a= —a— if, s=—a+ f, me =a — ff, 





where i= 7-1, p= Vu*=a, 








lu- 


ticle. 


11) 


22) 














The general solution of the system (22) has the form [3] 


Hy = Cy%q, + Coty. + carga + cans, 
Ze = Cy2y + coten + cores + core, 


(23) 


Di = CrtPay + catia + cathis + chic, 
Pa = cythar + cothan + cathes + catPac- 





Here 
21, = e~*! (a cos Bt + B sin Be), ay, = e~*' (a sin Bt — B cos Bt), 
43 = e*' cos Bt, 24 = =e sinBt, 
22, = e * ((6* — a?) cos pt — 208 sin 8), 22 = gs (2a8 cos Bt + (6* — a*) sin Bd), 
aay = e*! (+ cos pi — £ sinB t) ' rm = &'(F cos pt + = sin Br) ; 
— 3a 3a 
vu = *( 5) cos pt — 5 sim Be) vin = «(SE cos pe — F sin Br), 
is = 4aw%e™ cos pt, Pre = 4a w%e™ sin Be, 
>a = _ e~*' cos Be, Ya: = a e~* sin Bt, 
pos = e*! (4a cos Bt — 4afs sin Be), ‘bas = e*' (4a cos Bt + 4a* sin Bt). 


In order that x, and x tend to zero for t + «, it is necessary to choose the initial values » (0) and »2(0) so that 


3a 
p1(0) = errr (0) + catia (= — Far + So. 2 (0) = cya; (0) + capes (0)= — x Cy. (24) 


Actually, since the determinant 





Dra (0) Pre (0) | ie reg 0 
Pas (0) Hoe (0) 4a* 4aB 


is not zero, it follows from the last two equations of the system (23) that oc, = c,= 0, 
In this case, the equalities 
xy (0) = ¢y273 (0) + cary (0) = ac, — Bea, 
5 (0) = exems (0) + csvm(0) = (B* — ate, + 2afier she 
will hold, 
If we eliminate c, and co, from (24) and (25), we obtain the relations between » (0), p2(0), and x,(0), x_(0): 


— | _ i 
¥1 0) = Se [Gat + 6) 521 (0) + 220], HO = Sele O +z 2x). 


Hence the optimal control is obtained in the form 


f 
S1, 





55 21+ 57) for | —sea(e: + oes) 


u (x14, 22) = | i for —5o5(a1 + 2s) >i, 


me for — sas (21 + ges) ack: 
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THE CONNECTION BETWEEN S. A, CHAPLYGIN'S THEOREM 
AND THE THEORY OF OPTIMAL PROCESSES 


A. I. Averbukh 


(Moscow) 

Translated from Avtomatika i Telemekhanika, Vol, 22, No. 10, 
pp. 1309-1313, October, 1961 

Original article submitted April 11, 1961 


The connection between S. A. Chaplygin’s theorem on differential inequalities and the theory of 
optimal control is considered for the purpose of obtaining a necessary condition for the applica- 
bility of this theorem in a given interval, and also in order to obtain a statement concerning the 
dependence of the magnitude of this interval on the character of the variables in the problem, 


If two functions y(x) and 2(x), which coincide together with their first n—1 derivatives at x = Xo, satisfy the 
differential equation y ) = f{x, y, y’,..., y{@-1)) and the inequality z() = f(x, z, 2", « «« » 271) respectively, 
then clearly on some interval [X9, X;] they satisfy z(x) = y(x). The theorem of S, A. Chaplygin on differential in- 
equalities consists in indicating an upper limit for X, for linear equations. It is as follows: on the segment [X», X;]) 
the solution of the adjoint equation with zero initial conditions on the derivatives and positive initial condition at 
X = X; on the unknown function itself should be nonnegative [1, 5, 6]. 


In the case of nonlinear equations the matter is significantly more complicated. B, N, Petrov [2] constructed 
an example of a nonlinear equation such that for any ¢ > 0 there exists a comparison function z(x) for which X; < ¢, 
since for X = X, already z(x) < y(x); so that the upper limit of X, is equal to Xp». 


The problem® of finding conditions under which 2(x) = y(x) on the segment [Xo, X;], if 
y™ =f(z,...,.y¥*™), 2 ™>fi(z,..., 2), 
may be solved with the aid of methods used in the theory of optimal systems [3], 
In fact, if, on the segment [X», X,] the conditions 


y; = fi(Z, Yr «+ +s Yn)> 
z.> fy (Ze Bay + ~ +o Zn)s (i= 1,5... 2) (1) 
yi (Xo) = 2% (Xo) = Yo 


imply that 
BDPVYr (i =1,. <4 M)s (2) 


then for all u;(x) = 0 the solutions of the system 


y, = hi (z, Yi, + + +5 Yn)» 
zi, = fi (Z, Yay + « +s Zn) + (2), (3) 
y; (Xo) = 2% (Xo) = y? 


satisfy the inequality 


* In the following it will be termed Chaplygin's problem. 
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% (z) < 2 (z) (i mm i,.- +s n) (4) 
for all x€[X», X,}. 


To determine the conditions under which (4) is fulfilled, we use the fact that z;(x) is a functional of u(x) = 
= { us(x), . « « » Un(x) }, and y;(x) is the value of this functional for u(x) = 0. Consequently, it is necessary to deter- 
mine the variation in values of the functional when the control is varied from u = 0, 


In [4] a formula is obtained giving the increment of the value of the functional 


S = Xe,2, (t), 
where the functions x;(t) satisfy the system 


Wy = fy (Try + + oy Duy Ways + oy Une A) (id, oy ms HO) e,)- (5) 
If the control u(t) is changed from vj (t) to v(t) = v(t) + Su(t), then 


n T 
85. = Deda, (7) = —\ 1H (2, p, o, () — H (2, p, ¥, dt — a, ©) 
n=] 0 
| n 
H (z, p,u,t)= > Pil (z, u, t), (7) 


i=] 


and x4(t) and p;(t) are solutions of the canonical system 


0H (xz, p, v®, t) 
in RAEN SAGE, wah Fi (mj 6, 


(0) (8) 
p= — e- tala qu asi 8p \ico ays 





qi 





The quantity n (t) represents a residual term for which the estimate 


{In|<A 


ae 


>) dutdt 
1 


is deduced in [4]; here A is some constant which does not depend on the quantities 5u,,, and r is a segment of the t 
axis on which the quantities 5 u;,, are different from zero. 


The problem (3) considered is written with the help of formulas of the following form: 


vi (z)=0, vi) (z) =u (z), 
Sz, (2) = 2u (2) — yx (2), 
Cy = Lande; = 0 for j + k. 


In this cxse 


Xn 
tx (X) —yx(X) = —| Dipmdz — 0, (®) 
x 1 


y’ as he (y1, ++ +9 Un» z), 





n 
pao Ofe (Ya -- ++ Yn» 2), (10) 
P,; uP, ay; a (i arr n), 
Px (X) = — Cy, Pj (X) = 0 for j#k (A =i,..., n). (11) 
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The solution of Chaplygin's problem reduces to finding conditions under which for any u;(x) = 0 the right side 
of (9) is nonnegative. Carrying through exactly the reasoning in [4] which leads to the proof of the maximum prin- 
ciple, and considering the estimate for 1, we obtain a necessary condition for the nonnegativity of z(x)—y)(x) for 
all x€[Xo, X;] and u,(x) = 0 as follows: the functions p;(x) which satisfy the system (10) should be nonpositive on 
the segment [X, X;] for all systems of initial conditions (11), 


Thus the following is a necessary condition for the applicability of Chaplygin's theorem on the segment [X», Xq] 
for the system (1): for any k and any X €[X9, X,] the quantities pj, (x), which are solutions of n systems 


y; = fi (z, Yis- + +s Yn)s 


» 


’ of, 
Pin = — Li Pox dy, Dnt in ae (12) 


v=] 


Py, p(X) = — C4, nds, e (04,4 >0) (R= 1,.. 6m), 


should be nonpositive on [X», X). 
Here 5; ,is the Kronecker symbol (5; = 1 for i= k and 5; , = 0 for i #k). 


We note that the equations for Pi,k form a system adjoint to the system of variational equations associated with 
the original system (1). 


From the arbitrariness of | c, , | and the linearity of the system (12) it follows that for the applicability of 
Chaplygin's theorem on the segment [Xp, X,] it is necessary that the solutions p;(x) of the system adjoint to the varia- 
tional system possess the property that for any X €[X», X;], pj (x) < Ofor x€[Xq,X], if py(X) <0 for alli=1,...,0, 


In particular for n = 2 it follows from this condition that of the two functions p;(x), p2(x) which are solutions 
of system (12), only one may change sign (for any initial conditions), and then only once. 


In the case of a system linear in yj, it is shown in [4] that n = 0 and the condition mentioned is also sufficient. 


If it is required that the condition yj = zjbe satisfied not for all i, but only for some of them, then an analog- 
ous condition is obtained, but imposed only on some of the variables, In particular for the case of a single equation 
of nth order it is necessary that the solution of the adjoint equation with zero initial conditions on the derivatives at 
x= X not change sign on the segment [X», XJ for any X €[X», Xy]. These conditions coincide with the conditions in 
(1, 5, 6], In the case of a linear equation of nth order this condition is also sufficient. 


In a completely analogous manner one may investigate the following generalized Chaplygin's problem ("ex- 
tension of Chaplygin's theorem"): for which conditions do the inequalities 


uy, (z) > v, (z) > 0 (13) 
imply for any uj(x), v(x) the inequality 
%(z) > (2) (14) 
for solutions of the system 
Yj, = fe Yas - + +s Unt) My 2” = fy (By, + 0, Say 2) Oy (EA, 24, 0) (15) 


with initial conditions 24(X») = yg(Xo) = y$- 
A necessary condition for the satisfaction of (14) on the segment [X», X] with (15) holding is the nonpositivity 
on (Xp, X ] for any X €[Xo, X;] and any k of the quantities p, ).(x), solutions of the n systems of equations 


y= fi Yay + + +s Uns x) + uy (z), 
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Pp 3 p oe (i=1 k=1 ) 
= - — t= 1,..+,%,R =1,..., A), 
k et 


(16) 
Pr p(X) = — Ce Sin (Cr, e > 9), ys (X,) = ¥? 


for any 
u, (z) > 0. 


In this case the applicability of Chaplygin's theorem depends on the behavior of the collection of all solutions 
of the system (16) obtained for all admissible values of uj(x), whereas in the usual Chaplygin's problem only the be- 
havior of the solutions for us(x) = 0 is important. This fact is noted in [7], 


With the help of the method indicated one may investigate the question of the dependence of the interval of 
applicability of Chaplygin's theorem on the variables of the problem, In other words, one may solve the problem of 
the possibility of expanding this interval by means of a suitable change of variables y;, .. . , yn. 


We introduce new variables 
Ynys = Pi Yr -+ +r Yn 2%) (i= 1,..., 2). 
Then instead of (12) we obtain n systems 
y¥,= h (Yi, oe 09 Yn 2) oe Uj, 
= + a) 4 
= 2m, % 
Paste = O (i=1,..., a; Rai,..., nr), 
n af, n 
P. --{2 Pik Oy i+ Di Pats, w X 


x | Sage et uy) + pee + zal 


(17) 


with initial conditions 


Py (X) = 9, Putin (X) = — 84,2, 
yi (X_) = ¥?, Ynit (X,) = 9: (y?,---» ¥%, XQ). 


The conditions for the applicability of Chaplygin's theorem will be the nonpositivity on [X», X) for all 
X €[Xq, Xq] and all k of the quantities 


r= Puet 3% By, Pot, ks (18) 


formed from pj, Pp +4 and y; and satisfying system (17). This follows from the expression similar to (9) for the system 
(17), For the case of Chaplygin's theorem in the usual sense it is necessary to take uy = 0 in (17), and for the extended 
theorem of Chaplygin, to consider all uj = 0. 


From (17) it follows that the quantities w; ; satisfy n systems of differential equations: 


of, 
Ww, = — Wie ar (i=i....,a;k=1,..., n) (19) 
a Oy; 


with initial conditions 


wi,x(X) = —>- (19a) 
(since pn +i, k =~ 54,4) 
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Comparing (19) with (12), one sees that the equations for w; , ate the same equations adjoint to the system of 
variational equations, as are the equations for p; , in (12), From this it is clear that if the functions ; are such 
that 39; / dy), do not change sign for any i and k, then it is not possible to expand the region of applicability of Cha- 
plygin’s theorem by transforming variables, The case when 9 9/0 y, may vanish for certain y is more complicated 
and is not considered here. ff 


In conclusion the author expresses deep thanks to L, I. Rozonoér for his interest, great help, and valuable coun- 
sel as regards making the bases of this article sounder. 
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POSSIBILITIES FOR SUPPRESSION OF DISTURBANCES 
IN A CERTAIN CLASS OF DYNAMIC SYSTEMS 


M. V. Meerov 


(Moscow) 

Translated from Avtomatika i Telemekhanika, Vol, 22, No, 10, 
pp. 1314-1323, October, 1961 

Original article submitted March 29, 1961 


The connection between the structure of a dynamic system and the intensity of disturbance 
effects is discussed as depending on the points of entrance of the disturbances into the system. 
It is shown that for a class of structures which remain stable with unlimited increase of am- 
plification, disturbances can be suppressed with any degree of accuracy regardless of the point 
where they have first been applied to the system —except at the main input. A procedure for 
carrying out the design of structures incorporating this property is given. The problem has 
been solved for systems of multiarticulated controls, as well as for systems with a single con- 
trolled function. 


The most important problem in the contemporary theory and practice of automatic control is choosing the design 
and parameters of the systems in such a way that the effects of disturbances will be at a minimum. A good deal of 
attention has been and continues to be given to this problem in modern literature [1-7]. In these papers, basically 
the case is investigated when a signal together with some disturbances is applied to the input of the system, and the 
problem consists in separating out the useful signal from the background of disturbances. In many automatic control 
systems, and in a number of tracking systems the useful signal is not accompanied by any disturbance at the input. 

In these cases, however, disturbances are applied at various intermediate points of the control circuit. They consist 
of irrelevant loads and of perturbations of various kinds, some of which may be also of a random nature. 


In this paper just such a case is discussed, It is shown that for one class of structures (namely those which ad- 
mit limitless increase of amplification without harm to their stability), under certain conditions specified below, sup- 
pression of disturbances consisting of irrelevant loads or of perturbations can be carried out by increase in amplifica- 
tion of several components of the control circuit. It is made clear how the effect of a disturbance depends on the 
point of its application, as well as on the structure of the entire control chain, The problem is solved for systems 
having a single controlled quantity or function,and also for systems having several interrelated controlled elements 
(the latter being sometimes called multiarticulate systems) [11, 12]. In the case considered in this paper the dis- 
turbance can be applied at any place except at the input of the signal, 


1, System with a Single Quantity to be Controlled 

In Fig. 1a block diagram of an automatic control system of this kind is presented. A useful signal x,, (standard 
test signal) is applied to the system, and the latter is required to reproduce it. It is assumed that xg; is free of any 
disturbances. To simplify the following discussion, we shall assume that our system consists of four dynamic com- 
ponents, This does not make the conclusions which will be reached below any less valid in the general case of many 
components: Indeed a derivation which is valid in the general case has been given in paper [13], 


The transfer function of the system component i not subjected to the effect of a disturbance is designated by 
K,Ry(p)/D;4(p). With respect to components upon which disturbances are acting, it is assumed for the sake of generality 
that the points of application of the disturbances do not coincide with the points of application of the useful signal, 
and the transfer functions for the disturbances are designated by K' ,R';(p)/D';(p). 


A condition which must be fulfilled without fail by all systems of the class under consideration is that the 
first component, that is,the one to which the useful signal is applied, must not be subjected to the action of any dis- 
turbances. It is assumed that the amplification factors of the components which are not directly subjected to the 
action of disturbances can vary within wide limits and assume sufficiently large values. 
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Fig. 1 


Let us find the transfer function for the complete system shown in Fig. 1. We have 





%,= ae (zst — Zour], 
v= Bi 4 SO, " 
Pout = “pote as + Se ie 


After exclusion of x3, %, and X, from equations (1) to (4), and a few transformations, we obtain the following 
expression for Xoyt: 


4 
[] 2: (2) Dy (P) Di (P) Krk aKaKery, 
i= 





Tout = 4 4 bie + 
[] 2s) Dy (>) D, (p) + Dy () Dy (P) [] Ry) BiksKK, 
i==] i=1 


4 
D (p) Ds (p) Dy (p) Ro (Pp) Ra (p) Ry, (p) KyKsKati + [] Di(p) Dy (p) Ry (P) Kyte 





+ ; ~~ (5) 
[] 2: () Dye) D, (p) + Dj) Dy (p) T] Re) KiKsKoKs 
i==1 i=] 


We divide the numerator and denominator of (5) by KyKs. Designating 1/K,= 1/K, = m, we then obtain the 
following expression for Xgy¢: 


‘4 4 
| m: I] Di (P) Ds (p) Dap) + Da (p) D4 (p) TL RiP) KsK,| Tout = 
i=1 


i= 


a 
= D;(p) Di (p) [] Ri (p) K2K ater + 


i=] 


+ mD, (p) De(p) Dy (p) Rs (p)Ra(p) Ro (p) KoKaht + 
+ m® |] D, (p) D3 (p) Ri (p) Kalo. 


t=1 


(6) 
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From here on we shall assume that the disturbance and its derivatives are limited in modulus; otherwise, they 
may be any functions of time, including random functions. 


Let us consider K, and Kg to be sufficiently large figures—theoretically infinitely large. If the system remains 
stable when m -~ 0, then we obtain from (6) 


lim Xout = *st° 


Putting it differently, when the system remains stable with unlimited increase of amplification in the com- 
ponents which are not subjected to direct action of the disturbances, then it is possible to obtain as exact a reproduc- 
tion of the useful signal as desired. Formula (6) shows also that disturbances are the more strongly suppressed, the 
farther the components to which they are applied are from the beginning of the system. The basic principle of this 
system structure is that the disturbances are suppressed by increased amplification in the components which are located 
in the control circuit between the useful signal input and the component which is acted upon by the disturbances. 


2. Physical Feasibility of Class of Structures Under Consideration 

As was already stated elsewhere [13], and as follows from the results obtained in the preceding section of this 
paper, the system considered here can be carried out in reality, provided it belongs to the class of systems which re- 
main stable when an unlimited increase of the amplification factors takes place. This question requires special con- 
sideration in this case: Stability or instability will depend upon the location of the poles of equation (5) or, which 
amounts to the same thing, upon the roots of the characteristic equation 





4 4 
(o'T] D.(p) +1] iter KK] D5 (p) Di(p)=0 at m->0. nm 
t=1 i=1 


In order that all roots of (7) shall be on the left from the imaginary axis on the plane of roots, it is necessary 
and sufficient that the roots of each of the factors in (7) be on the left from the imaginary axis on the roots plane at 
m-»+ 0, But since the roots of the equation D',(p)D’, (p) = 0 do not depend upon m, these expressions can be left out 
from consideration. And we will assume that they are on the left of the imaginary axis. Thus, the stability will be 
determined by the roots of the following equation: 


mT Di(p) + KK. II Ri (p) = at m-— +0. (8) 


i=1 
It is known [10] that the roots of equation (8) will be located on the left from the imaginary axis, provided the 
following conditions are satisfied. 
Let us designate the degree of the polynomial II D;(p) by Ne, and that of the polynomial II R; (p) 
by Ny. Then the conditions of stability will be as follows: 


a) Ne-N, =2; 
b) the equation Il Ri(p)K,K,=0 must satisfy the stability conditions; 
t= 
c) depending upon whether Na~Ny is equal two or one, a certain definite relationship must be valid between 


the coefficients of the polynomials Il D,(p) and K,K, Il R; (p). 
i=l i=1 
Let us discuss the cases which are most difficult from the point of view of their being carried out in practice. 
Suppose that all components in the circuit of the system shown diagrammatically in Fig. 1 are elementary dynamic 
sections. For the sake of simplicity we shall consider the circuit of the system represented in Fig. 1 to be stable by 
design [14], and that all its components are of an aperiodic character. 


If there is a total of N such components in the circuit, and v of them have large amplification factors, then 
equation (8) can be rewritten as follows: 





he 





N n N 
mT] d+Tp)+][ x [] xi=0. 
i=1 i= i=n+o+1 


If we introduce the designation 


T[K: I] Ki=Kaeg: 


i=1 j=n+o41 


then we obtain the following equation: 


N 
m° || (14+ Tip)+ Kieg = 0. (9) 


t=1 


It is obvious that equation (9) will satisfy the stability conditions only when N =2. This is a trivial case which 
is of little interest. Let us investigate the possibility of stabilizing the system when N > 2. 


From paper [10] it is known that to stabilize a system of the type represented by formula (9) it is necessary to 
introduce the derivatives at least from the (N-2)th to the first. Let us introduce into the system N~2 amplifiers 
whose amplification factors can be made sufficiently large. Let us correct these amplifiers by negative feedback 
paths having the transfer functions p; / (1 + Typ) (Fig. 2). Concerning the rest of the system we shall assume that 
among the N aperiodic components, there are « which are not directly subjected to the influence of disturbances, 
and that the amplification factors of these can be varied within wide limits. To simplify the computing operations, 
let us arrange the components which are directly acted upon by disturbances in two groups as shown in Fig. 2. 


K; Kk 
Ko = ky, = Dip) (Pp) 


% >h-A >a > FH> 





























Fig. 2 


For the first N-2 amplifiers — which use feedback — we have the following formulas: 


N—2 


Ki g(t + Tp) 
Ki 

"= Dey“ (11) 
(12) 

_ Ks K, 

= D (p) aT rw hy» 

_ Ks 13 
sa a ow 

a K, 
Z out D, (p) Ts + D, (p) he. (14) 
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Excluding X;, %, %, and Xj. from equations (10) to (14) we obtain 
N-2 


(Tow TI Tl (Tw + 1 + Kiss) Dy (p) Di (p) + 


4 ay 
+ TK: [1 kit + 79) D5(p) D.(p)} tout = 
toms] t==} 
4 N-—2 


= []K: Tl Ktol4 + Tw) Ds (p) Di (p) me + 
i=1 t=] 
N-2 


+D, (p) Dz (p)Di(p) T] (Tip +4 +Kiow) KeKoKsh + 


f=] 


‘ Na 
+ [] Di (Ds (>) [] (Tp + 1 +K joy) fr (15) 
t=} t=} 
N-2 
We divide both the right and left side of equation (15) by Il KigK,Ke, and assume that the quantities 
1 
Kio , Ky, and Kg are all of the same order of magnitude. Then, designating a + m,after a few ele- 
oO 1 
mentary operations we obtain the following expression: 
‘ N—2 
{m* 112.00) [m+ TT It + Pup) + mew 2 IT + Pur) + 
t—1 j= - 
N-2 | 
+m*—y DY [] (t+ Te) +... +e] 4 
jaa i=1 
tej 
i*j+1 
N-2 N—2 


+ KK, Ta +7) D, (p) Dy (P) tour = K2Ks [] (4 + Tip)D3 (p) Da p) 2+ 


i=] 


N—2 N—2 
+ m| mes T] (1 + 1) + mn I] @+7)+.. + ur |x 


j=1 4 


N—2 


x KgK,D, (p) Dy (p) Di (p) fy + me [mr I] ¢ + 79) + map 


N—-2 4 
xb G+ Te) +. +HN] TI Ds (2) Ds (p) fr: (16) 


i¢j 


It is evident from equation (16) that the magnitude of the amplification factors of the component amplifiers 
using feedback has no effect upon reproduction of the disturbances: the disturbances will neither be more amplified 
nor more suppressed by changes of the amplification factors of these components; the only effect of these com- 
ponents upon the disturbances consists in changing their amplitudes by a factor of wN-2, When u < 1, the distur- 
bances are correspondingly decreased. Substantially, however, like before, the disturbances are suppressed by a 
proper choice of amplification factors in the components which are not directly subjected to the action of the dis- 
turbances and are not provided with feedbacks. 


Let us now examine the left side of equation (16) which when equated to zero becomes the characteristic equa- 
tion of the system and determines the stability of the system, and hence the feasibility of itsrealization. We thus have 
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N—2 N—2 N-2 


m® TT (i+ Tw) I] Di(p)+m*—p 5) [] A+ rll] D,(p) + 


t=] f—1 f= i=] i= 
iej 


N—2 
feet mby®—* TT De(p) + KeKe [] t+ Tw) = 0. 


fom] i=] 


(17) 


The difference in degree of two consecutive polynomials in (17), except for the two last, ones, is one, and the 
difference in degree of the last two polynomials is two, since the degree of the polynomial Il D,(p) is equal 4 


according to the conditions of the problem setup (because the case of N aperiodic circuit sections is considered). 
Thus, the structural conditions of stability are satisfied [10], The degeneracy equation in the case under considera- 
tion is 


N—2 


Ky [] (1 +7) =0, (18) 


i=1 


which always satisfies stability conditions. Consequently, according to paper [10], for stability (feasibility) of the 
system it is necessary and sufficient that an auxiliary equation of the third kind [10] satisfy the stability conditions. 
And this can be secured by proper choices of values for the parameters T; and for the amplification factors K,K,. 


Thus the possibility of realization of the described structures has been proven. In the case just considered 
there were additionally N—2 amplifiers with large amplification factors introduced into the system. If N-1 am- 
plifiers are introduced, instead of N—2 amplifiers, then an auxiliary equation of the first kind [10] will be obtained, 
which in the case under consideration satisfies the stability requirements always, because it is reduced to the form [10] 


[] (4 + 719) = 0. (19) 
i=1 


The number of amplifiers can be cut down to N/2, Then, the auxiliary equation will be of the second kind, 
and the feedback parameters (which in this case will be links of the second order) must be so chosen that the auxiliary 
equation will be satisfied. 


3. System with Several Quantities to Be Regulated 








Let us now examine a system of multiarticulated (interlocked parameters) control: it is assumed that in the 
controlled object there is not one quantity to be regulated, but several quantities all connected with each other 
through their serving to control the same object. Let us 
examine the case when each disturbance acts directly upon 
the object being controlled, It will be possible to extend 
the results of examination of this case to any other case with- 


out great difficulty, 


In Fig. 3 we have presented the structure of a system for 
the first of the several regulated quantities in the case when 


there is a separate perturbation acting on each of the regulated 
quantities. 


Let us write down the transfer function for the ith regu- 
lated quantity using the symbols and designations shown in Fig. 





KigRio(P)  —_Kighig (P) 4 KioMio (P) KigRig (P) 
% out = ~ Dy (P) oP) y¥i— ~D,, (P) > WP) > Ain (p) Ly + Di ) + fin + Mae (p ) is (20) 


ofl 
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E = Zist — Ziout, (21) 


Ky Rip (Pp) ie KipRip (Pp) 


Fas Dip (P) Fist “Dw Ot: (22) 


Excluding y; and € from equations (20)-(22), we obtain, after the necessary transformations have been carried 
out, the following expression: 


{Dio (p) Dip (p) + hip Rio (p)] Dio (Pp) Zout + Dio (Pp) Dip (P) KioRio (p) 


x >) Hiv (P) e= KipKioRio (p) Dio (p) tt + Dio (p) Dip (p) KiRio (Pp) fi + 


vo=1 
v¢i 


+ KioRio (p) Dio (P) fin: (23) 


In the above equations the subscript "o" is used to distinguish those members which pertain to the controlled 
object; the letter p (= Russian r) when used as a subscript designates that the quantity to which it is a subscript per- 
tains to a regulated quantity; fj is a disturbance applied to the ith regulated quantity; fj;;—the load in the ith regu- 
lated quantity; n—the total number of regulated quantities; and ajy(p)— an operator which determines the connec- 
tion between the ith and the vth regulated systems. 


The complete system of all equations for the multiarticulate system of regulation is obtained by giving i in 
equation (23) all values from i= 1 toi=n. 


We divide equation (23) by Kip and substitute m for 1/Kjp. Then we obtain the following: 


(mDio(p) Dip (p) + KicRio (p)) Dio (P) Zour+ 
+ mDigDip(p) KioRio(p) >) tiv (P) Z = KioRio (p) Dio (p) Zs + 
v=1 


vei 


+ mDio (p) Dip (p) KiRic (p) hi + MK icRioD ip (p) Dio (p) fin (24) 
It is evident from (24) that if the system remains stable at m - 0, that is,at Kjp ~ ©, then 


Bm Sag = '% gt: (25) 

Thus, in this case also a disturbance in the ith circuit loop can be suppressed by an increase of the amplifica- 
tion factor of the ith regulator. It follows from (25) that together with the disturbance also the load interaction fj}; 
is removed, and thus the multiarticulate (interlocked) system is resolved (with an accuracy to ¢) into a system of n 
independent regulation systems, The disturbances in the other circuit loops exert no direct effect upon the circuit 
loop under consideration, only an indirect effect by way of the other regulated quantities x». Obtaining of autonomy 
in this case simultaneously removes the influence of “foreign” disturbances. 


It remains to examine the question of stability of the resulting system of multiarticulate regulation under the 
condition that all Kj, ati=1,...,mtend toward infinity. The determinant of the system will be as follows: 


Ay, (p), By (p) G12 (p), By (P) ais, . -, Br (p) ain (Pp) 
By (p) 3 (Pp) Ave (p), Bs (p) aes, - - -» Bs (Pp) Gan (P) (26) 


By, (P) Om (P), Bn (P) Ona (P), - - -» Ann (P)s :; 


A= 


Au (P) = [mDyo (p) Din (p) + KioRio ()] Dio (P), 
B, (p) = mDio (p) KioRio (p). 
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The characteristic equation which determines whether the entire system is stable or not is as follows: 
A= 0, (27) 


Equations of the type of (27) were encountered already in [10]. It was shown there that, if the systems for 
every regulated quantity remain stable at m - 0, then also the complete system in its entirety will be stable, provid- 
ed that the degeneracy equation of the complete multiarticulate system satisfies the stability conditions and that 
auxiliary equations of the first, second, or third kind too satisfy the stability conditions [10]. 


It follows from the results of the preceding section, that it is possible to stabilize every separate system by 
introduction of additional amplifiers with sufficiently high amplification factors provided with negative feedbacks 
having transfer functions of the type 


pl or ae a 
1+T,p ap*+bp+1 


It is not difficult to see that in this case the stability of the degenerate system will fully depend upon the 
parameters of the components which are included in the feedback lines of the amplifiers. Hence, stability of the 
degenerate system can always be secured. 


For the case when the auxiliary equation is of the first kind, it satisfies the stability conditions automatically; 
if the auxiliary equation is of the second or third kind, it is always possible to secure satisfaction of the stability con- 


ditions by the auxiliary equation by a suitable choice of the parameters in the regulators and in the amplifier feed- 
backs. 


Thus, in this case it is possible to secure stability, and hence feasibility of real applications of the described 
system. 


SUMMARY 
On the basis of the investigation carried out by us, the following conclusions can be made: 


1. If in a closed dynamic system not all components are subject to the influence of disturbances from the 
outside, then the suppression of the disturbances can be achieved by an increase of the amplification factors of those 
components which are not directly subject to the action of the disturbances. And this is brought about by the fact 
that the disturbances from any given component are suppressed by the amplifications in the system components pre- 
ceding that particular component, starting the count from the input to the system of the useful, or control, signal as 
the beginning. 


2. It is feasible to turn such systems into realities, if their structure is of the class which remains stable when 
an unlimited increase of the amplification factors takes place. 


3. The obtained results cover not only systems with a single regulated quantity, but also systems of multi- 
articulate regulation, And at that, as should be expected [8], the multiarticulate system, simultaneously with being 
freed from the effects of disturbances, is also dissociated into separate autonomous systems (with an accuracy of 
autonomy up to €). 


LITERATURE CITED 
V. S. Pugachev, Application of Canonical Resolution of Random Functions to Determination of a Linear System 
Optimum. Avtomatika i Telemekhanika, 17, No, 6 (1956), 
A, A, Fel'dbaum, Electrical System of Automatic Regulation. Oborongiz, State Defence Press (1954). 
N. Wiener, Extrapolation, Interpolation, and Smoothing of Stationary Time Series, John Wiley (1949). 
Ya. Z, Tsypkin, Pulse Systems, State Phys. and Math, Press (1958), 
A. A, Fel'dbaum, Computing Devices of Automatic Systems. State Phys. and Math, Press (1960), 
A. A, Rizkin, Fundamentals of Amplification Circuit Theory. Soviet Radio Press (1954), 
B, Ya, Kogan, Experimental Electronic Models and Their Application for Testing of Automatic Regulation 
Systems. State Phys. and Math. Press (1959), 
8.  M. V. Meerov, Principles for Design of Automatic Regulation Systems with Small Stationary Errors. Avtomatika 
i Telemekhanika, 10, No. 2 (1949), 


_ 
>. 


IAAP OW 


1189 








10, 


11, 


13, 


14. 


1190 


M. V. Meerov, On Automatic Regulation Systems Stable with Amplifications as High as Desired. Avtomatika 

i Telemekhaniks, 7, No. 4 (1947), 

M, V. Meerov, Synthesis of Structure of Automatic Regulation with High Accuracy. State Phys. and Math. Press 
(1959). 

M, V. Meerov, Systems of Multiarticulate Regulation Stable at Small Stationary Errors. Proc, of 2nd All-Union 
Conference on Automatic Contro] (1955) vol. 2. 

M. V, Meerov, On Autonomy of Multiarticulate Systems Which Remain Stable with Unlimited Increase of 
Established Accuracy, Avtomatika i Telemekhanika, 17, No. 5 (1956). 

M, V. Meerov, On Stability Against Disturbances of a Certain Class of Dynamic Systems. Doklady Akad. Nauk 
SSSR, 139, No, 4 (1961). 

M, A. Aizerman, Certain Structural Conditions for Stability of Automatic Regulation Systems. Avtomatika i 
Telemekhanika, 10, No, 2 (1948), 





All abbreviations of periodicals in the above bibliography are letter-by-letter transliter- 
ations of the abbreviations as given in the original Russian journal. Some or all of this peri- 
odical literature may well be available in English translation. A complete list of the cover-to- 
cover English translations appears at the back of this issue. 











THEORY OF CONTROL SYSTEMS WITH LIMITED-SPEED 
SERVOMECHANISMS* 


L. S. Gol'dfarb 


(Moscow) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No, 10, 
pp. 1324-1332, October, 1961 

Original article submitted March 15, 1961 


A paper [1] concerning investigations of the stability and self-oscillations of automatic control 
systems having a nonlinear element — a limited-speed servomotor with rigid negative feed- 
back — was submitted to the IFAC First International Congress on Automatic Control. The 
analysis was performed by using the phase plane method as well as the describing-function ap- 
proach. In this, the nonlinear element was considered as an elementary factor, whichrendered - 
the describing -function approach more complicated and reduced its effectiveness and also im - 
paired the accuracy and clarity of the investigation results. 

In the present article, the detailed structure of a nonlinear element is considered, which 
makes it possible to simplify and elucidate the investigations, to obtain some new results, and 
to consider more complex cases of practical importance. 


Paper [1] is concerned with an automatic control system (Fig. 1) which consists of two inertial elements and a 
limited -speed servomotor, which is encompassed by rigid negative feedback. The input quantity E of the nonlinear 
element is related to its output F by the relationships 


E=F fa |S |< Fin 
aF 
dt 


dF 
at 


F,, for FE, (1) 


= — Fs for | ~F 


where Fin is positive and corresponds to the maximum rate of change in F that the servomotor can secure. 


In [1], the system's motion is investigated by using the phase plane method and the describing -function ap- 
proach. The article provides a comparison between the values of self-oscillation periods for equal time constants 
of the inertial elements in dependence on their gains SM, These values were calculated by using both methods, It 
is stated that the describing -function approach secures a satisfactory accuracy. 


The method used by the authors of paper [1] for the mathematical description of a limited-speed servomotor 
with rigid feedback leads to the fact that the nonlinear element's equivalent gain vector (1) contains an imaginary 
component. This is connected with the fact that the element under consideration includes a nonlinear element with 
a “limited power” characteristic (NE IV in [2]) as well as a linear element — an integrating factor. 


The complexity of the equivalent gain complicates calculations and makes them less clear. 


This drawback can be avoided by separating the nonlinear element in the manner shown in Fig. 2. By con- 
sidering the nonlinear element's structure in detail, it is possible greatly to simplify the calculations and to make 
them clear in the highest degree for the case considered in [1] as well as for systems with more complex servo~ 
mechanism dynamics, whereby the investigation method can be extended to a larger number of practical problems, 


* The article was prepared for publication by E. B, Pasternak. 
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The block diagram in Fig. 2 apparently provides a more accurate description of the properties of a system with a 
limited -speed servomotor that is encompassed by negative feedback. 


The servomotor is characterized by the maximum rate F,, and by the width a of the linear zone. The gain in 
the linear zone is 


F .: 


The mathematical description of the rigid feedback loop Wz, which corresponds to NF in [1], coincides with 
(1) if # tends to infinity. 


Actually, for ~fm< F < F,,, the transfer function of the loop 


| J 
gah Fg ee 1)" 


vlfF 


tends to unity with an increase in #. In this, F = E in the limit. It is obvious from Fig,2 that the other conditions 
are satisfied, 


The consideration of the degenerate case where ! = » enabled the authors of paper [1] to reduce the analysis 
of a third-order system to the analysis of a second-order system, which made it possible to consider its motion in the 
phase plane, However, by determining the equivalent gain according 









































NF to (1), they also limited the describing -function approach to this case 
® — %[ 5 E ~- F of maximum idealization, which apparently was not necessary. 
@ T,p*! = In the present paper, such limitation is not imposed, and the 
limiting case is considered separately at the end of the article. 
Pi The transfer function of the linear portion of the system is 
inp : 
W; (p) => [1 + Wi (p)], @) 


Fig. 1. Block diagram of the nonlinear 

systemi of [1], where W;(p)= SM/[(Tmp + 1) (Typ + 1)] is the transfer function of 
the two series-coupled inertial elements, and Tm and Ty are positive 
real numbers. The corresponding complex gain is 


SM 


Wi (0) = GF jeT.) dF 10Ts) * 
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Let us consider the following parameters: 
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Fig. 2. Block diagram of the control 
system, 

Re Wi = (1 — 82") SM Im Wi= 2dS MQ) (5) 
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SMVT ale w 1: 
W, (j2) = SMV Pao wy, (tay, 


1 
Wo (2) = — a7 + Ga (7) 





is the normalized complex gain, the hodograph of which (Fig. 3) represents the hodograph of the normalized com- 


of the series-coupling of two inertial elements, which is shifted by a 





1 1 
plex gain —=-W1 (8) = =o) F jean 
to the right of the imaginary axis. 














Fig. 3. Hodographs of Wo(jQ) and Wgo(ja). 
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d=1; SH=10; T,, =Q) 


12) 


: 1 
Fig. 4. Construction of thew, (j)and ~ Freie eoeraree 


The vectors W ,(jQ) are rotated through —1/2 with respect to the vectors Wg(jQ), and 


Re W, (ja) =SM¥ Tals tm W, (j0), (8) 


Im W; (ja) — = “4 ¥?oTs Ro w, (2). 








The vector of the nonlinear element's equivalent gain Jyféyy{[2] contains only a real component. Thus, the 
possible self-oscillation regimes of the system (Fig. 4) 


‘ i 
W, (jQ) = — (Q + 00) 


A 
Iyérv(F) (10) 
correspond to the roots of the equation 
Im W, (2) = 0, (11) 
whence 
Re W, (72) = 0. (12) 


Thus, the determination of the frequency of the possible self-oscillations can be reduced to the determination 
of frequencies Q, and Qz at the intersections of the Wy, hodograph with the imaginary axis. The frequencies Q, and 
M, are determined in an elementary manner. For this, the complicated calculations and graphical constructions 
presented in [1] are not necessary. The appearance of the Wy hodograph in the third quadrant of the complex plane 
is a necessary condition for the presence of self-oscillations only for p # «. 
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It should be noted that the frequencies 2, and Q2 of the possible self-oscillations are determined only by the 
parameters SM and d of the inertial elements and that they are independent of the servomotor parameters. Let us 
write (12) in explicit form: 

{ i—Q 
sm + (1 — Q*)? + 4a8Qs 





=0, 


(4 — Q*)* + 4d°Q? + SM (1 — Q*) = 0, (13) 


Let us introduce the parameter p =<... which was considered in [1]; it is obvious from Fig. 3 that o >Q> 1. 


Consequently, O<p<rf. 


Dependence (13) is represented graphically in Fig. 5. The minimum value of SM that secures the appearance 
of self-oscillations is given by 
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SMe, = 4d (d + 4). (14) 
In this, 


Qi = QF = 2d + 1. 


If SM > SM, there are two different intersection frequencies 1 <Q;< 2d +1 <Q: 








Q,=Y “+1 —2— lVISM —M dF il (sm — Md d—M), (15) 





(16) 





Q,=Y +1-2+ VSM — d+ iisM —M@d—ipl. 
The hodograph of the complex gain W; of the linear portion of the system intersects the real axis at Q = Q, 
and 2 = Q2. Let us determine the length of the thus obtained segments (Fig. 4). 
According to (5), (7), and (8), 


2dSM YT,T, 
(i — By + hor 





Re W; (jQ) = — 


If we take into account (13), the denominator of this expression at the intersection points is equal to SM(Qj,2— 1), 
whence 


Re Wy; (j2,) = — a te 
ei? pop ae (17) 
while 
Re W; (32) < — VT,,T,<Re W7(jQs). (18) 
If SM increases, 


Re W, (j2;) +—oo, Re W; (jQ,)~0. 


In the linear approximation, according to Nyquist's criterion, the system is stable for 


—1<pRe W, ((Qy), 
and for 


— 1 >p Re W, (jQ,). (19) 


The meaning of (19) is obvious, The transfer function cP + 1)" of the servomotor loop corresponds to an 


inertial factoz with the time constant 1/ and a unit gain. In the linear approximation, the system is equivalent to 
a system consisting of three inertial elements, which is stable for relatively small as well as relatively large time 
constant values of one of these elements. 


With an increase in , the system becomes equivalent to a system consisting of two inertial elements, which 
is stable for any parameter values, 


In the presence of limitations, the system loses this property and becomes in this case a self-oscillation system 
with a severe excitation regime. The stability of the system's equilibrium position “in the small” does not ensure its 
operating ability, which once again emphasizes the necessity of taking into account the nonlinearities of the charac- 
teristics in determining the stability of even such simple systems. 
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Let us write (10) in explicit form by using (17): 
= VF nls ; (Q =Q, Q= Qs), 


er cc ae 
F 
Jnéav = —~J Eve (0<J nov, S 1). 


According to [2], 





4fi i 
JNEive = <|[ya- + |, 


7 a = arc sin; A; is ae of self-oscillations at the input of the nealing element. The hodograph 
Ay 
~fzrc corresponds to ER Ee section of the real axis (Fig. 4); the Siig curve is given in Fig. 6. 
B — InEiv, — 


The equation of the system's self-oscillations is 


ar eVit. 1 


ome te : 0. | 
1,2 a(@,—1) meg a= (20) 





Three regular cases are possible (Fig. 7). 

I, Kg < Ky < 1. As in the case where SM < SM¢,, the system is stable in the linear approximation and also 
if the speed limitation is taken into account. The hodograph of the complex gain of the system's linear portion inter- 
sects the real axis to the right of the point 


Il, Ke< 1 < Ky. In this case, the hodograph of the complex gain of the system's linear portion intersects 
only once the hodograph of the nonlinear element's equivalent gain. 


| 











Fig. 8 Fig. 9 


The intersection point K,; = Sng” CTesPonds to a table limiting cycle [2]; the equilibrium position in the 


linear approximation is unstable, A self-oscillation regime with soft excitation is observed in the system. Figure 8 
shows the phase pattern of the system's model for this regime, which was obtained by means of an MN-7 electronic 
simulator. 

II, 1 < Ky < Ky. In this case, there are two intersection points, of which the point Ky=—~ corresponds to 


a stable Limiting cycle, while the point Ky = 5—— corresponds to an unstable limiting cycle. On the basis of (19), 
i) 
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is to 








the equilibrium position is stable, The system operates under the regime of severe excitation of self-oscillations. 
The corresponding phase pattern for the system's model is shown in Fig. 9. By using Fig. 7, the self-oscillation am- 
plitude can readily be determined from Eq, (20). 


In conclusion, we shall consider the particular case where a = 0, which was considered in [1], For this, the 
system becomes a relay system, The transfer function of the linear portion 


SM 
W, (P) = + sa TFT 





contains the term 1/p, where the order of the numerator is only by unity lower than the order of the denominator. 
In calculating such systems, the describing-function approach in its usual form can result in considerable errors, which 
can be avoided by taking into account the “jump” phenomenon [3]. 


A more accurate complex gain for the servomotor's feedback circuit can readily be determined analytically 
by taking into account the jump and the exact signal shape: 





W 5. =a (—Zi+ y ep) = a4 (21) 
The more accurate complex gain of the system's linear portion is 
Wg. = bial AF Wo.» (22) 
where 
Veo = t (t—sw tam - (Faw tot Wo. (23) 


Thus, the oer of the jump leads to a somewhat greater — — — ofa) shift of the nor- 
malized hodograph => um W100) to the right of the imaginary axis and to a bh 
the real axis (Fig. 3, a line). 


im sm Mt f pee 
Since, in this case, 

aan te eu 4 J. 
Jné= lim JNavy = =". 


the frequency and the amplitude of the limiting cycles can readily be determined analytically from equations 
similar to (13) and (20): 


(1 — 29° + 4a" — © sm (Q*— 1) = 0, (24) 
Fin WT ml ndQ 
Ay =—* SG" * (1+ goa): (25) 


In [1], the describing -function approach was used for determining only the half-period p =< for d= 1 and for 
six values of SM. 


The diagrams of the half-period and amplitude of self-oscillations in dependence on the system's parameters 
were plotted (Figs. 10 and 11) on the basis of (24) and (25), 


For the sake of comparison, Fig. 10 also shows the results from [1], which were obtained by using the phase 
plane method and the describing-function approach. 


It should be noted that the exact results for the limiting case can be obtained in a much simpler way by using 
the frequency method of investigation relay systems [4] instead of using the phase plane method. 
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results obtained in [1] by using the describing-function ap- 
proach, 
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It is demonstrated that an investigation of unitary-code automatic systems of the type under 


study can be reduced to an investigation of analog systems that contain equivalent nonlinear 
elements. 


A method is given for analyzing unitary -code automatic systems. 


An analysis is made of the accuracy of a specific unitary -code system, and recommenda - 
tions are given for increasing the accuracy of such a system. 


With the development of mathematical machines digital computers are finding an ever wider application in 
automatic control and regulation systems. 


Digital automatic systems in which digital data correspond not to the total values of the physical quantities 
but to their increments shall be called unitary-code automatic systems (UCAS). Such systems are finding the widest 
application in industry. 


Notwithstanding the wide usage of UCAS, the problems of the dynamics in such systems have not as yet been 
adequately studied. The solution of the problem is complicated by the fact that the methods for analyzing conven- 
tional pulse systems are inapplicable to UCAS, 


In this paper we shall demonstrate that the investigation of UCAS of a given type can be reduced to an in- 
vestigation of continuous systems containing equivalent nonlinear elements. We shall estimate the accuracy of a 
specific unitary -code system and point out a method for increasing its accuracy. 


The Equivalent Block Diagram 


The block diagram of the UCAS is shown in Fig. 1, where 1, 5 are coding units, 2 is a reversive counter, 3 is 
a decoder, and 4 is the analog section of the system. 





In this setup the coding units 1 and 5 convert analog data into unitary codes (Fig. 2), and therefore the signals 
Uy and U, represent trains of pulses with an identical height and different time intervals between adjacent pulses. 


The reversive counter with the carry circuit is a comparison block. The comparison operation in the counter 
consists of the following: the increments of the quantities are applied in the form of pulses to the reversive counter 
where the increment pulses Uy are added and the increment pulses Uy are subtracted. As a result of this the number 
indicated by the counter is equal to the difference rUy — LU, at each instant of time. 


The decoder converts the binary code (the outputs of the reversive counter triggers) into a voltage whose mag- 
nitude is proportional to the number of pulses EUy —DUx. 


The block diagram shown in Fig. 1 can be converted to the form shown in Fig, 3, where ¥[y] and ¥{x] are 
equivalent nonlinear elements of a special type whose characteristics are shown in Fig. 4a; W(s) is the transfer func- 
tion of the linear analog section of the system. 


It is evident that the equivalent nonlinear elements ¥{y] and ¥[x] can be treated as a parallel connection of 
linear and nonlinear sections (Fig. 4b) [1]. Therefore the block diagram shown in Fig. 1 can be reduced to a dif- 
ferent form (Fig. 5) when the conversion coefficient K, is taken into account. 


Thus the investigation of unitary-code automatic systems can be reduced to an investigation of continuous 
systems containing equivalent nonlinear elements of a special type. 
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Estimating the Accuracy 
We shall use k(t) to denote the pulse transient function for the limiting linear analog system. Then in accordance 
with Fig. 5 we can write 





oo 


z(t) = ty (t — 0) + Yale, y, t — Ok (x) de. (1) 


For the system error we obtain the expression 


co 


e (t) ={y (t) —\ y(t — t) k(x) ax} + | Mn (z, y, t — t) k (x) de. (2) 


ou"8 


The first term in the right side of Eq. (2), which determines the dynamic error ¢4(t) of the system, can be 
computed in the conventional manner. The second term, which determines the error caused by level quantization, 
is defined in accordance with an upper tolerance. 
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-From Fig. 5 it is evident that ©,,,(x, y, t)= x, t)—@(y, t), From this we obtain the expression for estimating 
the absolute value of the function ,,(x, y, t): 


1Madz,¥, |< 1/9 (z9/ +19, 4). 





ating 


But the difference between the absolute values of quantized and unquantized quantities cannot exceed one-half 
the quantization step 4/2; i.e., | (x, t)| = 4/2; | ey, t)| = 4/2 (cf. Figs. 4a and 4b). As a result we obtain the 
relationship 


Paz, y, | [| D(z, )|+|(9Y, | <A. 


Therefore the maximum error caused by level quantization is estimated as 


t t 
| ask!) maxl <\| Pant, y,t—D| lk(t*)|dt< al k (t) | dr. 


In a steady-state mode we shall have the inequality 


co 


| aro) max! <4{ [k (t) |r. (3) 


From this it follows that if the pulse transient function k(t) is non-negative, then the maximum error caused 
by level quantization will not exceed a value equal to the value of the product of the quantization step and the value 
of the transient function at the same instant; the maximum steady-state error produced by level quantization will not 
exceed the value of A. 


In the case where the pulse transient function changes sign we have 
l€af°)max| SNA - = (n> 1). (4) 


We shall use €, to denote the error which arises during the process of preparing the program (i.e., the ap~- 
proximation error during the preparation of the information). Then we shall obtain the following expressions for the 
maximum total error of the UCAS: 


| &z (t) mand < | &a(t)| + | fan(t) max! + | eal, (5) 
Jes (0°) max! S| ea (oo) | + | fart) maxi + | ea). 


Analysis of a Certain Unitary-Code System 








We shall study a unitary-code system for program control of metal-cutting machine tools; its block diagram 
is shown in Fig. 6 [3]. The transfer function for the analog section of the open-loop system is given by the formula 


NW (s) = KyW, (3) We (8), (6) 


i+7; 
W, (s) a 


N 


Kra Kud*tok rc* 
Wy (s) = —EPasd + Toe) C4 Tos) 
ae 1+ KRa KT GK recs 
(1 + Tys) (1 + Tas) (1 + To) 














Here N = 1/47 is the transfer number of the reducer, T, = 0.024 sec ; Tz = 0.2 sec ; Ts = 2,0 sec ; Kpg = 8.8; 
KMo = 2.0; Krg = 3.4; Kpc = 42. 


The following quantities are given: the conversion coefficient for the unit which converts analog quantities 
into unitary -code quantities Kc; = 100 units/mm; the conversion coefficient for the device which converts code 
quantities into a voltage (the decoder) K,2 = 2 v/units; the value of a pulse 4 = 0,01 mm; the approximation error 
&g = 0,005 mm; and the input signal with respect to the Y coordinate y(t)= lt. 
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In accordance with the graph for the transient response of the reduced limiting linear analog system (Fig. 7a), 
1°) 
[|k(e)|de = 1.68, 


0 


| Gar 0) max linv< A\ | k(x) Jde = 0.04168 uw. 
0 


As we know [2], the error € 4 (@) is determined from the expression 





& (0) inv= Coy (t)+ ¢, 20 : (7) 
In this case 
Cp, =0, C, = - an aie == 0.013 200 , 
NK, Kara ® Mo 7 
dy (t) 


y (t) = it, <a = 1: 


Substituting Cp, Cy, y(t) and dy(t)/dt into (7), we obtain €g (@)iny = 0.013 mm. 
Therefore 


| €z (00) maxlinvS | & d()inv| + | fark°o)max linvt | &a| = 0.0348 nm. 


We see that the total error |é>(~)maxliny is large. It consists chiefly of the errors €an()max and eq (~). 
Therefore a reduction of these components is of great significance. 


From Eq, (3) and (7) it is evident that a reduction of the oscillatory nature of the transient response leads direct- 
ly to a reduction of the error &47(t)max: am increase in the gain leads to a reduction of the error € q(t). 


Since the gain Kp = KraKMo K7GKgc in a given system is sufficiently large, it can be assumed over a wide 
frequency range that 
Lm | W, (jo): ~ 0, arg Wy (jo) = 0. (8) 
This is evident from the graph of the logarithmic frequency response for the equivalent open-loop analog 
system (Fig. 8). 
Therefore in the intermediate frequency range which characterizes the transient response of the system the 
transfer function for the analog section can be written as 


NW (5) = Kad; (8) Ws (8) = Kas (8) = Kan te, (9) 





and the transfer function for the equivalent analog section of the system can be written as 
1+ T. 1+T 
Koy Keg NW (8) = KerKe Kana =B a, (10) 
The characteristic equation. for the equivalent closed-loop system will be 


# + BT,5s + B = 0. 
The roots of the latter equation are 


a7? 
dg th Oe gy —B=—p+VB — B. (11) 
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Fig. 8 
In our case 


B= Key KeoKan= 0.297 < 4, 


and therefore the roots 4, Ag of the characteristic equation are complex, 


Thus it follows from expression (11) that the coefficient 6 increases (this causes a more rapid attenuation of the 
oscillations) as the angular frequencies of the oscillations decrease. Therefore as gain 6 increases, the error 
€an(t)max decreases. 
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In that case the roots of the characteristic equation will be real when 8 =1, 


It should be noted that Eq, (9) and (10) will not be valid for 8 > 5.02. This fact can be established from the 
gtaph of the logarithmic frequency response (Fig. 8). In fact, for 6 > 5.02 the phase and amplitude margins become 
inadequate. Therefore the system will be either less stable or unstable. 


In accordance with the graphs for the transient responses (Fig. 7b): 
for 8 = 1 





| fal}max | <A { | k(x)|dv = A (1.435 + 0.135) =0.0127 nx, 


eg (001) = Cu 2 — 0.0034 mm; 


for 6 = 5,02 


|fanlolmax fe < A | | (x) | de = A (1.04 + 0,04) = 0.0108 sux, 


&q (co): = Cu EO ) _ 0.00068 us. ; 


The overall regulation for 6 = 5.02 is insignificant; here the errors eg(t)z and €,,,(t)4x are much less than 
Ed(t)inv and Ep iny(t)max- Therefore the total error | € 5, (t)max | is less than | 65 jny(t)max | - 


In a steady-state mode: 


| €x1 (9°)max | < | €ank-) max + | &q, (00) | + | &a | = 
= 0.0211 sm < | €x5,,(0°) maxl, 

| €x, (0°) max] <| &ane(°)max | + | fae (00) | + | &2 | = 

=0.01648 wm < | ex: (0°) max | < | exp f°)maxl: 


The transient response time T, for 6 = 5.02 is appreciably smaller than Tiny [Tiny = 13 sec , T; = 5 sec , Tz = 
= 1,2 sec (Fig. 7)). 


Based on the above, we conclude that in the system under study it is expedient to choose the gain within the 
limits 1 < 6 = 5,02, 


Thus in our unitary-code automatic system, which has internal acceleration feedback, increasing the gain of 
the equivalent open-loop system is an effective method for improving the performance of the system. 





SUMMARY 
The investigation of unitary -code automatic systems of the type under study can be sedueed to an investiga- 
tion of analog systems containing equivalent nonlinear elements of a special type. 


In the case where level quantization can be neglected, UCAS can be treated as linear analog systems, The 
method cited above for estimating the accuracy, and the method for increasing the accuracy can be treated as a 
first step in the investigation of unitary -code automatic systems. 
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In the majority of works devoted to the determination of self-oscillation parameters it is as- 
sumed that the parameters of the members in the linear part of the system stay constant with 
respect to time. However, many relay systems contain a member in which the parameters 
take on different values for the closed and open positions of a relay (for example, systems 
with relay servomechanisms, vibrating regulators of electrical machines), Methods of deter- 
mining the frequency characteristics for several particular cases of such systems are given in 
{1, 6}, 

In the present article a method is presented for determining self-oscillation parameters 
based on the principle of harmonic balance [2] in a relay system containing a member of the 
first order with variable parameters, which is connected immediately after the relay element. 


1. Statement of the Problem 
In Fig. 1 the block diagram for the relay system of an automatic control is represented. In this illustration RE 
is a relay element having the characteristic depicted in Fig. 2; 1 is a section with variable parameters, and the 
transfer function of this section has the form 


4 z y 4 
“ee FL PL REF! 12 | 
W, (p) = Tet for jz|=86 


Fig. 1 W,(p) = aT for z=0; 








2 is the “invariable” part of the system with a transfer function W(p). 


Assuming that self-oscillations take place in the closed system, we 
| will determine their parameters—the period, the amplitude, and the dura- 
r tion of the relay closure. 

a 














-a 
t A It must be noted that despite the linearity of section 1, in the closed 
8 automatic control system it cannot belong to the linear part of the system 
= os because its parameters depend on the parameters of the self-oscillations 
(the length of the period and the amplitude). 
43 We will limit ourselves to a consideration of symmetrical self-oscilla - 


tions, and for this case we will find the equivalent complex amplification 
coefficient of RE and of section 1. 


2. Determination of the Output Signal of Section 1 
~ The output signal of the relay element has the shape depicted in Fig. 3. 
We will designate the period of the self-oscillations by Tp, the duration of the relay closure by yT», The 
Output coordinate of section 1 is the steady reaction of it to a series of rectangular pulses with alternating signs. 








r To determine this reaction the results of pulse system theory can be used. 
A similar approach to the analysis of relay systems is used in [5]. In the 
present case section 1 has variable parameters, and so we will take advan- 
tage of the results in reference [3] for the determination of its output co- 
ordinates. First we will convert to dimensionless time f = t/T». We will 






































t introduce also the designation ¢ = t—nT», where n= 0, 1, 2,3,...,0 = 
se=1, 
— 1 ~~ On the basis of [3] the discrete transfer function of the pulse system 
v0 formed by section 1 is written in the form 
yaad Wi (qe) = &( 1-2 = ew) 
,o) = —_ ' ‘ 
1(9 ( a BE (<e<71) (2.1) 


(i —eFr) eP'e—v) 69 





Wi, (gq. &) =k rae (1<e<i), 
T pot ' 
B=7, B=, 4= Br +B (1 — 7). 


The coordinate x can be represented as a multistage function 
a(n) = Bcos an. (2,2) 
Based on (2,1) and (2.2) we get for the output coordinate of section 1 for steady operation the expressions 


(2.3) 


y! (n, e) = Bk ( EMRE Se Bt 


—*) — 1)" 0<e<y), 
ite FE , 


ea e~ Fr) e Few) 


H(n, e) = Bk (— 1)»—! 1). 
y" (n, @) (— 1) p07 E (r<e<1) 





3. Complex Amplification Coefficient for the Non-Linear Part of the System 
The equivalent complex amplification coefficient of the RE in common with section 1, according to [2], has 
the form 





I (@, A) = g (A, @) + jb (A, @), (3.1) 
where w = «/To, 
B’(i—v) 
g (A, @) =3(\y¥ \y i (n, e)e * sinnede + 
0 
I+y 
2 B (i—y) 1 B’(i—v) 
oly \ y(n,e)e * sin aede + \ y(n, e)e * sin neds | : On 
i= +r 
2 
1-1 
2 B’(l—y) 
b (A, @) = =| \ y(n, 2) e * cos nede+ 
0 
itv 
2 B (1—y) 1 B’(1—+) 
+ \ y'(n,2) e * cosnede+ \ y!!(n,e)e * cos xe de |. (3.3) 
Ir ity 
2 2 





1) 


2) 


3) 


2) 





Carrying out the integration in (3.2) and (3.3) we obtain 




















_ 2B k (x* — BB’) (B’ —B) (1 —e*) (1 4 e-F'A— 
a RA TPT | 14+e 8 a: a cos +! + 
By) o—P'I— 
+ (x*B* + xB* + 28°84) sin St — a os! Sashes: +]. tun 
2B A P PF —v) 
a FBV RS | — 88" — ts rene x sin “} — 
— (x? + BB’ ) (p ob 8’) sin 2 aa f (1 —e~?r) (1 +e ene n® cos 2% ° 
ite 2 (3.5) 


Taking into account that when the signal 4= A sin —t is acting on the input of the RE we have the relation 
0 


2 
7*2 arc Cos 























“A? 
1 ’ pa 
and also that 8 = —, B’ = SF We obtain 
n(1—Y) 
2B k @ (w*TT’ —1)(T —T") (1i—e “FY (1467 re 
§ = 5A (i +-o°F)(1 4-09") * as" 
i+e _ “TE 
+ (T?@* + To + 2))/ 1-4 - 
= SS) 
— w* (7? — ry 1a te e a . |; (3.6) 
i+e oT; 
_2Y = __ afi—v) 
b = 28 k @ (@*T'T’—1) (T—T") (. Ti, oF 
™ 9A (i +o" \(i+0°7"*) ——* x 
i+e oTE 
xV 1-5 —ett+ TT’) (T +7 VV 1— + 
(oF) FY 
+ i—e i+e = (T* — T"*) @ 5 (3.7) 
1i+e TE 
TT’ 
where Te = FET U=) * 
We will consider a series of special cases of these expressions: 
a) o — 0, g~ Vis oe bx~0. 


Moreover, it is easily seen that I(w, A) is the equivalent complex amplification coefficient of the RE. 
b) @ + oo, g +0, 6 0, consequently, / (w, A) — 0. 
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c) <1, 7 +0,/ (@, A) +0. 


. = 
d) _ 4Bk V = Aa 














A ~%st=1. 6 =—aarerm ~ 
b- —2Bky/ 4S _ oT _ 9 
e) sa nA At a oo Fh. 
’ ” = 4Bk am 4SBk oT a* 
T=! .8=ZAeTTT srt V 1- Re o=—SA writ i—A - 


This case corresponds to a series connection of the RE and the section with constant parameters (k'= k, T* = T), 


a*® a ) 
4Bk V <6 ye oT th th ( are cos fo 


T’ — ov, a= -—— 








f) nA 1+ oT? : 
El Ae S+4 + tb (recon * fur) 
Ts 1+ #7 . 


This result concurs with the result in reference [1]. 


Expressions (3.6) and (3.7) obtained for the equivalent complex amplification coefficient permit the problem 
of determining the self-oscillation parameters to be solved by the method of harmonic balance. 





























Fig. 5 





The starting expression for this calculation is the 
well-known relation [2] 





1(@, A) = Woo - (3.8) 


A solution of equation (3.8) is produced graphically by the following means, 


In the complex plane (g, jb) a family of root loci I(w, A) is plotted, each of which corresponds to a fixed value 
of w (Fig. 4)? The construction is carried out by means of calculations of the components of Kw, A) according to for- 
mulas (3.6) and (3,7) when A is varied from a to + for every value of w = const, A root locus of the function 
—1/W(jw) is constructed on these same coordinates, The self-oscillation parameters are the values of A and w corre- 
sponding to the point of intersection of the root locus —1/W(jw) with that one of the root loci Iw, A) for which the 
value of w corresponds to the value of w of the function —1/W(jw) at this point. 


The stability of the self-oscillations with the parameters determined in this way is verified by the same meth- 
od as in the case of a non-linear system with constant parameters [2], 


If the characteristic of the relay element has the shape shown in Fig. 5, the corresponding root locus I'(w, A) 
can be obtained from the root locus I(w, A) by rotating all vectors through the angle 


* Note that from cases c) and d) it follows that the root loci 1I(w, A) are closed ior all values of w. 
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= T). 


(3.8) 


| value 
to for- 


corre~ 
h the 


meth- 


), A) 





Ay ¥ (oxy — cg) / 2w, 


where a, = arc sin(a/A), %, = arc sin (b/A)and by substituting ag = A sin (a+ «2) /2 (see [1). 


Thus the problem of finding self-oscillation parameters in a relay system containing a member with stepwise 


changing parameters comes to the construction of a family of root loci 1(w, A) and a root locus of the inverse am- 
plitude -phase characteristic of the constant part of the system. 
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cover English translations appears at the back of this issue. 
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The paper studies the functioning of finite automata in media that respond to their actions in 
a random fashion, The concept of the expediency of automaton behavior is introduced and 
examples are cited for automata which have expedient behavior. 

The behavior of automata in media whose probability properties depend on time in ac- 
cordance with a Markov chain is studied. ° 


Introduction 
In this paper we study the functioning (behavior) of finite automata in media that respond to their actions in 
a random fashion. 


It is assumed that an automaton perceives any response of the medium as a reproof ("penalty") for its action or 
as the absence of such a reproof (a “non-penalty"); it is also assumed that the medium “penalizes” the automaton 
with a probability that is specified for each of its actions. For these assumptions it will be demonstrated that the 
behavior of a finite automaton in a stationary random medium is described by a finite Markov chain; the magnitude 
of the mathematical expectation of the penalty serves as the measure of expediency for the behavior of the automaton, 


The paper cites examples of automata that have expedient behavior. In particular, an automaton construction 
(an automaton with a linear tactic) is indicated for which the mathematical expectation of a penalty decreases as 
the memory capacity increases and coincides in the limit with the minimal possible expectation in the specified 
random medium. Such an automaton performs (with a probability close to unity) that action for which the probability 
of a penalty is minimal in the specified medium. 


Then we study the behavior of finite automata in a random medium whose probability properties depend on 
time in accordance with a Markov chain. We investigate the behavior of an automaton with a linear tactic in this 
random medium. It is demonstrated in particular that the mathematical expectation of a penalty in this case reaches 
the minimum for a certain fixed memory capacity and increases as it is increased further. The existence of an op- 
timal memory capacity in a given medium for automata with a linear tactic makes it possible to design automata 
which have the most expedient behavior in the specified random medium. 


The paper cites graphs that indicate the dependence of the mathematical expectation of a penalty on the 
memory capacity and parameters of the medium; we also provide a table which specifies the optimal values of 
memory capacity and the corresponding values of the mathematical expectation of a penalty for various random 
media. 


1. The Behavior of Automata in Stationary Random Media 
Assume that the finite automaton A is described by its canonical equations*® * 





gp (¢ + 1) = D@ (d, s(t + 1)), (1) 


* A brief presentation of the results is given in [1]. 

** The definition of a finite automaton used here coincides with the definition of a non-primitive network cited in 
(2, 3}. It is not difficult to arrive at the same definition if we start from the definition of a restrictedly -deter - 
ministic operator in accordance with the analysis performed by N. E, Kobrinskii and B. A. Trakhtenbrot [4, 5] and 
then study the set of restrictively -deterministic operators which can be distinguished solely by their initial states. 





m or 
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f®=F@ Ww) @=1,2,...). (2) 


We shall assume that the input variable can acquire only two values: 0 and 1, The value s = 1 shall be called 
a “penalty,” and the value s = 0 shall be called “non-penalty.” 


We shall assume further that the output variable f(t) of the automaton can acquire k different values f;,....f,. 
The values of the variable f(t) shall be called the “actions” of the automaton, and we shall say that at the instant t 
the automaton has performed the ath action if f(t)= f, (a =1,..., k). 


The values of the variable ¢(t) are called the “states” of the automaton. We shall assume that the automaton 
A can have m states 9, ..., 9p (m = k), and we shall say that the automaton is in the jth state at the instant t 
if y(t) = yj. The number m shall be called the “memory capacity" of the automaton. The action f, is called an 
action corresponding to the state 9j if fy = F(9j). 


Equation (2) describes the dependence of the automaton actions on its states under these conditions; Eq. (1) 
describes the transitions of the automaton states due to the effect of the input variables. Since the input variable 
s(t) can acquire two values, Eq. (1) describes a pair of representations of the set of states for the finite automaton 
A in itself; one of these representations is specified for s = 1, and the other is specified for s = 0, 


For purposes of our subsequent analysis it will be convenient to describe these representations using a special 
matrix of states for the automaton: A(s) = |lajj(s) Il (i, j= 1,...,m). 


Each row in this matrix contains exactly one element equal to unity for any fixed value of s; the remaining 
elements of this row are zeros, Such matrices are called simple matrices in [2, 3]. The matrix of states A(s) for 
the finite automaton A determines the transitions of states in the following manner: if at the instant t the automaton 
was in the state j;, then at the instant t + 1 the automaton will make the transition to a state gj for which a4; [s(t+ 
+1)J=1. 


From the matrix A(s) it is possible to restore Eq. (1) single-valuedly in such a way that the matrix of states 
and Eq. (2) fully define the finite automaton A.* 


The transitions of states for each value of the input variable s(t) can be conveniently depicted by graphs of 
state.** These graphs are related to the matrix A(s) in the following manner. Each state gj of the automaton A 
is juxtaposed with a vertex i of the graph, and each non-zero matrix element ajj(s) is juxtaposed with an arrow 
dirécted from the vertex i to the vertex j. From the simplicity of the matrix A(s) it follows under these conditions 
that exactly one arrow emanates from each vertex of the graph, For automata which are studied in this paper the 
transitions of states are described by a pair of such graphs (for s = 1 and for s= 0), From the graphs of states it is 
possible to restore the matrix A(s) in a single-valued manner. 


We shall now proceed to study the functioning of automata in media that respond in random fashion to their 
actions, 


We shall state that the automaton A functions in the stationary random medium C = C (py, . . . , Pyd if 
the values of the input variable and the actions of the automaton are related as follows: the action fy (a = 1, ...,k) 
performed by the automaton at the instant t generates the value s = 1 (a penalty) at the instantt + 1 with the proba- 
bility p,,, and the value s = 0 (a non-penalty) with the probability q, = 1—p,. 


Assume that at the instant t the automaton was in the state gj (i= 1, ... ,m),which corresponds to the action 
fay = F(9j)(aj= 1,...-.k). Then the probability pjj that the automaton will make the transition from the state 
#j to the state yj is determined from the formula 


Py = Pa, a; (1) + qa, ais (0) (i, j = i,. . «,™m), (3) 


The matrix P = ||p4jl] of the transition probabilities is stochastic; this can easily be demonstrated by making 
use of the simplicity of the matrix A(s). Thus the functioning of the automaton in a stationary random medium is 
described by a finite Markov chain, Without any substantial limitation of the general nature of our analysis we can 


* For more details cf. [2, 3}. 
** Graphs of this type were used in the paper by V. 1. Shestakov [6], 
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assume that this chain is ergodic [i.e., we can assume that there exist final (limiting) probabilities for the states of 
the automaton in the specified medium which do not depend on its initial state). 


We shall use rj to denote che final probability of the state gj (i= 1,...,m), andog (a=1,..., k) shall 
denote the sum of the final probabilities for such states ¢; which correspond to the action f, [i.e., states for which 
F(%) =f). 


Then the mathematical expectation M = M(A, C) of a penalty for the automaton A in the medium C is ex- 
pressed by the formula 


in 
M (A, C) = >} pace (4) 


a=] 


Henceforth we shall make use of the substitutions 





Pit... +P, 


M max= max (pi, .--, Per), Mmin= min(pi, ..-, Px), Mo = : 


Then it is obvious that 


Mnin << M < Mmax: (5) 


The expediency of the automaton behavior resides in the lowering of the quantity M; it can be defined as the 
closeness of the quantity M to M,,jp. An automaton for which M = Mp is naturally called an automaton which does 
not have expedient behavior. When an automaton with expedient behavior functions in any stationary medium the 
average penalty M meted out to it will be less than Mg (excluding the case where the probabilities of penalties for 
all actions are equal), 


1K pt s=1 a SOF op es SE 
Or 2D 5-0 Operas. 
Fig. 2 Fig. 2 


E le 1. The automaton l,,2 has two states g; and g and two actions f; and §; F() = fy, F(¢2) = &. The 
matrix of states for this automaton will be 


01 10 
A(t)=(49)s 4 @=(93), 
i.e., the states of the automaton will change for a penalty and will remain unchanged for a non-penalty; the graphs 


of states for this automaton are shown in Fig. 1. We shall study the functioning of this automaton in the medium 
C=C (p;, pg). Having formulated the matrix of transitional probabilities in accordance with (3), we obtain 


Pa(hh . 


From this we find the following equations for computing the final probabilities r; and r: 
r= "1 + peta, fs = pir + gare- 
The normalization condition is given by 


™+hm=1. 








Having computed the values r; = p2/(py+ Pz), m = Py/(P; + Pe), we find the mathematical expectation of a 
penalty in accordance with (4): 





2 
M (L2,2, C) = pity + para = a . 
It is not difficult to verify the fact that M(l2,., C) = Mg (equality occurs only for p; = pz). 
Thus the automaton lL,» has expedient behavior. 
Example 2, The automaton len (an automaton with a linear tactic) is a natural generalization of the auto- 
maton Ig; it has 2n states g,..., ¢gn and two different actions f;, &. Under these conditions 


F (1) = F (gs) =~. = F (Gq) = far F (ngs) = F (ngs) = «++ = F Gen) = fr (6) 


The graphs of states for the automaton Ign,z are shown in Fig. 2; the matrix of states for this automaton has the 
elements a4j(8) which are determined according to the following formulas: 


ay; (0) = 1 for i = 2,3,...,m,m + 2,...,2n, i = i— 45 a3; 0) = ayy s ny ( = 1; 
aj (1) = 1 for i=1,2,...n—4, n+ 4,..,2n — fandj =i + 4; (7) 


@n.on (1) = Gonn (1) = 4. 


The remaining elements are zeros. 


We shall compute the mathematical expectation M(Ign,, C) of a penalty for the automaton Ign. which 
functions in the medium C = C (py, ps). 


Having determined the transition probabilities from (3) and (7), we obtain the following system of equations 
for computing the final probabilities mG =1,..., 2n): 


ry = 91% 1 its, Pati = 9? n4i + 9? nies 

ra = Piri + Qirs, Tate Pints + pis wed 

H“es6 6 « whee <p (8) 
= Pile t 91" k+41, Patk ~ Pi?n4x—1 - eras, 

Tn = Pin + Pa?en, Ton = Pion, + Piln 


and the normalization condition is given by 
+m +..-t+hn=1 


Writing the solution in the form rj = a, ~ r tnek = ade k= 1, +++»), we Obtain the following charac- 
teristic equations for determining the eigennumbers A, and A,: 


Gah — 4, + Pp, = 90 (a = 1, 2). 


4 Dorermining A1? » 24 r»f = 1, »P) = ps /ay = cy, »$*) = pe/ce = a, from these equations, we shall seek solutions 
in form 


ry = Arak + By, tyyy = Asay + Bs. 


Making use of the equations for r; and tp +,, we prove that B, = B, = 0; from the equations for tp, mp it follows 


that 
oy oy 
Atfpa~ “Tha: i. 








We now compute the sums 


7 — 1 


n n ay 
5 = Di te = 41 Dat = Ais. 
k=l k=1 





oo Drage = Ae Sob ot, 


From the normalization conditions o, + 6, = 1 and Eq, (9) we find the values of the coefficients A;, A, and 
compute the mathematical] expectation M (lgn,2, C) of a penalty on the basis of formula (4): 


n 


1PA—-% ft 9 
Py Pt py Pa — % 





M (Lyy,gs ©) = (10) 





1Pi—% 1 Ae 
Py Pn pp Ps 





It is important to note that M(I2gnz, C) is a decreasing function of the memory capacity and that for M,,;,, = 
= 1/2 


‘lim M Lon» ©) = Mri (11) 


This relationship means that for a sufficiently large memory capacity n the automaton Ign» performs almost 
exclusively that action for which the probability of a penalty is minimal, Automata which have this property shall 


be called asymptotically optimal, 








k-1)n+2 \ n-Ined ~~. 
(k-1)n+1 K-I)n+1 


Fig. 3 


Example 3. We shall study the automaton Lyn,» which has kn states oeyet » ?kn and k actions f;,..., fy. 
Assume 
F (1) = F (@2) =... = F(9,) = fi, 
F (Pnis) = F (Png) = -- + ¥ (an) = fe 


F (Puapats ) = ¥ (u—iynte) = +--+ = F (nn) = fe. 


The gtaphs of states for this automaton are cited in Fig. 3. It is not difficult to note that for k = 2 the auto- 
maton Lyn, becomes the automaton with a linear tactic described in example 2. 


In the medium C = C(p;, . . . , Py) the mathematical expectation M (Lyn, C) of a penalty for the automaton 
Lyn,k 18 expressed by the formula 
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5 eae { 
Pa—% pr" 

M (Ly, C) = j (12) 
kn,k sy fentet 
Pa 4a p” 








a=] 


The automaton Lyn is also asymptotically optimal, 
Example 4, We shall now describe the automaton G2p 2 which has 2n states and two actions f, and &: 


F (9) = F (2) = .. - = F(Q,) = fi, 
F (®n41) = F (nya) = +--+ = F en) = Io 


The graphs for this automaton are shown in Fig. 4. Such an automaton is naturally called an automaton with 
a hysteresis tactic. 


When this automaton functions in a medium C(p, 1—p) the mathematical expectation of a penalty M(Ggp 3, C) 
is determined from a formula 








a (a*—1)*+ na"? (a —1) P 
M (Gan,2, C) =S 4 (a™ — 1) (a"t? — 1) ana f° ): 


It can be demonstrated that the automation Ggp 2 is asymptotically minimal but that 


M (Den, C)> M (Lon. C), 


so that an automaton with a hysteresis tactic has lower expediency of behavior than an automaton with a linear 
tactic. 


2. The Behavior of Automata in Composite Media 

In the preceding section we studied the behavior of automata in media whose probability properties are in- 
dependent of time. The expediency of the automaton behavior in such media is manifested in a lowering of the 

mathematical expectation of the penalty in those cases when 

7 n-t.n 1 n-f the probability properties of the medium are not known in ad- 
ob... Cy vance. The a priori knowledge of the constant p;, ...,» Pk 
makes it senseless to design an automaton which is capable of 
Riedl o—sow.sn-eomef ) performing a number of different actions depending on its state. 
in Se-tw°E Bf gn dnt a8 atl In fact, knowing the values p;, . . . , Py» it is possible to design 
an automaton which has just one state and performs one action — 


precisely that action for which the probability of a penalty is 
minimal, 


In this section we shall study functioning of finite automata in a medium whose probability properties vary 
in random fashion, Then even if we know the possible probability characteristics of the medium we cannot design 
an automaton with expedient behavior that would have just one action. ; 


We shall assume that the time dependence of the probability properties of the medium is specified in a special 
manner: it is assumed that the medium in which the automaton functions is composed of stationary random media 
whose switching is achieved by a Markov chain. 


Thus we shall study the simplest composite medium K= K(C;, C,, 5) specified by a Markov chain with two 
states C, and C, and a transition probability matrix 


A=(*Z° an 











Fig. 4 
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The state Compegneet te he aationary, random medium Cy = Copp, cee » Bios the state C» corresponds to 
the medium C; = C(p??, eee PR). 


We shall state that the automaton A functions in a composite medium K if at each instant it functions either in 
the medium C, or in the medium C; (i.e., if at each instant its actions and the values of the input variable are related 
in the fashion described above for one of the media), Under these conditions if the automaton was functioning in the 
medium Cg (6 = 1, 2) at the instant t, it follows that at the instant t + 1 it will function in the same medium with 
the probability 1-8, and in the other medium with the probability 5.° 


It is not difficult to prove that the parameter 6 is the average frequency with which the states of the com- 
posite medium are switched. In fact, the average number T of cycles during which the state of the medium remains 
constant can be computed as follows: 


T= > m(i—dns=—. (13) 


m=1 


We shall say that the system S consisting of the automaton and the composite medium is in the state ¥, (8) 
(i= 1,...,m; 6 = 1, 2) at the instant t if at that instant the automaton is in the state gj and the composite medium 
is in the state Cg. The probabilities for transitions between states of the system are described by the cell matrix 
m= lngil G.i=1,...,m): 


re (1) + gay; (O)1(1 — 8) = [PM ais (4) + gag; (0)] at 
sf 


[par (1) + gai (015 = [pags (A) + gay (0)] (1 — 8)" si 


Here pe is the probability of a penalty in the medium Cg for the action fo; = F(g;). Assume that the system 
under study is oy ie Then the final probabilities ry” of the states vO Mfor the system S will not depend on its initial 
states. The mathematical expectation M = M (A, K) of a penalty for the automaton A in the medium K can then be 
given by the formula 


k 
M (A, K) = 2 (po + po, (15) 
a=] ‘ 


where oo 5) ( = 1,...,k; 68 = 1, 2) is the total probability of those states of the system for which the automaton 
performs the action fg, and the composite medium is in the state Cg. 


We shall study the behavior of an automaton with a linear tactic in the composite medium K (Cy, CG, 5 )—for 
example, the automaton Ign 2 in example 2 of the preceding section. 


We shall recall that this automaton has 2n states y;,.... zp and two actions f;, &. It is specified by the 
formulas (6) and (7), and its graph of states is shown in Fig. 2. 


We shall compute the mathematical expectation of a penalty for the automaton Ign 2 under conditions where 
we have assumed for simplicity that 


Ci=C(p,i—p), C2=C (i —p, p). (16) 


In order to determine the final probabilities of the states for the system we have a set of equations analogous 
to (8); however, in this system the probabilities r; (i= 1, . . . , 2n) of the states of the automaton Ign » in a stationary 
medium are replaced by the vectors R, = (r!), rf’) (i= 1, ..., 2n) for the probabilities of the system states; the 
quantities p;, Pe, 43, %@ are replaced by the corresponding matrices. 


* Of course, it is possible to study the composite medium K = K (Cy, ..., Cs, D) that has s states Cy,..., Cs 
corresponding to different stationary random media,and the transition probability matrix D. However, in this paper 
we limit ourselves to the simplest case where the matrix D is symmetrical and s = 2, 








The set of equations is written in the following manner when condition (16) is taken into account: 


Ry = QR: + QRs, Rasy = SRags + SRare 
Ra = SRy'+ QRs, Rags = Ung, + SRars, 
Roth. thy, ovis: _ 
Rg dn + Og. ha oe ae 
where 
Qa (Oy de S(O a ch yg) PI. (18) 


We shall seek the solution in the form 
R, = Red, Rap = RA * (k= 1,.- om) 


where Ry = rf), rf?) is a constant vector. 


In order to determine the eigennumbers \ and the eigenvectors we make use of the characteristic equation 
(QA* — 1B + S) Ro = 0, (19) 


where E is a unit matrix of the second order. 


From the equation Det (Qd”—AE +S) = 0 we find four values of the eigennumbers; A = Az = 1, and Ag and A, 
which are roots of the quadratic equation 


(a + 1)* (1 —38) 
fala 22( : a— Wy —1) +1 =0. (20) 





Then, finding the eigenvectors from (19), we write the solution of the system in the form 


R, = ARP) + BRE + crR® + Dre, 
where ! 


RY = (1,1); 

R® = (n (1 + 0) > + (1 — av, te (1 +) p55 — (1 —a); 
RP = (— Ad (e+ 1) 5g, ad (203 — Ge + 1) ng + 20) 
RO = (— at @ + yp. att (on — (a + 1) OP ng + 28). 


Then, we determine the coefficients A, B, C, D from the equations for R, and Rp, 4 1 and write the solutions 
for Ry = (re, 1°) (k= 1... m): 





i an ate (z —a) (1 —az) + 2°™—*+1 (2x —a)(a—1) +-2* (1 —az)(a—1) +(2x—a) (4 — az) 
k a (1 — 2) [x*" (x — a) — (1 —az)] a 





png et + Oe) + ee 
t a (i —2) [2" (@—a) —(i —az)] 


where x is any root of Eq. (20); d is a normalizing multiplier. The solutions Rp +, = (¢ ) ». @) k= 1,...,0) 
are found from the equations 


Om rh, hy = MD, dai 
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Making use of expression (15), we arrive at the following for the mathematical expectation of a penalty for 
the automaton Ign» with a linear tactic that functions in the medium K: 














4s h ny —1 
M (Le: K) 4 SS : 
2 2 2 = (a + 1)¥chny + (a—1)'eth J shny o) 
where 
— (ita? 1—s = ,- 
ae te oe oO eee 


Figure 5 shows the curves which characterize the dependence of M on n for various values of 5 and p = 0,33, 
Figure 6 shows the same graph for a fixed 5 = 0,01 and various values of p. 














Fig. 5 Fig. 6 


. It is not difficult to prove that M(Ign2, K) reaches the minimum m for a certain value of the memory capacity 
= ng; the cases n = 0 and n = » produce no advantage in comparison to an automaton which does not have ex- 
edient behavior. 


The presence of this minimum can be explained on the basis of the fact that for a memory capacity which is 
too small the information on that state of the composite medium in which the automaton functions is not adequately 
used; for an excessive increase of n there occurs, as it were, an averaging of the statistical properties of the two 
states of the composite medium. It is natural therefore that as5 decreases we observe an increase in the value of 
Ng and the value of m decreases.* For 5 -» 0, formula (23) becomes (10) so that ng », m +M,j,, = min (p, 
1—p). As & increases, the reverse process occurs. Thus, for 


1—28 (co + 1) 
6(1— 4) S a 





the minimal mathematical expectation of the penalty is reached for n= 1, 


Formula (21) permits us to construct automata with a linear tactic that have the most expedient behavior in 
a specified composite medium. 


In choosing the memory capacity for such automata it is useful to employ a table which cites the values of 
ng and d= "4 -m for several values of p and 5. In each “cell” of the table we write a pair of numbers: the first 
number is no, and the second number is d. 


* Note that,based on the meaning of the parameter 5 [formula (15)}a decrease in 5 is equivalent to increasing the 
rate at which the automaton operates. 
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Values of 6 























P 0.001 0.010 0.032 0.100 0.320 0.450 
0.10 3; 0.396 2; 0.372 2; 0.336 1; 0.256 1; 0.115 1; 0.032 
0.20 5; 0.294 3; 0.266 2, 0.223 2; 0.157° | 41; 0.065 1; 0.018 
0.25 6; 0.244 4; 0.212 3; 0.172 2, 0.116 1; 0.045 1; 0.012 
0.33 8; 0.158 5; 0.125 3; 0.091 2; 0.055 1; 0.020 1, 0.006 
0.40 11, 0.089 6; 0.056 4; 0.037 2; 0.020 1; 0.007 1; 0.002 | 
0.45 16; 0.036 7; 0.017 4; 0.010 2; 0.005 1; 0.002 1; 0.001 | 














Note also that in our case My* ¥2, and the difference d='/2 —m can serve as a measure for distinguishing 
between the random media C, and C, for a finite automaton with a linear tactic. 


In conclusion I express my deep appreciation to I. M, Gel'fand, D, S, Lebedev, and O, B, Lupanovy for their 
interest and attention to this paper, and to B, D. Efremov for his great assistance in performing the computations. 
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THE SYNTHESIS OF CONTACT NETWORKS 
WITH "REAL" CONTACTS*® 


V. N. Roginskii 
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pp. 1355-1359, October, 1961 

Original article submitted February 24, 1961 


The paper presents a method for synthesizing contact networks which do not interrupt operation 
during transients, A mathematical device based on three-digit logic is cited. 


In the process involved in the operation of a multicycle relay network the changes in the states of an individual 
relay do not occur instantaneously, and in the "transient period” contact circuits may be disrupted due to a non- 
simultaneous closing and opening of individual contacts of the same relay. An analogous picture holds for the case 
where several relays change their states during transition between one state of the unit to another (the so-called 
"competition between relays"). As we know [1-3] the existing contact network algebra, which is based on two- 
digit Boolean algebra, treats contact networks in static states by eliminating transient periods from the analysis. 
Recently a number of papers have appeared [4-8] which study the operation of contact networks during transient 
periods, These papers have chiefly been devoted to analysis, During such periods Gr, Moisil [6] proposed the use of 
three-digit Lukas'evich logic. 


In this paper we demonstrate that this device can be used to synthesize contact networks which do not produce 
violations during transients. 


The static states of the relay, as usual, correspond to 0 and 1; the transient period corresponds to 1/2. The 
relay contacts at each specified instant can be either open (zero) or closed (1). The states of normally open (a) and 
normally closed (4) contacts of the relay A are defined uniquely in static states, but it is impossible to indicate their 
states during transient periods because the contact may close or open either at the beginning or at the end of the 
period. Such an indefinite state of the contacts a and @ is denoted by 1/2. One normally open contact plus one 








State of relay | 0 1/2 1 State of relay | 0 172} 1 
a 0 1/2 1 a” 0 1 1 
a 1 1/2 0 a" 1 1 0 
a’ 0 0 1 a 0 1 0 
a’ 1 0 0 x 1 0 1 























normally closed contact form a switching contact. We shall study two types of switching contacts [9]: the "throw- 
over” contact (Fig. 1a) for which the normally closed contact first opens when the relay operates, and the normally 
open contact only closes after this operation has been completed (i.e., there is a time interval during which both 
contacts are open); and a “transitional” contact (Fig. 1b) for which closure occurs first and is then followed by open- 
ing (i.e., there is a period during which both contacts are closed). These contacts are respectively denoted by a'—" 
and a"—a", Moreover, we shall study "temporary -make" (Fig. 1c) and "temporary-break” (Fig. 1d) contacts which 
close or open the circuit only during transient periods. These contacts are denoted by # and ¥, respectively. The 
states for all the contacts during different periods of relay operation are characterized by Table 1°* (Fig. 2). 


: Paper read at the Fourth All-Union Mathematical Conference on July 7, 1961.. 

** Note that a’, 3, a" and a” correspond to the Lukas'evich modalities [6] va, na, a, and ya. The symbol & 
for a single relay coincides with the operator "Da" [10], but for multiple-element units they do not correspond to 
each other. 














In order to transform networks with such contacts we shall make use of the three-digit Lukas’ evich logic [6] 
in which the operations of addition (+), multiplication (-), and inversion (N) are defined as follows: 





x+y 0 "le f xy 0 "ly f 
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The letters x, y, and z may correspond to any contact symbols in Table 1 which are treated as variables, Addi- 
jual tion and multiplication are interpreted as parallel and series connection of contacts, and inversion is interpreted as 
finding the circuit with the “reciprocal” structural admittance, 


e In this algebra, which is a distributive structure, the following relationships are valid: 
z+2=2, w= 2, (1) 
. z+y =~y+2, zy = yz, (2) 
z+ (y+ 2) = (e+ y) +2, x (yz) = (zy)z, (3) 
7 z(y+2) =2y +22, 2+ ys = (z+ y) (x + 2), (4) 
a zr+ay=2, c(z+y) = 2, (5) 
heir z+0=2, 20=0, (6) 
' z+i=+1, z-1 =z, (7) 
NNz = 2, (8) 
N (x + y) = Na-Ny, N (a-y) = Nz + Ny. (9) 


Based on these definitions, we establish the fact that the following relationships hold between functions corre - 
sponding to different contacts: 



































x aja a’ |a’la"\la"|\ ala (10) 
ow o - - v a 
ally Nz\| a|a|a" \|a"\a@’\a’\a ia 
h 
spen- é=a’"-@’,i=a'+d, (11) 
"7 
fhich a+a4+4=1, a-4-4=0, (12) 
he 

a’"=a+4, @#@=4+4 4, (13) 

$ a’ =a-d, @ =4-d. (14) 
to 
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Without treating any other of the possible transformations [6] we shall use the algebra cited above to formulate 
contact networks which include “make” (a), “break” (_), “throw-over"(a’, a") and "transitional" (a", *) contacts. 
If contacts a or @ are present in the structural formula as a result of synthesis, then this means that no requirements 
are imposed on them with respect to their operation in transient periods. 


We shall use the switching table as the basis for synthesis. First we shall study a table which allows transitions 
only into neighboring states (i.e2, in each transition from one cycle to the next a change in the state of just one relay 
is allewed). 


By analogy with the synthesis of networks for static states [1-3] we shall juxtapose the transient period during 
which the state of a certain relay A changes with a circuit that contains a series temporary-make contact # and 
those contacts of the remaining relays that correspond to their states in the previous and next cycles. Thus, for 
example, if a network consisting of three relays A, B, and C makes the transition from the state 001 to the state 101, 
then we shall juxtapose this transition (as well as the reverse transition) with the circuit abc,which will be closed only 
for certain transitions, 
‘ ! 
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In order to obtain a circuit which will not be broken during the transient periods it is necessary to sum not only 
the gransfers corresponding to the static network states in which the specified circuit must be closed, but also the ex- 
pressions corresponding to the transient periods. Under these conditions the general solution [3] can include as con- 
ditional components expressions corresponding to all transitions which do not appear during the process of operation 
of the given network (“unused transitions"), In this case the latter include all double (of the type 6c), triple, etc. 
transitions, We can assume that all transitions which directly precede or follow the period during which the circuit 
must be closed are conditional transitions, 


As an example we shall study a network with three relays which operate according to the Gray code and in which 
it is necessary to obtain a circuit that is closed during the cycles 3 to 6, The operation of this network is described 


by the following switching table: 


No, of cycles 012345678 
A 011001106 
B 0011441000 
Cc 00001414110 
f 

No, of states 013267540 


If we neglect the unused transitions, then for the circuit in f the structural admittance will be 
t = abe + abe + abe + abc + abe + abc + abe. (15) 


Performing transformations by factoring and making use of the relationship (12), we obtain 
f = ab (e+ ¢+ ¢) + abe + ac(b+ 6 + b) = ab + abe + ac. (16) 





hich 





The presence of the term &bc indicates that we must create a circuit in the network which closes during the 


transient period corresponding to operation of the relay A. In order to eliminate 2 we combine the third, fourth, 
and fifth terms in (15); as a result we obtain 


f = ab + be + ae. (17) 
The corresponding network (Fig. 3a) does not produce breaks irrespective of the sequence in which the contacts 
operate during the transient period, 


In this example the unused transitions will be all the double and triple transitions as well as the transitions 
000-010, 100-101, 110-111 and 011-001, When we take into account the “extreme” transitions 110-010 and 101-001 
the general solution for the circuit f will be 


abe + abe + abe + abe + abe + abe + 

abc, abe «= abe «=, abe =O, be 
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Discarding all of the conditional terms which do not result in simplifications we obtain 4: 
f=ab+ac+a—. (19) 
Making use of relationship (13), we transform formula (19) so as to eliminate 
f=ab+ac+a(b+c)=(@+a)b+(a+a)c=a'b + a"c. (20) 


This expression corresponds to a network with a transitional contact belonging to the relay A (Fig. 3b). 


Inversion in accordance with expressions (8)-(10) leads to networks that have the reciprocal structural admit- 
tance not only in their static states but also during their transient periods, For the networks 3a and 3b indicated 
above we obtain 


N (ab + ac + be) = (a +b) (a +¢)(b + ©), 
N (a"b + a"c) = (a’ + b) (a’ +0). (21) 


Expressions (21) correspond to the networks in Figs. 3c and 3d. 


The method cited above can also be used to obtain a circuit that closes only during a certain transition. Thus, 


for example, if it is required to obtain a circuit g at the instant when the relay A operates for the first time we can 
write 


a-=- ——_-— 


g=abe = a*a"be. (22) 


Now we shall study the case where the transition is not to a neighboring state (i.e., the case where the succes- 
sive cycles differ in the states of several relays). Under these conditions we shall assume that the sequence of the 
lag in the change of state is identical for all the relays. During the transition process it is possible for all combina - 
tions of operating and non-operating relay states to appear for relays which change their states. 

Thus, for transition of a device with two relays A and B from the state 00 (with the number i = 0) to the state 
11 (j = 3) the transition may take three forms: 00-11 (both relays change their states simultaneously), 00-01-11, 
and 00-10-11, In order to obtain a circuit (we shall denote it by 24) which is closed during the entire transition 
time it is necessary to sum the expressions corresponding to all of the states which may appear; i.e., in this example 


24j= ab + ab + ab + ab+ ab +ab+ ab = 
a(6+6+6)4+6@+ a4+4+a)+ab+ab- 


a -— = (23) 
=a+6+4ab+ab. 

















have the following form for the transition from the state i (which corresponds to the transfer k;) to the state j (with 
a transfer kj): 


during the process involved in synthesizing relay networks; it also makes it possible to determine the necessity of 
using “throw-over” or “transitional” contacts to assure normal operation of a network during the transient period. 


1, 
2. 
3. 


9. 
10, 





By analogous reasoning for the case where s relays Ay, Az, . .. , A, change their states the circuit 24; will 


a 
y= Diar+ > hn. (24) 
Ist 


li,j 


The method cited above permits us to take into account not only the static states but also the transient periods 
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The paper studies the problem of finding the optimal technological mode of operation for 
complex production processes which do not have a mathematical description. The block 
diagram is cited for a learning automaton of the tabular type which accumulates and 
evaluates the information arriving from the object in its natural mode of operation, and 
then develops recommendations to the operator on the optimum mode for the process. An 
approximate computation of the learning time is made. 


1. Formulating the Control Problem 

In the chemical, metallurgical and cement industries there exist very complex processes which can be charac- 
terized by the following features: 1) a large number of parameters which affect the course of the process; 2) the 
presence of uncontrolled external perturbations which affect the process in such a way as to possibly alter the time 
characteristics of the object; 3) the absence of any information on the form of the function that interrelates the 
process parameters; in particular, this function may be discontinuous, may have several local extremums, or may 
have no extremums at all. 





We shall introduce the following classification of object parameters (Fig, 1): 


1. Thesetof primary process parameters is denoted by the vector x ae Xm . These are the con- 
trolled parameters whose magnitudes depend on external factors and cannot be altered by the operator. The primary 
parameters include the indices for the input components of the process, certain 
conditions governing the flow of the process, etc. 














2. The set of secondary process parameters which characterize the quality 
| of the output product fs denoted by the vector k{ ky, ke, ..., ky}. 
t——-4Object [=F We shall assume that there is a certain region K’, where the values of the 
| vector k satisfy specified requirements governing the quality of the output product. 

z 3. The index of the process effectiveness is denoted by Il. The index can 
Fig. 1 be measured at the output of the object or it can be computed in an arithmetic 


unit from the known values of the other parameters. 


4. The set of controlling parameters which the operator can affect by altering the process in one direction or 
the other is denoted by the vector z { 2, %,..., 2}. 


5. The set of uncontrolled perturbations is denoted by the vector y {yy Ye,+-+»Yn}+ We call those per- 
turbations uncontrolled for which the magnitudes cannot be measured for some particular reason. In particular, un- 
controlled perturbations include variations in the characteristics of the technological equipment due to aging, un- 
controlled quality indices of the input components, etc. 


The problem of controlling the object reduces to the following: for specified x and y it is necessary to deter- 
mine the value of the vector z corresponding to the maximum value of the index II; i.e., it is necessary to find the 
maximum for the function 


II = Fy (2) 


for the condition k €K° and z€Tz, where I, is the region of allowed values of z. 








In the case where k is outside the region K° it is necessary to determine the value of the vector Ag which re- 
turns k to the region K®, 


In view of the absence of complete information on the process and on the equations interrelating the para- 
meters of the object, the operator relies to a considerable extent on his experience and intuition in establishing the 
technological anode, The effect of subjective factors on the control process often makes the mode of operation far 
from optimal; this leads to a lowering of the efficiency and the incomplete utilization of the technological equip- 
ment, etc. 


2. Certain Concepts Concerning the Realization of the Automatic Unit 

The problem of designing an automatic unit for achieving an orderly control and finding the optimal techno- 
logical mode for complex objects is associated with great difficulties. A necessary condition for realizing such a 
device is the presence of a fleet of primary instruments which provide ob- 
jective information on the existing factors that determine the technological 
mode with a sufficient frequency. A preliminary investigation of the process 
must establish the minimum number of parameters which must be in- 
troduced into the control system. For a large number of parameters the 
introduction of information into the automatic control unit should be 
realized from a scanning monitor system which scans the process para- 
meters in sequence. For purposes of introducing them into the digital con- 
trol unit the analog quantities from the outputs of the primary instruments 
are quantized with respect to level and with respect to time. Time quan- 
tization of the input quantities is achieved during the time (tq) required 
for a cycle of operation of the scanning monitor system, 








Fig. 2 


In Fig. 2 curve a reflects the actual variation of the parameter x;, 
and the line b represents the transformed dependence of the parameter 
x4 On time which was obtained at the output of the scanning monitor system, Here o is the magnitude of the level 
quantization, and tq is the time required for an operating cycle of the system which derives the information. 


In studying the problem involving the principle according to which the automatic unit is designed it is neces- 
sary to take the following factors into account: 


1) the controlled object has a lag t, determined by the spread between the points at which the primary and 
secondary information is extracted; here the quantity ty may reach several tens of minutes; 


2) all the primary parameters are measured at the object input; the secondary parameters are measured at 
the output; the controlling parameters may affect the process at any instant from 0 to t,. 


The device which makes it possible to determine the optimum technological mode of operation for a complex 
object of the type under study may be a system that performs storage and statistical manipulation of the information 
arriving from the object when it is in a natural mode of operation. The control algorithm is found by the system as 
the result of manipulating information from the object. The control algorithm is perfected to the extent that the 
quantity of manipulated information is increased. 


It is evident that at each instant the effectiveness of the control is determined by the degree to which the con- 
trol algorithm determined by the automaton corresponds to the optimum control method that applies to the state of 
the object at the instant under study. 


When this correspondence is violated it is necessary to adjust the control algorithm; here the algorithm for the 
actual setting is determined as a result of storing and evaluating information derived from previous settings. A con- 
trol system of this type can be called a learning system, since the control algorithm which it formulates is perfected 
with the passage of time to the degree that information arriving from the object is stored and processed. 


3. The Block Diagram for a Learning Automaton of the Tabular Type 
We shall study the simplest learning system —a learning automaton formulated according to the tabular meth- 
od [1], 


The block diagram for an automaton of the tabular type is shown in Fig. 3. The extraction of information 
from the object is performed by the scanning monitor system. In extracting the information the system takes into 
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account the lag and the spread between the points at which the control signals are applied. Information is fed into 
the automaton at time intervals equal to a cycle of operation tq for the scanning monitor system, The automaton 
contains two basic logic blocks and an arithmetic unit AU which is designed for computing the index 1 and for 
calculating the difference Az=2zj;—z;-, and the sum 
z + Az, where z; and 2; -, are the values of the vector 
: Object K z at the ith and (i—1)th instants of information extraction, 
respectively, The memory of the block 1 is a certain table 
where it is possible to write up to m different vectors (z;.) 
pe ee % from the addresses x;, whereby anew vector is written in 
the cell with the address x; for the condition 
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- morro F (2;, x, t) = ky @R®. 
Zz f 
, 2 Ae Determining the mode of operation is achieved 
Az cyclically by choosing the vector 2; corresponding to the 
. maximum value of the index II from the cell with the 
z+Az address x;. Thus the mode of operation is found before the 


information on the secondary parameters is extracted; this 

is necessary when the values ty are appreciable, If the 
appearance of uncontrolled inputs causes a violation of the 
process (if in the mode found by the automaton the vector 

k is shifted outside the region K°), then the contents of the memory for block 1 must be changed, This change is 
performed by means of the block 2,whose memory contains the table in which the vectors Az are written from the 
addresses ky (the vectors Az represent quantities by which it is necessary to vary the controlling parameters in order 
to return the process into the region K°), 


Fig. 3 


We shall study the problem involving the organization of the memory for block 2, 


We shall assume that definite sets of values of the vector y correspond to a definite probability interrelation- 
ship between the magnitudes of the performance deviations and the magnitudes of the variations in the controlling 
parameters. This interrelationship, for example, may be reflected in the probability sequence governing the use of 
the different vectors Az or in the probability frequency governing the use of the vectors Az for a specified Ak. 


In the latter case it is possible to obtain a certain distribution of the vectors 4 z which for specified conditions 
makes it possible to realize search in the minimum time, 


In order to realize different distributions of a finite number of vectors Az we can proceed as follows, We 
shall write the vectors Az with frequencies proportional to the required probability in the cell with the address Akj. 
When any number can be chosen from the cell Ak; with an equal probability, the probability for choosing each 
Az ; is proportional to the number of Az; which have been written; therefore it is proportional to the probability of 
Az; in the distribution which is reproduced. 


If a variation of the vector y causes a variation of the frequency with which the specified vector Az; is success- 
fully used when a performance deviation Ak; occurs, then this leads to a variation in the probability with which the 
vector Az; is chosen from the cell Ak;; i.e., it leads to a variation in the distribution of the probabilities with which 
various Az are chosen, 


A change in the contents of the memory for block 2 changes the algorithm according to which the search for 
the region K° is performed; a change in the contents of the memory for block 1 changes the control algorithm which 
is formulated for the process, The adaptation of the contents of the memories for both blocks to the state of the ob- 
ject which exists at a specified instant constitutes the learning process in the automaton under study. 


The value of the vector representing the controlling parameters z ; which is found by the automaton is given 
to the operator in the form of a recommendation with respect to the mode of the process. The degree 
to which the automaton has “learned” can be defined as the ratio between the number of recommendations that have 
been successfully used by the operator during a specified time interval and the over-all number of recommendations 
produced. 





In principle the automaton under study can also operate without any participation of the operator. For this 
Purpose it is necessary to write n vectors Az of equal probability for each address in the memory of block 2. Since 
the automaton loop is closed by the object, the automaton sets the memory of block 2 and fills the memory of 
block 1. 


4. The Approximate Estimation of the Learning Time for the Automaton 

We shall perform an approximate quantitative estimate of the required learning time for an automaton into 
which information is introduced at time intervals tg. Since the components of the vector x which are introduced into 
the automaton are quantized with respect to level, it follows that the vector x can acquire N different values. Assume 
that the appearance of the uncontrolled vector yj; has made it necessary to change all the recommendations written 
in the memory of block 1 of the automaton; here it is necessary to change m, recommendations for each address 
X,. We assume that k= 1,2,...,N. Assume further that the appearance of any of the N addresses at the discrete 
instants of time itg(i = 1, 2,..., S) is an independent random event with the probability p= 1/N, This assump- 
tion is valid for tq >> MT,,where MT, is the mathematical expectation of the random quantity T, that represents 
the time interval between the appearance of two different values of the vector x. 
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We shall compute the probability Pyys for the appearance of each address x}, no less than m, times during S 
cycles: 





nN +My 
ogtmogtine CEM aa, 
Nk 
ze ~ (1) 


Here N° is the number of all possible sequences of the values of the vector x at S instants of time itg. 


The numerator of expression (1) is the number of all possible sequences of the values of the vector xin which 


each value x, is encountered no less than m,, times; CP1* ™1 is the number of combinations of S taken n, + m, at a time; 
N 


(mg, Mg, . - - » My» + » » My) is any combination of the numbers 0, 1, 2,..., (S— Xm), for which the condition 
k=1 


N N 
> hy = S — > my 
k=1 k=1 


is valid, 
Simple transformations of the formula (1) yield the following expression for the probability P),<: 


S! i 
Pus = ne (ny + my)! (ng + ma)! .. «(my + m,)!... (Ay + my) * (2) 





Figure 4 shows the dependence of Pjyjg on S for N= 5 and my = m=...= my = 1. By specifying the probabiliy 
Png it is possible to use formula (2) to determine the number S of cycles tq required for assuring the appearance of 
each of the N addresses x}, no less than my, times with the probability Pys: 
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N 
S = ¥,(Pys, N, > ms) (3) 
k=1 
In using Eq. (3) it is assumed that during the learning process the choice of the correct value of the vector 
Az by the operator or by the automaton is achieved with a probability equal to unity at each instant of it,. In fact, 
what occurs is the process involving the setting of the memory in block 2 or the process involving the adaptation of 
the operator himself; as a result the choice of the correct value of the vector Az occurs with a probability close to 


unity only after a certain period of time 5. Assume that the process of setting the memory in block 2 is exponential 
in nature (Fig. 5) and is described by the equation 


_ 3st 
p(t)=i-e *, (4) 


where p(t) is the probability for the proper choice of the vector Az; 6 is the time required for setting the memory 
in block 2. 


We shall compute the required number of cycles tg for learning when the setting of the memory in block 2 is 
taken into account. 


We shall determine the probability Psq that learning will require A cycles tq; i.e., we shall determine the 
probability that for A cycles the vector Az will be correctly chosen exactly S times, The solution of this type of 
problem has been treated in [2]. The unknown probability is computed from the formula 


Ch 
Psa = >) Ri, 
i=1 
where 
R, = p (Arty) p (Aa ta). - - P (Aste) [1 — p (Asyate)) . . 11 — p (Aat,)), 
Ch is the number of combinations in A taken S at a time, p(Ajtq) is the probability that the correct value of the 
vector Az will be chosen at the instant Ajtq; [Ay, Ag. . . - » Ajy » + » » AA) is any transposition of the numbers 1, 
2, 3, eoeere A, 
The probability Pa that at least A cycles tg will be required for learning is determined from the formula 
a—s eat? 
P,= > > R;. (5) 
j=0 i=1 


Figure 6 shows the dependence of P, on A for S = 20 and 5 = 0.3 Stg. 


Specifying the probabilities Pa, it is possible to determine the required number A of learning cycles using 
formula (5): 


A = Fs [Pa, S, p (t)). (6) 


A simultaneous study of expressions (3) and (6) permits us to compute the over-all learning time with the 
probability P = P)jsP4 when the time taken to set the memory of block 2 is taken into account: 


Ty = F,[ Pa, F, (Pw, N, > mx), P(t) | ta. (7) 


k=l 


In the example treated above, when N = 5 and m; = mg =... = My = 1, and the probabilities Pyyc = 0.88 and 
Pa = 0.997 are specified, it is possible to state that Ty, = 23 tq with the probability P = 0.87, 


For successful operation of the automaton it is necessary to satisfy the inequality 








Th <-% Twnin for B > 1, (8) 


where Tymin is the minimum time between substantial variations of the vector representing the uncontrolled inputs 
that is manifested with a sufficient degree of reliability. 


A learning automaton of the tabular type which is designed according to the scheme described above can be 
used both for ordering the process involved in the control of operating objects and for formulating technological 
operating -mode charts for objects that are being made operational. 
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The paper describes a white-noise generator which is designed for use in statistical investigations 
of infralow-frequency objects. 

The principle governing the operation of the generator is based on using fluctuations in a 
stream of particles from an ultralow-power source of radioactive radiation. 


The use of statistical methods to investigate automatic control systems is in a number of cases associated with 
the use of noise-signal generators. 


The usual requirement applied to noise generators is uniformity of the spectral noise power in a specified fre- 
quency band and a definite law governing the distribution of the signal amplitudes. An important requirement is a 
stationary noise signal and the absence of any parasitic regular components. 


We are thoroughly acquainted with noise generators that are used in the audiofrequency and radiofrequency 
ranges and use a gas-discharge device (for example, a thyratron) or a vacuum device (for example, a noise diode of 
the type "DSh-2,” etc.) as the primary source of the noise emf. The spectrum for the noise emf generated by such 
sources is usually sufficiently uniform over a wide frequency band. Such sources are low-power devices. In a num- 
ber of cases it is necessary to use noise-signal generators operating in the ultralow-frequency range when statistical 
investigations of automatic control systems are made. The generation of noise emfs in this range by means of the 
primary sources indicated above is an extremely complex engineering problem, since a portion of the energy from 
the noise-signal source that is confined to the required low-frequency band represents an infinitesimal part of the 
overall energy in the noise signal produced by the primary source. In this case the use of amplification of the order 
of 10*—10° produces distortion of the energy distribution of the original signal over the spectrum; it also introduces 
an additional nonstationary factor and parasitic regular components. 


It is especially difficult to obtain a stationary noise signal at a sufficiently high level with the required statis- 
tical characteristics in the infralow-frequency range 1 - 10~ to 1 cps. It should be noted that this range is charac- 
teristic for a large number of automatic control objects, in particular for the majority of industrial objects. 


The literature describes several generator networks which produce low-frequency noise signals [1-3]. The 
essence of the operation of these generators resides in using a band-pass filter to isolate a sufficiently uniform sector 
of the noise generator spectrum (for example, the spectrum produced by a gas-discharge device) at a frequency of the 
order of 1 to 10 kc and then converting this isolated noise into low-frequency noise by heterodyning. The structure 
of a generator based on such a method is shown in the block diagram of Fig. 1,which is described in [3]. Such meth- 
ods do not eliminate the necessity of employing large gains; this, as we have pointed out above, is undesirable, In 
[2], for example, the degree to which the noise signal is stationary is improved by making use of a so-called “noise 
AGC* circuit which is sometimes used in radio receivers. The use of AGC greatly complicates the noise generator 
circuit. 


The generators described in [1-3] produce noise with a uniform distribution over the spectrum in a band of 
the order of 0 to 30 cps. In order to obtain noise in a band covering lower frequencies (in a band ranging from 0 to 
3 cps) the method [4] based on using fluctuations of an electrical resistance formed by overflowing metal spheres be- 
tween the conducting axis of a rotating drum and its conducting cylindrical surface has been proposed. The incon- 
venience of using such a generator resides in the presence of mechanical and electromechanical elements, as well 








as in the dependence of the statistical characteristics of the generated noise on the structural features of the device 
or their variations. 


The method described below for obtaining a noise signal in the infralow- frequency range is based on the use of 
a radioactive-particle counter and permits us to use a simple circuit to obtain a noise signal with an adequate level 
(dispersion) and the uniform spectral density in the required frequency band, The law governing the amplitude dis- 
tribution is practically the same as the normal law. 
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Fig. 1. Block diagram for obtaining low-frequency noise by the heterodyning 
method: 1) Audiofrequency noise generator; 2) band-pass filter; 3) mixer; 
4) heterodyne; 5) low-pass filter; 6) output. 


In this case we use the fluctuations in the radioactive decay process that occurs in an ultralow-power source 
as the primary noise signal.” 


As we know, such processes are characterized by the Poisson distribution law 
P(k)= te, (1) 
The quantity At which appears in expression (1) is the average number of particles produced during the time 


interval t, 


The pulses obtained at the output of the gas-discharge counter can be converted by means of a multivibrator 
with two stable states into a so-called random telegraph signal; this signal is a voltage that acquires one of two 
possible values (u or —u) in random fashion (cf. Fig. 2). 


From expression (1) it is not difficult to obtain the expression for the 














- spectral density G(w) of a random telegraph signai in terms of the cor- 
" relation function for the process (cf., for example, [5, 6]): 
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Fig. 2. Random telegraph signal. 1+(+) 


The graph for the spectral density G(w) is shown in Fig. 3a, From expression (2) and the graph in Fig, 3a it 
is evident that in the low-frequency range it is possible to obtain a sufficiently uniform energy distribution. The 


uniformity of the spectral density G(w) can easily be seen to depend on the relationship between \ and w. By 
varying it is possible to obtain the required uniformity of the curve G(w) in a pass-band with the specified upper 
frequency limit w,;. In the majority of practical cases it is possible to use the simplified expression for the spectral 


density in the low-frequency range: 
Gw)=s. @) 


* For example, such a source may consist of a granule of glowing material which is used to coat the numbers on the 
scales of instruments, clocks, etc. 
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A precise estimate of the required value of ) in the frequency range with the specified upper frequency limit 
whi and a definite degree of nonuniformity n can be obtained from the relationship 


V i 1 
A =pj = Yeh woes WE , (4) 


where w pj is the upper frequency limit for the required noise band, n is the degree of nonuniformity of the spectral 
density : 


G (0) — GWyi ) 
n= CW) , (5) 





Expression (4) makes it possible to estimate the value of the upper frequency limit w},; for a specified degree 
of nonuniformity n and an experimentally determined value of A. The method for determining \ is indicated below. 
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Fig. 3. Graph for the spectral density G(w) and the probability den- 
sity f(x) of the generated noise, 


The block diagram of the generator is shown in Fig. 4, As the counter 1 we can use one of the types of stan- 
dard halogen radioactive -particle counters (for example, the "STS-5" model). The multivibrator 2 is used to con- 
vert the chaotic pulse train arriving from the output of the counter 1 into a random telegraph signal (Fig. 2). The 
time constant for the low-pass filter 3 is chosen from the condition 


>. (6) 


where \ is the mean pulse frequency required to assure the specified n in a pass-band with the upper frequency limit 
Whi. The buffer stage 4 is designed to eliminate the effect of the load on the mode of operation of the generator. 


The measurement of 2 is performed by switching the multivibrator input (and the input of the integrating cir- 
cuit 6 with the dc meter 7) to a calibrated pulse generator 5. The value of the measured ) corresponds to the reading 
on the generator scale (in cps) for which the dc component of the voltage at the integrator output (measured by the 
meter 7) is equal to the dc component obtained when the integrator is connected to the output of the particle counter 
1. A variation of is easily achieved by controlling the position of the 6 -particle source relative to the counter 
frame. 


The generator circuit contains two dual triodes (of the bantam or ultraminiature type). One of these forms 
the multivibrator with two stable states that is triggered by the pulses from the gas-discharge counter. One-half of 
the other dual triode forms the buffer stage (a cathode follower). The second half of this tube is used as a charging 
diode in the integrating circuit 6. When it is necessary to control or subtract the de component we use a divider 
which is fed from a common plate voltage source. 


The supply source for the generator consists of the simplest rectifier circuit,which is used for supplying the 
plate circuits, a step-down winding for supplying the filament circuits, and a voltage doubler rectifier circuit for 
supplying the counter. The generator operates from a 220/127 v, 50 cps line. 


The level of the noise signal (mean-square deviation) can be controlled over the range 0 to 6 v by varying the 
amplitude of the pulse voltage taken from the output of the multivibrator tube. 


In performing an experimental investigation of the generator operation we shall devote our major attention to 
clarifying the law governing the noise amplitude distribution. Based on an analysis of the principle governing the 
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Fig. 6 


operation of the generator, it is not difficult to see that the amplitude distribution of the signal before the low-pass 
filter differs from normal distribution. As we know, a low-pass filter whose pass-band is appreciably less than the 
input noise band has a strong normalizing effect on the amplitude distribution at the output (cf. ,for example, [7]). 


In order to clarify the law governing the amplitude distribution of the output signal we investigated experi - 
mental data and formulated a statistical distribution function [cf. [5]). The graph for the statistical distribution 
function is shown in Fig. 5. The same figure shows the curve corresponding to a normal distribution law with a dis- 
persion equal to the statistical dispersion of the investigated noise. The closeness of the curves proves that the nor- 
malizing effect of the low-pass filter is adequate for practical purposes. When necessary, inequality (8) can be 
strengthened. Thus, the probability density of the generated noise can be approximated with a sufficient accuracy 
by the expression for the normal (Gaussian) distribution law 


—(x—m,)? 
i 
f(x) = Vin ¢ ™ (7) 


where m, is the mathematical expectation; 0, is the mean-square deviation. 
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The graph of the function f(x) is shown in Fig. 3b. As an example of the resulting signal Fig. 6 shows one of 
its realizations. 


The resulting data makes it possible to draw the following conclusions concerning the proposed method of ob- 
taining infralow-frequency noise. 


1, The generated noise signal has a uniform spectral density G(w) which is determined by the expression (3) 
in the infralow-frequency range beginning at zero frequency. 


2. The probability distribution for the noise amplitude can be approximated by the normal distribution law 
with a sufficient accuracy for practical application. 


3. It follows that the noise signal is stationary from the properties of radioactive decay. As we know, any 
physical (or chemical) conditions cannot affect the decay rate. The stationary nature of the primary noise process 
is not violated in any other elements of the generator circuit in view of the appreciable level of the primary noise 
signal and the absence of amplifier elements in the circuit. 


4. The appreciable level of the primary noise signal also excludes the effect of parasitic regular components. 


5. The uniformity of the spectral density in the band with the upper frequency limit w}, is easily controlled 
by varying X. On the other hand, by measuring \ it is possible to determine the value of w}j for a specified value 
of n. 


6. The generator circuit provides the possibility of controlling the quantities A, Mg, O¢, and also provides a 
circuit for measuring the value of i. 


71. The simplicity of the generator circuit assures reliable operation without any need for setting or regulating 
during the operating process. 
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The paper studies proportional-integration and proportional -differentiation elements based on RC 
circuits with magnetic amplifiers; these elements are used in controllers. A method is given for 
choosing the parameters of the elements. Data based on an experimental investigation of highly 
stable proportional-integration and proportional-differentiation elements is cited. 


In order to design proportional-integration (PI) and proportional-differentiation (PD) controller elements it is 
customary to make use of RC circuits with electronic amplifiers. The reliability of PI and PD elements of con- 
troller units can be appreciably increased by using magnetic amplifiers. 


PI and PD elements are defined as elements with the transfer functions 
1 
Wp, (p) =k + 
PI T 
and i? 
W pp (p) = k (1 + Tap), (1b) 
where k is a proportionality coefficient, T; is the integration time constant, Tg is the differentiation time coristant. 





é (la) 


‘Proportional-Integration Elements Based on RC Circuits with Magnetic Amplifiers 
"Figure 1a shows the block diagram of a magnetic amplifier with an RC circuit. The control winding W; 

is supplied with dc current i,. A load resistance R; is connected in series with the operating windings W, (these are 

not shown in Fig. 1a) that are supplied with alternating current. In order to obtain integrating elements based on 

magnetic amplifiers a capacitance C (Rj; is the leakage resistance of the capacitor) and a resistor Ryg, are con- 

nected in the circuit of the negative feedback winding Wpfp. An increase in the amplifier gain is achieved by in- 

troducing positive feedback through the appropriate winding Wop: 


From an analysis of the block diagram shown in Fig. 1 we obtain an expression in operator form which relates 
the input voltage to the output voltage: 








OTR tee Fy + Hs bop* +- dip +- bs 
Mout = Yinks Rap tap bas * . 


™%  S T, Tp. e 
bale (Be): b= Te + ERE B) bm tt 
a= 7; (T +71) + Tet. (1+ 72) + 72TH (1+ 4), 
an=Te T+ Ti+ q [1 +7,(14+ 2) 4 p_1,(1 4 E): 
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a, = A. + 1 
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is the current gain referred to the control winding of the magnetic amplifier MA; 
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is the current gain referred to the negative feedback winding; 
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Fig. 1 
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is the current gain referred to the positive feedback winding; 
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is the time constant for the control winding circuit (L, and R, are the inductance and resistance of the control winding); 





sl = Woe F 
mnt) 


is the time constant of the negative feedback winding circuit (lz, and Re are the inductance and resistance of the 
negative feedback winding); 


T= We k 
(A + Ry) (1— WL L) 





is the time constant of the positive feedback winding circuit (Ig and R, are the inductance and resistance of the 
positive feedback winding); 
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TL = RC, TH =(R pe + A)C, T= 





are the time constants of the load circuit; T},=Rj,C is the time constant of the capacitor leakage circuit. 


Usually the time constants of PI and PD controller elements in units designed for the automatization of techno- 


logical processes are equal to several seconds; therefore T;, Ty, Ty « T and Tg, (Tj, >» T and Tp). In accordance with 
this formula (2) can be simplified considerably: 


a R,+R, pPTp +14+T7g, /T 
Lt “s P* fb fb / 7 ik 
out = Min hy Ro+R, PTH +T)+1+T/Ny * 








(2a) 


The dependence of the load current on the current in the control winding becomes 


pe - PT m +14+TH /T HY 
t= th 5}. +hFitT IT, 3) 





under these conditions. 
The transfer function for a magnetic amplifier with a real RC circuit (Ry, # ©, Rpg = 0) is written as 


PT g, +14+Tp /Ty, 4) 
1p(Tm +7T)+14T7/Ty,’ 





W (p) = k 


In order to lower the effect of capacitor leakage on the operation of units with RC circuits we attempt to 
choose Tg, and T « Ty. 


It should be noted that the magnitude of leakage resistance Rj), depends on the temperature of the ambient 
medium; therefore as the temperature increases T},, decreases. This can lead to instability in the operation of units 
with RC circuits. 


When the resistance of the negative feedback winding is appreciable or when there is a series resistance in the 


feedback circuit the magnetic amplifier with an RC circuit can be treated as a proportional-integration element 
(a PI element) since Tg, < T + Tp. 


We shall determine the parameters of this element on the basis of analyzing the transient responses for a unit 
step variation of the control current (Fig. 1b fort > 0, i, = 1). Here the load current is determined from the ex- 
pression 





Soe fi 4a TL orp (3) © 


‘L=hpe tt 1+TITy T +T 


We shall define the proportionality coefficient k of the element as the ratio between the current in the load 
of the element and the control current at the instant when a step variation of the control occurs (for t = 0): 


— Mitlimo _ T 6 
~ Bim TFT ™ 


The integration time constant T; can be characterized by a quantity which is the reciprocal of the rate of rise 
of the load current at the instant t = 0 when the current in the control winding varies stepwise: 


(Tm +TY 


HI. sss) ; 


T; = 
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TABLE 1 

















Pasamieter of PI elements 
| 1+T/T 
With an RC circuit when k ? A ee 
capacitor leakage is taken , hike kg (1 +7 / By) + hike 
into account 














The statistical error coefficient k,, shall be defined as the quantity which is the reciprocal of the value of load 
current in a steady -state mode for stepwise variation of the control current: 


ee ee ee (8) 
Tikes AG +TH IT * 





ks 


In contrast to the ideal PI element (1a), the PI element under study (4) has a static error that is inversely 
proportional to the coefficient k,. 


In order to obtain low static errors it is necessary to increase the gain ky of the magnetic amplifier with respect 
to the control circuit. However, this leads to an appreciable decrease in the integration time constant of the PI 
element. Therefore a magnetic amplifier with an RC circuit can be recommended as a PI element only for systems 
with low integration time constants. 


When there is no series resistance in the feedback circuit (Tg, « T) and Ry, = «, the transfer function for the 
magnetic amplifier is written as 


, i 
W (P) = ST 


A magnetic amplifier with an ideal RC circuit connected as shown in Fig. 1a is an aperiodic section; here 
k = ky is the transfer coefficient for the section and T is the time constant for the section. 


A PI element can be designed using two magnetic amplifiers (Fig. 2), The control signal in this circuit is 
applied simultaneously to the inputs of the magnetic amplifier with the RC circuit (MA,) and a summing magnetic 
amplifier (MA;), The output of the magnetic amplifier with the RC circuit is connected into the second control 
circuit of the summing magnetic amplifier. 


The transfer function for this element will be given by 


k kT k 
W (p) = hat ot = PE ee (9) 


If we take into account the leakages in the RC circuit of the amplifier MA,, then the transfer function for the 
PI element will be 





where k; and T represent the gain and time constant of the magnetic amplifier with the RC circuit, ky and k, rep- 
resent the gains of the summing magnetic amplifier. 
It is obvious that 
ksT <f 
kg +T/hy + &k, : 





and therefore the element under study (Fig. 2) can be used as a PI element in controllers. 
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Fig. 2, a) Circuit diagram of PI element in magnetic amplifiers, 
b) block diagram of PI element, 


at 
~ 


——— MMA AA A 
u Fete um 


iOr Wy My "i Me 





















































hy VK 
Wpfb1 Wpfb1 
fb Wpfb1 
- wi 
W, 
nfbi : nfb1 i 
Wilt Ubi +A, w’ 


Ww. w'. 
cl ? “ine cl 


Fig. 3. Circuit for the magnetic amplifier with an RC circuit. 














The parameters of PI elements connected as per the circuit in Fig. 2 are shown in Table 1, 


From Table 1 it is evident that the leakages of the RC circuit in the amplifier MA, of the PI element do not 
affect the parameters k and T;. In order to reduce the static error in a specified PI element it is necessary to in- 
crease the gains (ky, ks, ky) of the magnetic amplifiers with respect to their control circuits. This can be achieved 
without reducing the integration time constant if the coefficient kg is increased. The integration time constant can 
be increased in comparison to the time constant of an aperiodic section (for kyk, < 1). In the magnetic amplifier 
MA, it is also expedient to choose T « Ty. 




















Fig. 4. Circuit diagram of a PD element. 


Thus the latter circuit can be recommended as a PI element for systems with an appreciable integration time 
constant, 


The basic requirement which is imposed on magnetic amplifiers with RC circuits is the necessity of achieving 
an appreciable time constant T for a high null stability. 


The time constant depends on a load resistance R,, the capacitance of the capacitor C, and the winding ratios 
Wnfh/W 1 and Wppp/W,. The assurance of high null stability for the amplifier when the positive feedback coefficient 
is large presents well-known difficulties, 


We performed experimental investigations of the stability of differential magnetic amplifier circuits with 
simultaneous use of internal and external positive current feedback based on the current in each saturable reactor, 
and amplifiers with simultaneous utilization of internal positive current feedback with respect to the current in each 
saturable reactor and external positive voltage feedback. When the supply voltage varies by 210% the null stability 
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TABLE 2 


























Magnetic amplifier with Summing magnetic amplifier 
RC circuit 
winding data 
diameter, diameter, 
mm np mm 
Wii 550 0.2 Wh2 300 0.35 
Wi 550 0.2 Wi 300 0.35 
Wpfb1 275 0,15 Wnfb2 250 0.2 
Wfb1 0.15 Wnfb2 250 0.2 
Wnfb1 000 0.15 Wis 500 0.2 
Wpil 600 0.2 Wee 1000 0.2 
We 100 0.2 Wes 1000 0.2 
Wee 200 0.2 Wea 1000 0.2 
parameters 
ky 13 kg 220 
kg 650 kg 56 
Power gain 95000 Power gain 380000 
T (computed 100 sec 
value) 











of such circuits proved to be highly inadequate (withrespect to the power applied to the magnetic amplifier 
input the null drift was 10~ to 10~" watts). 


The maximum stability was achieved in a differential magnetic amplifier circuit which used external positive 
feedback with respect to the difference between the saturable reactor currents (Fig. 3). We shall indicate certain 
parameters of this amplifier: up; = 100 v, Ry; = 5.1 kohm, CLy= 1 wf, C = 30 uf, Rpfp.= 350 ohm (the internal 
resistance of the winding Wpf)1). This circuit is recommended for magnetic amplifiers with RC circuits. 


The null drift when the supply voltage varies by + 10% is equal to 0.8 - 10~° watts with respect to the power 
applied to the input; it is 0.6 v with respect to voltage (for t° = 20°C); when the temperature of the medium varies 
from 20 to 50°C and the value of u,, is nominal, the null drift is 0.5 - 10° w; the null drift for continuous operation 
over a period of 8 hours is 0.1 - 10°" w. 


The maximum power and voltage at the output of the magnetic amplifier are equal to 2.1 w and 100 v (the 
maximum output power, normalized to the input, is equal to 1.2 - 10 w). The null drift of the MA is insignificant 
when the supply voltage varies by + 10% with variation of temperature from 20 to 50°C or for continuous operation of 
the unit over a period of 8 hours. 


In a summing magnetic amplifier (MA; in Fig. 4) we use a differential circuit with internal positive feedback 
with respect to the current through each saturable reactor. The null drift, expressed in terms of the power normalized 
to the input, is 0,5 - 10~* w when the supply voltage varies by + 10%; it is 1 v at the output with respect to voltage 
(when the temperature of the ambient medium is 20°C), When the temperature of the ambient medium varies from 
20 to 50°C the null drift is 0.15 - 10* w (for a nominal supply voltage). The null stability for continuous operation 
over a period of 8 hours is 10-® w. The maximum load power is 0.66 w for a voltage of 20 v. 


The magnetic amplifier with an RC circuit and the summing amplifier have wound toroidal cores made of the 
alloy “79NMA" (ribbon thickness 0.1 mm; core dimensions 50 x 32 x 10), We used wire of the type "PELShO" for 
the windings. Rectifying silicon diodes of the type "D204" were used in the circuits. 


The data for the amplifier windings is given in Table 2 in accordance with the notation used in the figure. 
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Figure 5 shows the transient responses for the magnetic amplifier with the RC circuit shown in Fig. 3. The 
descending curves correspond to a stepwise variation of i,, from the corresponding value to zero, From Fig, 5 it 
follows that for the indicated circuit parameters a time constant of T » 100 sec is obtained (in the linear range of 
operation). The difference between the experimental and computed values of the time constant T is explained by 
deviations of the actual values of R;, and C from the nominal value, 
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Fig. 5 


The parameters of the PI element (computation is performed using Tables 1 and 2 for Ty, » T) are the 
following: the proportionality coefficient is k = 220; the integration time constant is Tj = 30 sec; and the static 
error coefficient is kst = 1.05 - 107, 


Proportional-Differentiation Elements Based on Magnetic Amplifiers 

The block diagram shown in Fig. 6a consists of two magnetic amp: MA, and MA; here MA, is con- 
nected into the feedback loop of MA,. The amplifier MA, represents the previously described aperiodic sections; 
MaAg is an adder. In order to derive the transfer function we shall study the block diagram of the element (Fig. 6b). 
Here the following notation is used: k, is the current gain of MA, with respect to the control circuit; T is the time 
constant of the aperiodic section; kg is the current gain of MA, with respect to the control circuit; k, is the gain of 
MA; with respect to the negative feedback circuit; kg = RL2/Rc, is the transfer coefficient of the divider (in order 
to eliminate the shunting effect it is desirable to make kg « 1). 


The transfer function of the element is written as 














ke hs (TP +1) 
W (p) = u Lar Tp+itkk eg (11) 
P+ 


We write the expression for the transfer function of the element in the following form when the effect of the 
capacitor leakage resistance in the amplifier MA, is taken into account: 


* By increasing kg and reducing k, and ky it is possible to increase the integration time of the PI element without 
altering kg. 
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(12) 


T T T 
The inequalities T> T+ kikakg and T+TINk > TFT Ty hahaha are always satisfied, and therefore 


the element can be used in controllers as a PD element. 





Ik? 





Fig. 6. a) Simplified circuit diagram; b) block diagram of the PD element. 


We shall determine the parameters of this element from a study of the transient response when the current 
icg in the control winding of MA, varies according to a linear law (Fig. 7). 


The current at the output of the PD element will be: when leakage of the capacitor in the RC circuit is 
neglected 


ius = ten phe [e+ (7 — spi) (1 —exe(— +p “* 9), 9 


and when the leakage of the capacitor in the RC circuit is taken into account 
red ka (i +-T / %) t+ ( T sia T ) x 
Ls TET Tt bike C+ TIN 1+ 71 Tit hikeka 


x (1 — exp (— ds LG) sth ‘))] . (13b) 








The proportionality coefficient k is equal to the ratio between the rates of change of the current in the load of 
the element and the control current in a steady-state mode (for t=): 


ri [# I.. | [2]... i, 


The differentiation time constant Tg determines the time by which the load current iL; leads the normalized 
current kicg in a steady-state mode, This time can be determined from the relationship (cf. Fig. 7) 


ins (t+ Ta) = igs (0) fort +00. (15) 
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TABLE 3 




















Parameter of the PD element 
Type of PD element k Tq kq 
ks kykgkg 
With an ideal RC 1+ hikka Tithhkka hykek 
circuit ’ 
With an RC circuit ke(i+ 7 / Tp r Kilkee kikakq 
taking capacitor PP 1 Tct krkekg) +P / Fy (l + 1 / Nyt kikekg) TTT A, 
leaka ; into 
acco 














The coefficient kq,which determines the quality of the differentiation, is equal to the ratio between the dif- 
ferentiation time constant and the transient response time constant [the time constant of the exponents in formulas (13)], 


The formulas for determining the parameters of PD elements with different sections in the negative feedback 


circuit are given in Table 3. 


From Table 3 it follows that all the parameters of magnetic amplifier PD elements depend on the transfer 
coefficient for the closed loop of MA, and MAj; this coefficient is equal to the product kykgkg. In order to achieve 
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Fig. 7. Transient response in the PD element for 
a linearly increasing control current. 


a better quality of differentiation one should choose 
kykgkg >» 1. Under these conditions the differentia - 
tion time constant Tg for a PD element with an 
ideal RC circuit will be a maximum and equal to 
the time constant T of an aperiodic section; the 
value of Tg for a PD element with an RC circuit 
will differ from T by a factor(1 + T/T},) in this 
case when the leakage of the capacitor is taken 
into account, The proportionality coefficient in 
the PD elements under study depends both on the 
transfer coefficient of the loop and on the gain 

kg with respect to the control circuit of the adder, 


The parameters and operational stability 
of magnetic amplifier PD elements will increase 
as the RC circuits which are used in them approach 
the ideal (i.e., as the leakage in the capacitors 
decreases), 
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The simplified circuit diagram of an experimental magnetic amplifier PD element based on the block diagra 
shown in Fig, 6 is given in Fig. 4. The parameters of MA, and MA; are cited in Table 2. The parameters of the PD 
element (the computation is performed in accordance with Table 3 for Ty), » T and kg = 0.05 ) are the following: 

the proportionality coefficient is k = 5.9 and less; the differentiation 
time constant is Tq = 97.5 sec and less; the coefficient for the quality 
[ sec : of the differentiation is ka = 36.4. 


- The designed PD element (cf. circuit in Fig. 4) was used jointly 
with the simulator unit "EMU-8." We obtained the dynamic charac- 
teristics of a PD element with its output supplying the integrator of 
the "EMU-8"; here the input control current of the PD element was 
varied stepwise (Fig. 8a). We investigated the characteristics of the 
PD element when the input control current is varied according to a 
79- NMA linear law at various rates (Fig. 8b). The differentiation time con- 
stants Tg in the linear range of operation of the PD element were 
equal to 80 sec in accordance with Fig. 8; i.e., they were approximate- 
ly equal to the time constant of the magnetic amplifier with an RC 
See circuit. In the nonlinear range of operation of the PD element the 
a} values of Tg increase. From the characteristics which are shown it 

is evident that the quality coefficient which characterizes the tran- 
Fig. 9 sient response in the PD element is close to 40. Thus the experi- 

mental parameters of the described PD element are in good agree- 

ment with the computed parameters (for Tj;,> 200 sec ). 


The characteristics under study were obtained for both directions of current in the control winding; their spread 
proved to be very insignificant, and therefore Fig. 8 shows the characteristics for only one direction of the control 
current. 


From Tables 1 and 3 it follows that the parameters of PI and PD elements are determined to a considerable 
extent by the quality of the magnetic amplifier with the RC circuit that is used in the networks. 


For specified core dimensions of the magnetic amplifier with the RC circuit, a specified current ratio in the 
load, and a specified wire diameter for the winding it is possible to compute the function T = f(Wpf/W 1). The 
corresponding curves for the core materials "79NMA” and “65NP" are shown in Fig. 9. The current ratio in the load 
of a magnetic amplifier with cores of "79NMA” is equal to 11; for cores of “65NP” it is equal to 6 (the dimensions 
of the cores are 50 x 32 x 10 mm and the diameter of the wire which is used for the operating winding is 0.2 mm; 
the diameter of the wire used for the feedback winding is 0.15 mm). 


It follows from the functions cited above that in a magnetic amplifier with an RC circuit and toroidal cores 
of “79NMA" it is possible to obtain a time constant equal to approximately 200 seconds for the relationships cited 
above (for C = 30 uf, when the ratio between the number of turns is Wnq,/W_ = 3). 


It is expedient to choose the value of W,4,/W , somewhat to the right of the maximum on the curves for 
T = f(Wpg/WL) (Fig. 9), since if this is not the case less power is drawn from the amplifier and supply sources with 
higher voltages are required. Therefore an insignificant error in choosing Wpfp/W can lead to appreciable reduc- 
tions in the time constant. An increase in the time constant of PI and PD elements can be achieved both by in- 
creasing the capacitance in the RC circuit and by increasing the dimensions of the magnetic amplifier cores. The 
expediency of the latter solution is determined in designing PI and PD elements by taking into account the engineer- 
ing requirements imposed on the control system. 
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Conditions are determined for obtaining a sinusoidal output voltage shape from frequency doubling 
magnetic modulators which have magnetic excitation and signal fields that are either parallel or 
mutually perpendicular. It is shown that these modulators have the highest amplification and sig - 
nal-to-noise ratio. A method is given for calculating such modulators. 


1. Higher Harmonics, Coefficient of Transmission,and Lower Threshold of Sensitivity 
Of all the known types of magnetic modulators those operating on the frequency doubling principle have poten- 

tially the lowest threshold ofsensitivity. The voltage appearing at a modulator output when a signal is present con- 

tains, as a rule, not only the second harmonic of the supply frequency but also higher even harmonics. For example, 

in a typical series output modulator which is supplied with sinusoidal current the amplitude ratios of the fourth, 

sixth, eighth, and tenth harmonics of the output voltage to the second harmonic are respectively 1.47, 1.01, 0.62 

and 0.43. 





A substantial content of the higher even harmonics in the modulator output voltage often makes it necessary 

to use special filters to suppress them.” The presence of these filters reduces the value of the signal that passes to 

the input of the subsequent electronic or semiconductor ac amplifier, increases the requirements on frequency sta- 
bility for the supply, narrows the passband of the circuit, and limits the possi- 

bility of introducing strong negative feedback. Of no less importance, the 

presence of a substantial amount of the higher harmonics at the modulator 








ee R output is linked with the reduction in amplitude of the essential double fre- 

‘ quency harmonic in the output voltage and a corresponding reduction of the 

w Ww, signal-to-noise ratio. One can easily be convinced of this if one considers 

| 900 the basic relations of a typical magnetic modulator designed, for example, 





there is no ac component in these windings under steady state conditions (i.e., 
ic = Ig = const ), We will assume also that the modulator is working into an 
—L 9 6 a open circuit, i.e., an infinitely large resistance is connected to the terminals 
ke of its output windings W2. In the analysis we will neglect the effect of 
Fig. 1 hysteresis. 


The current i in the exciting windings W, creates in the cores identical 
ac fields of intensity H = H(t). The control current I, creates a dc magnetic 
field in the first core: 








1 
OD) like the circuit of Fig. 1. 
We W, Wy We We will assume that there is connected in series with the control 
L a eo e— windings W,. of this modulator a large enough ac impedance Z, so that 
c 











AH =H, =-©°, (1) 


and in the second core a similar field of opposite sign. 





* When amplifying extremely weak signals it is often necessary to provide filters for the suppression of uneven har- 
monics due to imbalance, especially when using modulators with parallel fields, 








The value of the voltage at the modulator output is found from the formula 
d 
to = W,S < (B, — B,), @) 
where S is the cross-sectional area of one of the cores, and B, and B, are the values of the magnetic induction in the 


respective modulator cores. 
Confining ourselves to small signals it can be assumed that 


B, = B(H + AH) = B (H) + AH =B(H) + wAd, 
B, = B(H — AH) = B(H) —5, AH = B (H) — pad, 


where yt q = 4B/dH is the differential magnetic permeability. 
Substituting in (2) the expressions obtained for B, and B, we find 


d 
& = 2WaSAH —4, (3) 

The quantity #g has a maximum value (i gyax) when H = 0 and decreases monotonically with an increase in 
the absolute magnitude of H (Fig. 2a). If, for example, the excitation current i, and consequently the field in- 
tensity H, are varied according to a sinusoidal law (H= H,, sin wt, Fig. 2b), theng will vary at twice the frequency 
within the limits from ! dmin t© H dmax (Fig. 2c). 
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Fig. 2 
The differential magnetic permeability is an even function of time g(t) = 4 q(—t) and can therefore be rep- 
resented by means of the following series: 


bg= > Pan COs 2n at. 
n=0 


For the case under consideration, corresponding to Fig. 2c, it is not difficult to show that Han = 0. This is the 
result of the monotonic character of the change in 1 q during a quarter-period. 
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When wt = 0 


Ha (0) =Hdmax = po + pet pe tpet+.--, (4) 


and when wt = 1/2 


vg (+) =Hamin = Po — Ba toe — pe +... (6) 





Consequently 
Hdmax-Hdmin = 2 (ps +pe+po+. .). (6) 
For a sufficiently large excitation current (or field) Ugmax > Hdmin and 
Saad ae a Oe 1) 


We will now investigate the amplitude of the second harmonic of the output voltage, for which it is not dif- 
ficult to find from (3) and (7) 


E, = 40W,SAHp, = 40W,SAH [*asmar — (He + tio + Hie + ~)| (8) 


Obviously the maximum value of the second harmonic of the output voltage is equal to 


Eemax = 2@W Sud max AH, (9) 
which is achieved in the absence of the higher even harmonics of the output voltage. The relation is 
Ey 2Hs 


= Es (10) 
Emax Vdmax. ™+M+Pu+Put* 





The relationships obtained are confirmed experimentally as indicated by the results of measurements presented 
in Fig. 3. The modulator of Fig. 1 was constructed on spiral toroidal cores made of 80NKhS alloy. The amplitudes 


E2/Eamax 











of discrete harmonics in the output voltage (up to the fourteenth, inclusive) were measured with a voltage analyzer, 


which allowed the values of the corresponding components of the differential magnetic permeability to be calculated 
according to the formula 


seesaeaiideanadeae «1 
Hen= SaW SAH ~ ThaWWesl 
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The magnitude of Ez max is found by a measurement of the average value of modulator output voltage for 
which from (3) there can be obtained 


n/2 
2 
Bay =|= \ fo dat | — 
0 
Hdmin 
= |8fW, SAH \ dyg| = 8fWeSAH (uamax-Hdmin)- (12) 


HUdmax 


For #dmax > Hdmin We find from (9) and (12) 
_* 
Ez max * F Fay: 


The solid curve in Fig. 3 depicts the dependence of the ratio E,/Egmax = 242 /Hdmax on the amplitude H,, 
of a sinusoidal exciting field for the modulator as measured by the indicated method, At first, the ratio E2/Ezgmax 
grows with increasing H,,,, which is due to the diminishing value of Ugmin- When Hp, » 0.08 oe, Ugmax >! dmin: 
and further increase of Hm develops chiefly a greater content of the higher even harmonics in the curve Ug = Ug (w t), 
at the expense of a smaller second harmonic 2. 


The dashed curve in Fig. 3 is the function 





He a 
PatPe+Pithu 1(Hm), 


computed from measured values of 2, Hg, M49, and [1 44. 


As was to be expected, this function coming from formula (10) practically coincides with the dependence of 
E2/Egmax On Hm. The small discrepancy which is seen at the relatively large values of Hp, apparently is due 
still higher even harmonics which are beginning to develop at the indicated values of Hp. 


The transmission coefficient (amplification) of the modulator is reduced in the same degree as the second har- 
monic of the output voltage is reduced (with respect to the maximum possible value Ezmax): 


teats BE: ~ dae ,S “a 
Us The Rl : 





where R,, is the resistance of the control circuit. 


The lower limit of the threshold of sensitivity of the modulator is due to the magnetic noise [1], expressed 
by the formula 


e 


rm HGR 





watts, 


where E yy is the voltage of the magnetic noise in a given passband. By altering the shape of the exciting current 
curve of the modulator the ratio 2 / H#gmax can be changed. Assuming that the value of Ey4, moreover, is prac- 


tically unchanged, we find that P,, = cies? For example, if the excitation conditions of the modulator are 
2 


such that #12 is four times smaller than the maximum possible value of 0.5 Hdmax, then the lower threshold of the 
modulator sensitivity becomes 16 times higher than the value attained for a sinusoidal shape of the output voltage. 


It should be noted that such an output voltage shape is attained with a more uniform time rate of change in 
the magnetic induction of the cores than when the modulator is excited by means of a sinusoidal current or voltage 
(see below). Therefore the magnetizing increments (Barkhausen increments) also will be more uniformly distributed, 
which is apparently accompanied by a reduction in the absolute value of E,,. Considering also that a sinusoidal out- 
put voltage shape from the modulator often eliminates the need for an output filter together with the losses in it, 
in such cases a still more substantial increase in signal-to-noise ratio can be expected. 








2. Conditions for Obtaining a Sinusoidal Output Voltage Shape in Modulators with 
Parallel Fields 
Modulators withcores having magnetic exciting and signal fields in parallel (or antiparallel) directions are 


arranged in accord with the circuits presented in Fig. 1. For such modulators a sinusoidal shape of the output volt- 
age can obviously be obtained if 








Pd = Po + pe cos 2wt. (13) 


Considering that 4 gay = Ho + 2, it is possible to write 





Hd = Hamax(! —& + cos 2w2), (14) 
where 
a He val HPdmax™ Mdmin 
'= Hdmax 2 Hdmax SOM oe 


If the magnetizing curve of the modulator cores is given analytically in the form of a function B = f(H), then 
by solving the equation 


: (H) = Ug max! ral E + E cos 2a), (16) 
it is possible to find the law of the time rate of change of the exciting field intensity H or the magnetic induction B 


for which the modulator output voltage will be sinusoidal. 


Let us look at some examples, 


It is known that the average curve of magnetization (that is, the average between the rising and descending 
branches of the limiting hysteresis loop) for a number of materials used in modulators (for example, 80NKhS al- 
loy, ferrocart-2000, and others) can be approximated with suf- 
































Vv ficient accuracy by the arctangent [1] 
selon = = B, arc tgBH = /(H). (17) 
' 
ft 0 x = ie Toes From this expression it can be seen that when H ~ », the 
| induction B -» Bs, where Bs is the induction for saturation, More- 
over, 
8 
b . (4B 1m 
/ \ ; Hdmax = lim(T), = PB a8) 
aft P wi «x wt 
/ Substituting (17) in (16) and assuming for € the maximum 
é ud possible value, equal to 0.5, we find, taking into account (18), 
‘ee —_—_ = (1+ cos 2 @t) = cos*wt 
vA én # dmax “~~ 
i 0 sf x ar Jc wt or 
Fig. 4 BH = + tg wl. 


Substituting this value in (17) we obtain the required law for the change of magnetic induction: 
B=+— Bot = + 4/B,t. (19) 


Such a law for the change in magnetic induction is obtained when the modulator is supplied from a source of 
alternating voltage with rectangular shape and an amplitude of 


Vm = 8/W SB, (20) 
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provided that the voltage drop in the resistance of the exciting circuit can be neglected. In Fig. 4 there are shown 
diagrams of the change in supply voltage, in magnetic induction, and magnetic permeability for the case being 
considered, 


If the magnetization curve of the cores cannot be approximated with the arctangent, or if it is impossible to 
neglect completely the voltage drop in the resistance of the excitation circuit (for example, when B — B,), then 
the shape of the modulator output voltage will not be purely sinusoidal, However, the observed deviations in the 
magnetization curves for cores made of the magnetic materials most often used for magnetic modulators (for ex- 
ample, alloys 80NKhS, 79NM, 79NMA, ferrocart-2000), as compared with the magnetization curve corre- 
sponding to the arctangent (17), generally produce distortions in the shape of the modulator output voltage which 
can be easily compensated by a proper selection of the resistance R in the excitation circuit. By considering the 
curves in Fig. 5 one can easily verify this. 




















Fig. 5 
The modulator output voltage found from formula (3) can be also written in the following form: 


eo = 2W.SAH “hd 2B 








adB dt ° 
In the example considered above dB/dt = + 4fB,, and 
2BB, 
dud d n \ , B 
GB GB) igs & 8 ad a ad ai — 
TB / 


For toroidal cores made from the magnetic materials specified above the deviation of the relation dy g/dB = 
= f(8) from the sinusoidal (21) is of the type shown in Fig. 5a. In order to obtain a sinusoidal shape of the out- 
put voltage €, = (duiq/dB) - (dB/dt) (Fig. 5b), it is necessary that dB/dt should vary in accord with the dashed line 
in Fig. 5c. The variation dB/dt= f(wt) close to that required can be attained by connecting the proper resistance 
into the excitation circuit, In practice it is extremely simple to select experimentally the vaiues of a rectangularly 
shaped supply voltage and a resistance for the excitation circuit which provide a sinusoidally shaped output voltage. 


As an example there are presented in the table the results of measurements of the ratio E,/E; max.and a co- 
efficient of nonlinear distortion for the output voltage of a modulator 





VEG + B+ By + Biot Fie t Bia 
Kp = EB; 





made with cores of 79NMA alloy, for various excitation conditions. For each case here the same average value of 
the ac excitation field intensity was established, equal to 0.1 oe. 


Note that when the field intensity of excitation is increased, the shape of the output voltage deteriorates 
sharply for a modulator excited by sinusoidal current or sinusoidal voltage, but wtien excited by rectangular-shaped 
voltage, a curve of relatively low distortion can be obtained by a suitable change in the value of a series-connected 
resistance. 
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Method of Sinusoidal Sinusoidal Voltage of rec- 
excitation current voltage tangular shape 
E2/Eemax 0.56 0.68 0.93 
Kp 1.43 0,97 0.2 











In a modulator having cores of 79NM alloy which was supplied with rectangularly shaped voltage,still less 
distortion was obtained in the shape of the output voltage: Kp = 0.06 for E,/Egmax = 0.97. 


B é “Qs 


20 


IS 











wt 


It must be noted that in some cases the fixing of a certain shape ofthe excitation current curve which will 
give a sinusoidal shape for the output voltage can be done as expedient, Let us assume, as often occurs for cores 


made of transformer steel, that the average magnetization curve can be approximated with sufficient accuracy by 
the hyperbolic sine 


H =ashBB. 


In this case it is not difficult to get 


dB HU dmax Bdmax Hd max 


Pa GH ~coshBB ~ Vi-sh*pe 1+(2) ’ 
a 


i 
Hd max = a 


Substituting this expression for f'(H) in (16) and solving for H we find 








H 4E sin? wt (1 — E sin*® wt) 
Fa Se ise (22) 


In Fig. 6 there are presented graphs of the required shape for the excitation field (or current) intensity(in one 
half-period), which were plotted from formula (22) for various values of E. When — » 0.4, the required shape of the 
excitation current is attained with sufficient accuracy if the modulator is supplied from a source of rectangularly 
shaped voltage. For higher values of € it is necessary to take supplementary measures in order to get the required 
shape of the excitation current. A satisfactory method of changing the shape of the excitation current within broad 
limits is to connect a saturating choke coil in the excitation circuit of the modulator. 
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3. a Modulators with Mutually Perpendicular Fields 

design and the basic circuit of the modulator are shown in Fig. 7. A ferrite core 1 is composed of two 
identical halves having ring-shaped grooves in which the excitation winding W, is located. The control winding W, 
and the output winding W, are arranged uniformly over the entire length of the core. 


In this modulator the output voltage can no longer be calculated from formula (3) because the increase in the 
magnetic field AH = H, due to the input signal does not add algebraically to the excitation field H, 








Fig. 7 


i,e., the ratio of the magnetic induction B) in the longitudinal direction to the field intensity H | acting in this 
direction is equal to the ratio of the induction B, in the transverse direction to the field intensity H, acting in the 
transverse direction. Since Hjj = ASH=H_ isthe signal field intensity, then for the output voltage of the modulator 
we have 


éB d dp dB 
€o= WS, —- = W,S, AH F = WS, AH ag S . ew 


where Sj is the cross-section area of the core for the longitudinal magnetic field. 


The principal difference between this formula and formula (3) for the modulator with parallel fields is that 
in expression (24) for the modulator with mutually perpendicular fields py is the static magnetic permeability of 
the core while 1 q in (3) is the differential magnetic permeability. Moreover, both 1g and pt are determined by the 
average magnetization curve,which does not have points of discontinuity. Therefore 


Amax = (Zone ip (FA ps = Bdmax- 


Confining ourselves to small signals it can be assumed that the function 1(t) is identically determined by the 
law of change of the transverse induction component B, (or H, ) and does not depend on the value of the signal. 
Therefore instead of equation (16) it is possible to write the following condition for the sinusoidal voltage at the 
output of a modulator having mutually perpendicular fields: 


B 
H == Hmax(1 —§ + § c0s 20%), (25) 


If the magnetization curve is approximated by the arctangent (17), then a sinusoidal output voltage will be ob- 
tained for different excitation conditions than for the case of the modulator with parallel fields. 

Actually, for the case under consideration it is possible in practice to get a sinusoidal output voltage shape by 
varying the induction according to a sinusoidal law 


B, = B,sinat. (26) 





. 


e 





From (17) and (26) we find for the change of magnetic permeability with time 


nm 
—~ sin wt 
emg a! OpGley aes ae (27) 


— — 7 


Bmax Umax Ay tg + sin wt 











The solid curve in Fig. 8 represents the required law for the change of 4 when € = 0,5, as established by the 
right-hand part of equation (25), The dashed curve was computed from formula (27) and corresponds to the change 
of induction in accord with (26). 
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In Fig. 9 curve 1 represents the function B = f(H) which is described by the arctangent (17). It is interesting to 
note that if the magnetization curve can be approximated by the function 


B 
oukinana @8) 
ies ter yuh mc 
B {+ cos 5— 


which corresponds to curve 2 in Fig. 9, then the modulator output voltage will be sinusoidal when the change of in- 
duction is directly proportional to time, in accord with (19), i.e., when the modulator excitation is from a source of 
ac voltage of rectangular shape. Indeed, substituting (19) in (28) we find 


BBB, 
ATT ss 





+ cos 2wt) = 0.54,,4,(1 + cos 21). 


From Fig. 9 it is seen that magnetization curve 2 found from formula (28) coincides in the area of low induc- 
tion with magnetization curve 1 which corresponds to the arctangent (17), but in the area of high induction it lies 
below the other. This type of magnetization curve is observed for cores of the kind represented in Fig. 7 as a result 
of non-uniform core magnetization by the transverse magnetic excitation field. The variation of magnetic path 


1255 








length over wide limits and also the presence of a varying area of core cross section for the transverse field lead 
to incomplete saturation, i.e., at a time when part of the core is saturated, another part proves to be unsaturated. 
Therefore the curve representing the ratio of the average value of magnetic induction over the cross section of the 
core will lie below the magnetization curve of the core material in the area of high induction. 


An experimental study of a large number of modulators having mutually perpendicular fields which were made 
with cores of ferrocart-2000 demonstrated that it is possible to get in practice a sinusoidal output voltage shape 
by an appropriate choice of the amplitude of a rectangular-shaped supply voltage and of the resistance in the excita- 
tion circuit [2]. In Fig. 10a, there is presented a typical oscillogram of the supply and output voltages obtained for 
such modulators. — 


4. Modulators with a Two-Phase Supply 

The most convenient source of rectangular-shaped voltage for a magnetic modulator supply is a magnetic- 
transistor ac to dc converter which is characterized by high efficiency and small size (see, for example, [1], ch. 7 
to 15). 


The use of such converters allows the frequency of the modulator to be increased, and thus its coefficient of 
amplification and signal-to-noise ratio are increased, Also among the advantages of magnetic-semiconductor con- 
verters are the facts that they serve as unusual stabilizers for the modulator supplies and they can be built with a 
negligible second harmonic component in the output voltage. Finally, it is not difficult to construct, on the prin- 
ciple of the magnetic-transistor converters, a two-phase supply source having two rectangularly shaped output volt- 
ages shifted in phase by 90° [2]. 

The use of two modulators supplied by a two-phase source permits the second harmonic to be eliminated in the 
control circuit formed by the series connection of the modulator control windings. Therefore the need of inserting 
a filter or an impedance which is rather large for the double frequency current in the control circuit is eliminated, 
The output windings of the modulators are connected oppositely with respect to the control windings. In addition, at 
the modulator output there will be no voltages with frequencies of 4nf, where n = 1, 2, 3, ..., which provides a 
further improvement of the output voltage shape [2], An oscillogram of output voltages shown in Fig. 10b was taken 
after combining two single-phase modulators having mutually perpendicular fields in a two-phase circuit. 


5. Design of Modulators Having a Sinusoidally Shaped Output Voltage 
~ For the approximate design of modulators having an output voltage shape which is nearly sinusoidal, it is suf- 
ficient to know only the value of the maximum magnetic permeability on the dynamic hysteresis loop: 


Bg = ia = 
max ‘Bae dH - 











This quantity determines the value of the maximum permeability on the average magnetization curve 


. B 
lim = . 
Bmax— H ~ “d max 


Substituting (17) in (3) we find for the amplitude of the second harmonic at the output of a modulator having 
parallel fields, when it is operating in the open circuit condition (infinitely large load resistance) 


20§2W We Sd max 





Ey, = 40€W.StdmaAd = i Icy (29) 
where 
Wel 
AH =——- 


For a modulator having mutually perpendicular fields from (24) and (25) we have 


14max = (30) 


Ey) = 20 








de 


t 





We notice that 





2W Wc Sitd max 
I = Mmax 


and 





W.WS max 
i = Hmax" 


where M,,,x is the maximum value of the coefficient of mutual inductance between modulator windings W,, and W,. 
Therefore for both the the modulator with parallel fields and the modulator with mutually perpendicular fields 


E, = 206 Max! c. (31) 
In the particular case when W, = Wg, or when one of the windings is used for an input or an output we get 
E, = 20§Lemax! ec (32) 


where Lo max is the maximum value of inductance of the input (output) winding. 


We recall that the maximum possible value of € is equal to 0.5, Moreover, the value of € in the case of 
parallel fields is always larger than the value of € for mutually perpendicular fields, In the first case it is possible 
to take € = 0.5 without a substantial error, and in the second case, as experimental studies and calculations showed, 
E ws 0.4 to 0.45. 


The constant component of induction By which appears in the core due to the signal is calculated from (16) 
and (25) to be equal to 


By = poAH =tmax (1 — =) AH. (33) 


This expression can be used for the calculation of transient processes in the control network, which are de- 
scribed by the following linear differential equation® : 


dB di 
2WeS = +La = +Heic =Ve 


(le + ta) $+ Relic = Ve" 


where Lg is the inductance of the choke connected in series with the control winding and L¢ is the inductance of 
this winding for input signals, equal in the case of parallel fields to 


2w*.S (4 —8) Ww 
Le = SO dma. = 51 —)Minax 





and in the case of mutually perpendicular fields to 





wis, (1—8) W 
—s at pe nman = FE (1 —8)M max” 


When a resistance Ry, is present in the output of the modulator, a current i; flows through its winding Wz creat- 
ing a field intensity Hy which has an effect on the value of the emf induced in the output winding of the modulator. 
The value of the load current is found from the differential equation of the load circuit 


* This equation is written for a modulator with parallel fields and two cores. For a modulator with perpendicular 
fields it is necessary to substitute 2S = S}j- 
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dughH 
2W,S 2 — +O, = 0. (34) 
For simplicity we will only consider the effect of the second harmonic of the current i;. Then 


i; = /, sin (2@t — g). 


This current creates a field intensity 


w 
ee pad = 2 sin (20t — @) = Hy sin (20t — 9), 


which adds to the signal field H, = 1.W_/! , Therefore we have for the value of AH in formula (34) 
AH =H, + H, sin (2ot — 9). 
Assuming as before that the time rate of change of 11 q is determined by the supply voltage and is described 


by formula (14), and considering only the second harmonic, we find from (34),after substitutions for the values of 
Hq, 4H, and i,, 








2oL| — 72; He sin Zot + Hy cos (ot — q)| +R, Hy sin (20t — 9) = 0, rm, 
where 
2w2s (t—E)p Ww? Ww 
17 #____ + dmax_ ate = We (i — &) Maa (36) 


is the average value of the inductance of the modulator output winding. 
Assuming in (35) that wt = 0 we find 


Substituting in (35) the value 2wt = 1/2, we obtain 


2o0L, 7: 
H, — 
VR + (2@L,)* 


H 
- (37) 








Now it is possible to find the amplitude of the second harmonic voltage at the output of a modulator with an 
impedance load: 


3 
2R, oL, -——= 
Hl °4.7—-F WW, 
ve = 1 =—-—R = me 
i =, Ve +oLy “* 








(38) 


Expressing Lj in the numerator of (38) in terms of Mmax, according to (36), and considering (31), we find 


mp 6) Scns a (39) 
VR + Col, 
This formula is the theoretical basis of the empirical rule established earlier for the approximate output voltage 


calculation of a double frequency modulator, in which the modulator is considered as a generator of an emf E, which 
has an internal inductive reactance (see[1], ch.9, sec. 5C, and [2}): 











Ww. 
Zi = 20l, = 2obery . 
c 


It is obvious that by this method it is also possible to calculate approximately the effect of a load having an 
inductive or capacitative character. 


With an operating load the maximum value of output power, which is achieved when Ry = 2wLj, is found from 
(39), taking into consideration (36): 





2 
nom wh 05( Egle ek, o 


For the power amplification coefficient we have 


K, = - = ~ ~ 0:5(72 5). (41) 


where R,, is the resistance of the control circuit. If the modulator control circuit includes, besides an inductance L, 
in the winding Wc, only a resistance Rg then the ratio L./R¢ = T¢ is the time constant of the control circuit and 


eoa(cight 2 


For a two-phase magnetic modulator the values of E,, Lj, and Re are doubled. Therefore P; is also doubled 
with respect to the value given by formula (40) while Kp, rc, amd Ky/t¢ keep the same values as for a single- 
phase modulator. 


SUMMARY 
The studies which have been made show a substantial possibility of achieving by simple means in practice a 
sinusoidally shaped voltage at the output of modulators operating on the principle of frequency doubling. To this 
end it is satisfactory in the majority of cases to supply a modulator from a source of rectangularly shaped voltage. 
In addition to this the amplification coefficient of the modulator increases and the signal-to-noise ratio is improved. 


For an approximate calculation of the static and dynamic characteristics of modulators having a sinusoidally 
shaped output voltage it is sufficient to know only the maximum value of magnetic permeability on the dynamic 
hysteresis loop of the core material. 


The conditions for producing a sinusoidally shaped output voltage for modulators with parallel and with mutual- 
ly perpendicular fields differ in principle. However, for the same cores and winding data these regulators have prac- 
tically the same static and dynamic characteristics. 


When two of the modulators having mutually perpendicular fields are supplied from a two-phase source of 
rectangularly shaped voltage, it is possible to set up a magnetic modulator which has an extremely low threshold of 
sensitivity without filters in the supply, signal, or output circuits [2]. 
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Expressions are derived for the transmission rate of telemetered data, taking into account both 
equipment errors and the effect of interference in the communication channel. 
The concepts of resultant precision and the criterion of interference stability of tele- 


metering are presented. 


In the first portion of article [1], in the analysis of the efficiency of information transmission without taking 
account of interference in the communication channel, it was considered that the telemetering precision Nag is deter- 
mined only by the equipment error of the conversion transducer _ used (which was also considered identical to the 
discreteness error), that is, it was considered that 


Na= z ; (1) 


where 5, is the equipment error and Na is the equipment precision of the telemetering. 


In the general case, taking into account the effect of interference arising in the communications channel, 
we can say that the information received reproduces the value of the telemetered parameter with a resultant error 6, 
determined both by the equipment error 54 and by the error 5; caused by interference, which are, in general, quan- 
tities independent of each other. Before defining 5 we call attention to the fact that the equipment error 5, is 
usually taken as a maximum. It represents the maximum of the absolute deviation (averaged over the set of values 
at a given point of the scale), referred to the range of measurement. The error 5; caused by interference is defined 
as @ mean-square error or a mean error (see [2-4)). 


We define the resultant error § as a maximum, similar in physical meaning to the equipment error. 


If 54 is a mean-square error, the square of which is proportional to the dispersion of the distribution of errors 
caused by interference, then by the rule of summation of the dispersions of independent quantities we obtain 


(2) = (a 
= EPA : 


where r and rg are the proportionality factors between the maximum value and the mean-square value of the resultant 
error and the equipment error respectively.” 





* Obviously, if & (or 5,4) obeys the normal distribution law, then the corresponding r (or ra) may be set equal to 3. 








If 5; is a mean error, then, adding the mean values of the errors, we obtain the following expression for the 
resultant error: 


$= F844 r'd; » 


TA 


where r' and r', are the proportionality factors between the maximum and mean values of the resultant and equip- 
ment errors respectively.° 


We now define N as the resultant precision or resultant number of quantizing levels into which we subdivide 
the scale of measurements in 1eception as a result of the effect of interference and instability of equipment: 


N= 1/6, (3) 


We shall now proceed to define the rate of information transmission in telemetering. 


For a uniform distribution of the measured parameter (information) and for the new interpretation of the tele- 
metering precision N, the information transmission rate R in the presence of interference in the communication 
channel is expressible by the same well-known general formula 


R — 's® 


— (4) 


which was used in the case with no interference [1]. The difference here consists in the fact that instead of the equip- 
ment error the formula contains the resultant precision N. In this formula T is the transmission time of one measure - 
ment. 


Let us consider the value of log N, According to (3) 
log N = — log 6. 
Subdividing this into two terms, of which the first does not depend on interference, we obtain 
log N = — log}, — log ~ : 
A 
We introduce the notation 
5, 
v= > (5) 


We shall call Q the criterion of interference -stability of telemetering transmission. The fact is that the error 
5; caused by interference is not yet an exhaustive criterion of interference-stability. Such a criterion must be based 
on a comparison of the equipment error 6 with the resultant error 5 or with the error caused by interference 5;.° * 
In particular, the ratio 
54 
== = (5a) 
BT 3; 
may also serve as a criterion of interference-stability, However, the criterion Q(1 > Q> 0 when 0< &:< @) is 
preferable in the present case, since it leads to simpler expressions for the transmission rate. 


Taking (1) into account, we obtain 


log N = log Na — log = log NaQ. 


* The factor r’ (or r'a) is defined by the nature of the functional dependence of the absolute deviations on the value 
of the measured parameter. For example, for a linear dependence we have a factor r° (or r',) equal to 2. 

** We note that an analysis and comparison of interference-stability of telemetering methods should also be carried 
out using as a basis a criterion of the type of the proposed interference-stability criterion Q, whose value depends on 
the relative values of the errors 54 and 6 (or 5, and 5;). 








We make use of the expression (4), The rate of transmission of information in telemetering in the presence of 
interference in the communication channel is equal to 


1 
log NV, log = 
R= > (6) 





In accordance with the structure of this formula — two terms on the right side — we can represent the trans- 
mission rate as consisting of two parts: 


where Rg is the transmission rate in the case when there is no interference, a rate determined by the equipment error 
(it was just this value of the transmission rate that was considered in [1]), AR, is the loss in transmission rate caused 
by interference in the communication channel. 


From (2), (5), and (6) it is clear that if the error 5; caused by interference becomes zero, then there will be 
no loss in rate (AR; = C), since the interference-stability Q becomes equal to unity. 


Let us consider another definition of the transmission rate of telemetered information, suitable in a number of 
cases for the analysis of methods of telemetering transmission (see [1]). The difference in this definition consists in 
the fact that the rate of transmission of telemetered data (which in the present case we designate by V) is expressed 
not by logarithmic units per second but by the number of quantizing levels per second 


V =N/T. (8) 
The resultant precision N may be represented as consisting of two parts: 
N = Na—AN;,, (9) 


where AN; is the loss in precision caused by interference. 
Making use of (1) and (3), it is easy to convince ourselves that 


AN = Na(i-—@),N =NaQ. 
On the basis of (9) we represent the rate V also as consisting of two parts: 


V =V,—AVj. (10) 


Here Va = Na/T; the second term AV; = AN/T is the loss in transmission rate caused by interference. 


In [1] we obtained expressions for the transmission time of one measurement T for the basic methods of trans- 
mission of information in telemetering. In [2-4] formulas were derived for the errors caused by interference in trans- 
mission under conditions of weak and strong interference. 


These data are sufficient to obtain on the basis of (7) and (10) expressions for the transmission rate R or V, as 
well as the transmission efficiency criterion R/W considered in [1] (where W is the bandwidth occupied in the com- 
munication channel) for each method of modulation. 


The values of R and R/W (or V and V/W) for concrete telemetering systems make it possible to estimate and 
compare the transmission efficiencies of different systems. However, since the expressions R and R/W (or V and V/W), 
obtained in this manner are functions of the basic parameters characteristic of the telemetering, it is of interest to 
carry out an analysis of these expressions for the purpose of finding the conditions of most efficient transmission. 








e of 


6) 


or of 
's in 


8) 


(9) 


(10) 


‘ans~- 
trans~ 


, as 
om- 


and 
| V/W), 


: 
7 








1, 


LITERATURE CITED 
N. V. Pozin,” The efficiency of information transmission in telemetering: 1. Analysis without considering in- 
terference in the communication channel,” Avtomatika i Telemekhanika, 22, No, 9 (1961). 
V. A. Kashirin and G, A. Shastova, "The interference~-stability of the transmission of telemetered signals in a 
channel with fluctuation interference," Avtomatika i Telemekhanika, 19, No. 8 (1958), 
N. V. Pozin, "On the interference -stability of pulse-frequency telemetering,” Avtomatika i Telemekhanika, 
19, No. 10 (1958), 
N. V. Pozin, "The interference-stability of pulse-width and pulse-time telemetering for strong fluctuation 
interference,” Avtomatika i Telemekhanika, 20, No, 2 (1959). 


All abbreviations of periodicals in the above bibliography are letter-by-letter transliter- 
ations of the abbreviations es given in the original Russian journal. Some or all of this peri- 
odical literature may well be available in English translation. A complete list of the cover-to- 
cover English translations appears at the back of this issue. 








POTENTIAL AND REAL NOISE STABILITY OF MULTI- 
CHANNEL RADIOTELEMETERING SYSTEMS WITH TIME 
DIVISION OF CHANNELS UNDER WEAK FLUCTUATING 
NOISE 


A. F. Fomin 


(Moscow) 

Translated from Avtomatika i Telemekhanika, Vol. 22, No. 10, 
pp. 1405-1415, October, 1961 

Original article submitted February 8, 1961 


The potential noise stability is determined for various methods of transmitting multi-channel 
telemetering information with time division of the channels. A comparison is made of the 
various transmission methods for real and potential noise stability. 


Introduction 

Although at present there is a sufficiently large amount of work both on potential [1, 3-5], and on real [2, 6, 7] 
noise stability, yet not all radiotelemetering systems (RTS) with time division of the channels (TDC) have been 
examined with respect to potential noise stability, and there has not been presented a comparison of potential and 
real noise stability for many systems. 


In the present work, there will be used as a criterion for evaluating potential noise stability with weak 
fluctuating noise the value of the mean square error measured at the receiver output, taken for the entire range of 
parameter values, which is found from the formula [1, 4] 


Cis 3 
o5 ee a? (1) 


8 ) (2 
—T/2 


As a criterion for evaluating the noise stability of real receivers there will be used the value of the mean square 
error of the measurement, also taken for the entire range of parameter values, which we will determine from the 


formula 
Vp \? 
P 
ees - 2 
(er) (2) 


In the present work the ratios (Vs/Vp)out are utilized for various systems, as given by Nichols and Rauch [2] 
and other authors [6], 


When determining the real noise stability of radiotelemetering systems with time division of the channels it 
is assumed that the commutator of the receiving apparatus is rigidly synchronized with the commutator of the trans- 
mitting apparatus, and that in the output of each channel a low frequency filter (LFF) with a frequency cutoff AF = 
= Fmax 18 applied. 


We use the following designations: A(X, t) is the signal which is a function of time and the parameter being 
measured; F; is the interrogation cycling frequency of the parameters being measured; AF is the bandwidth of the 
receiver channel; F, is the bandwidth of the video channel; F,,,x is the maximum frequency of the parameter being 
measured; fpis the maximum deviation of the carrier frequency; M is the coefficient of amplitude modulation of 











.7 
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the carrier; N is the number of channels; T is the interrogation cycling period of the parameters being 
measured 


Vg is the effective voltage of the unmodulated carrier; V,, is the maximum amplitude of the carrier voltage; A is 
the parameter being measured in the interval + 1 (the value in this interval is equally probable); ar is the utiliza- 
tion factor for the channel time (aT = T/Tc); 5 is the mean-square error; o is the unit intensity of the fluc- 
tuating noise (v/ V cps); Wy is the angular frequency of the carrier; wp is the maximum angular deviation of the car- 
rier; T» is length of a pulse; Ar is the maximum deviation of a pulse; rf is the length of the rising (falling) front 
of a pulse. 


1. Pulse Amplitude Modulation-Amplitude Modulation System (PAM-AM) 
Potential noise stability. For a system with PAM-AM modulation the signal can be written in the form 








{Ym [1 -}- MAJ sin (wot ++ @) for = <t < 


A (h, t) = 
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0 fe <tc 5 
For this signal we obtain 


A’ (A, t) = Vm M sin (wot + @) for —BcticS, 


whence 


™%/2 t%/2 
. . OA (A, t) 2 ” ; ; { } 
[= \ e ai | dt =Vin M? \ sin® (@,t -|- @) dt = > Vin M*%. 


—t,/2 —t,/2 





When calculating the integral it is assumed [1] that for 





tT >= sin*® (Wot + @) = 5 ‘ 


We will use this assumption in the future without additional reservations. 
Substituting I in formula (1) and taking M = 1 we obtain 


wa a 
PR ie =. (3) 
P 8/ Vinte 


For a comparison with other systems it is convenient to convert from the power (maximum amplitude) of the 
signal in a pulse to the average power (average effective voltage) of the transmitter. When M = 1 for an N-channel 
system we will have approximately 


T 
vio= 2V5 ow: 
Substituting this value in formula (3) we obtain 


oP anne (4) 
P aver 


Real noise stability. On the basis of an examination of the signal and noise transmission processes in receiving 
apparatus it is possible to obtain the following expression [2]: 








a ES ee (5) 
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The mean-square error of the measurements will, in addition, be found from the expression 
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The ratio of the measurement errors during reception in real and ideal receivers will be 


6? 
2 r AF (7) 
Pr, ="F = oo « 
- $5 max 
Thus a real receiver with the PAM-AM system provides the potential noise stability only in the case where 
demodulation of a pulse series is accomplished with the help of an ideal LFF with a frequency cutoff equal to the 
maximum information frequency (AF = F,)4x)- 


2. Pulse Amplitude Modulation-Frequency Modulation System (PAM-FM) 
Potential noise stability. For a system of the PAM-FM type the signal can be represented in the form 








eee ee eel for —% /2<t<%/ 2, Ps 


Vin ©9S (Wot + 9) for —%T/2< t< — %/2. 
The expression for the error in a multi-channel system (r= T/ Nat» Vy, = 12Vg) will be 


Ran — i 
P 2 wr? 


Real noise stability. The expression for the mean square error in a general form, obtained on the basis of 
Nichols and Rauch’s data [2], is 


N3 ot (9) 
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The ratio of the errors determined from the receiving condition in a real and an ideal (per Kotel’nikov) receiver 
is found from the expression 


6 (10) 
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(11) 


In the special case when F, = 3/ r» and AF = Frmax we get P, = 1. 


It should be noted that a calculation of video bandwidth according to the formula F, = 3/T, obviously gives a 
somewhat higher value for F,. If we take F,< 3/7», then we find that the real noise stability according to Nichols 
and Rauch is higher than the potential according to Kotel'nikov, which is very unlikely. However, a detailed dis- 
cussion of this question is beyond the scope of the present article. 


3. Pulse Amplitude Modulation-Phase Modulation System (PAM-PM) 
——Feealal netoe wabiiiy- The signal can be represented in the form 


areal Mae for —%H/2<t< %/2, (12) 
Vin ©08 (Wot + Go) for %/2<t< —-x,/2, 








where @p = fp /F, is the maximum phase deviation of the carrier frequency. 
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On the basis of formula (1), taking into account that r9= T/Ndy, = fp/Fs, and Vp, = /2Vg , we get 


Now" F (13) 


Real noise stability. The mean square measurement error at the output of a real receiver is found from the 
expression [2] 








No F AF 
=, (14) 
4VEfn 
The ratio of the errors found from formulas (13) and (14) will be 
br 
p= 7 = 2AFT. (15) 
Pp 
In the special case when AF = Fy,,y 
p.' = 1. 


From the comparison given it is seen that the noise stability of a real receiver for the PAM-PM system with 
the proper choice of video bandwidth in the receiver and with an ideal LFF having a cutoff AF = F,,,, in the PAM 
demodulator likewise achieves the potential noise stability. 


4. Pulse Position Modulation-Amplitude Modulation System (PPM-AM) 
Potential noise stability. The signal for such a modulation can be represented by the expression 








A(h, t) =VmF [é —(t.+ “ r)) cos (wet + 9). (16) 


We will find the mean square error by superposing on the signal a low intensity noise during reception by an 
ideal receiver. To determine the derivative A‘ (A, t) we will substitute the variables #* = ¢t — (t, + a h): 


; OF dt* OF A 
A’ (dy t) = 55 5- = Vinge pr 008 (wot +), 


T/2 
I =TA*(@, )= {Vi A? [a] dt. 
—T /2 


The mean square measurement errors are 
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It is evident from these formulas that the error will be less,the greater Ar and the greater the unit energy 
of the fluctuation F(t )/dt*, i.e., the error depends on the pulse shape, or more exactly, on the type of rise in 
the leading and trailing edges of the pulse. 


If it is assumed that the rise and fall of the pulse occur instantly, then OF*/ dt «and 65 = 0. 


We will consider the case of a triangular pulse of 27 duration and with a linear rise and fall of the wave- 
fronts. Then Tr 
oF i (% y wii \ 
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Substituting the expression obtained into formula (17) and converting from the signal power in the pulse to the 
average power of the transmitter, for an N-channel system we find 


is T so 


TPN’ 
As a result we get 
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(18) 
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If it is assumed that the envelope of the signal pulse is represented by a function of the form 


ie (19) 


* 
then the expression for the mean square error can be obtained in the form [1] 


(20) 


A comparison of formulas (19) and (20) shows that when the parameters o, N, a7, Vg, T, and Fs=1/ rp 
are identical the system using pulses with the envelope F(t) = sinQ,t /Q,t has a somewhat higher noise stability. 


Real noise stability. We will determine the mean square measurement errors at the output of a real receiver 


operating on the leading edge of the pulse, which increases according toa linear law. The receiver operates when 
the pulse amplitude attains a certain value. 





The mean square displacement value of a triangular pulse due to noise is found from the formula [5] 





ree o 
At, = . (21) 


(when A is changed by unity the pulse is shifted by A 1 / 2), 
Converting to the average transmitter power for a multi-channel system, we obtain 


o?N3q2 


Qa (22) 
"4 Err 


It should be noted that such mean square measurement errors will appear while receiving one-sided Pulse 
Width Modulation (PWM) signals on a real receiver. 


The ratio of the measurement errors during the real and ideal reception of PPM-AM signals on the assumption 
of a linear rise in the wavefront of the pulse will be equal to 


p2 2 
ey es (23) 


If F, = 1/rp, then P2 = 2, 





18) 


19) 


20) 


ver 


21) 


(22) 


ption 





The mean square measurement error at the output of a real receiver while employing limiting short pulses 
(19) according to the data of Nicholsand Rauch is found from the expression 


o* Na? F? 
8a =a — (24) 
FS , 
From the ratio of Eq. (24) to (20) we get 
R=? yr ott: 


max 


From the comparison given it is evident that the real noise stability of this reception method for various pulse’ 
shapes will be 72 times lower than the potential noise stability of the transmission method, 


These facts agree completely with the facts presented by V, A. Kotel ‘nikov in [1]. In this work it is shown 
that the potential noise stability of PPM-AM signals can be completely realized only when the detection of the time 
position of a pulse is at the maximum pulse amplitude. 


5. Pulse Position Modulation-Frequency Modulation (PPM-FM) 

In the work of Nichols and Rauch [2) it is indicated that the use of FM and PM with PPM gives no advantages 
over PPM-AM because the basic merit of PPM is that the transmitting apparatus is only switched on for a small time 
interval. We will examine this situation with respect to the potential noise stability. 


For systems of the PPM-FM type which utilize pulses with envelopes described by expression (19) the signal 
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For this signal we find 


A’ (A,.t) = 
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For the calculation of the error we determine the integral 
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Substituting the value of I in (1) and taking into account that V,, = /2V_ and Ar= T/ Nay, we obtain 
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We will make a comparison of the PPM-AM and PPM-FM systems for the condition of the same average 
transmitter power Vee the same speed of operation T, number of channels N and a-y, the same PPM signal para - 
meters, and unit noise intensity o in the communication channel. 


To this end we will take the ratio of Eqs. (20) and (26): 





6p (PPM~AM) _ sar(a2) (1) 
5 1{PPM- FM) Qs / Ta, 


From (27) it is seen that the PPM-FM system in comparison with the PPM-AM system provides a much higher 
potential noise stability only in case the frequency modulation index is larger, i.e., 


WY TE, 
\ >a (28) 





It is easily shown that the fulfillment of this inequality requires a high index of frequency modulation which is 
practically unrealizable when using short pulses. 


Thus the PPM-FM system is actually worse with respect to noise stability than the PPM-AM system for the con- 
tinuous transmission of signals. 


6. Pulse Width Modulation-Amplitude Modulation (PWM-AM) 

In radiotelemetry systems one-sided PWM is the most widely used. Therefore we will compare systems with 
respect to the real and potential noise stability on the assumption that a one-sided PWM is used in which the leading 
edge of the pulse is rigidly clamped during reception. 


Potential noise stability. For such a system the signal can be represented in the form 
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for te << # < te + SE (4 +a), 
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for te > t > te +55 (1 + A). 





For simplicity we will assume that the trailing edge of the pulse falls off according to a linear law and the 
duration of the edge is rp, In this case 
T 
8 ng 2 vm 4 
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The mean square measurement error will be determined by the expression 


52 er 


P Vm 4 


(29) 


It can be shown that for two-sided PWM the error will be 72 times greater. 
Converting to the average effective signal voltage of the transmitter 
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(where 41/2 is the average pulse length) and substituting Ar = T/Na7T, we obtain 
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Real noise stability. The signal-to-noise ratio at the output of the ith channel of a multi-channel! receiver is 
found from the expression [2] 





a F, VE 
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Moreover, in a real receiver there will be the error 
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The ratio of the errors for real and ideal reception will be equal to 
AF 
Pe = FTF r. (32) 


If Fy = 0.5/rp, AF = Fray, then P2 = 1, 


Thus it can be considered that a real receiver of PWM-AM signals achieves in practice the potential noise 
stability of this method when the leading edge of the pulse is rigidly fixed and 


Fs =F AF = Finax- 


7. Pulse Width Modulation-Frequency Modulation (PWM-FM) 
Potential noise stability. If it is assumed that the frequency change takes place according to a linear law, 
then the signal of such a system can be represented in the form 
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where t = te + TF + a+ A)s te is the beginning of the leading edge of the pulse. For this signal (considering 
that Vin = 2Vg) we get 
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The mean square measuring error for multi-channel transmission (4 r= T/ Nay) will be determined by the 
expression 
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Real noise stability. The signal-to-noise ratio at the output of a real PWM-FM receiver is determined by the 
expression [2] 
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from which we will find the mean square error 
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The ratio of the errors for reception on real and ideal (per Kotel'nikov) receivers will be equal to 
AF 
Pe = 1.8 FsT rp a . (36) 


If 


P= Fei AF = Fmaythen PL = 1.1. 


Obviously, as with PAM-FM systems, a calculation of the bandwidth according to the formula F, = 1,8/ rr 
gives a somewhat oversized value for F,. 


8. Pulse Width Modulation-Phase Modulation (PWM-PM) 
Potential noise stability. Proceeding from the same premises as for the development of the type PWM-FM 
telemetering system, the signal for a PWM-PM system can be represented in the form 











( Vin Sin (Wot + 9:1) fr 0S t<t, 
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Vmsin [ost +9 — FPO —)] for h<t<t tp 


where = to + Tp +" (1+ A), t, is the beginning of the leading pulse edge. 
Making a change of variables and substituting V,,, =/2Ve. we get 


A’ (A, 0) = V2V_ 2? Ax con[ ust + 9 — FP — tc) for ty St < tet t. 
For the desired integral we have 
1m EOS 
“al 


F 


Moreover, the integration is taken over the limits from —1,/2 to +1 ,/2 because in the remaining time 
A'(A, t)= 0, 








The mean square error when receiving on an ideal receiver will be found from 


2_ N* at, 
Pp sv. oO : (38) 


Real noise stability, To find the error in a real PWM-PM receiver we will use the signal-to-noise ratio at the 
receiver output presented in [2]: V; V3 ®, (A \" Ve 


V. ~ Wa, \F,) VieVar ° 


p T 
biieas a o* Nta® FAF 
r 8? Os 





(39) 





The ratio of errors when receiving on an ideal and a real receiver (for identical signal parameters) will be 


ot APT AF is 
- s*F Fit pe max (40) 


a a 2 
If F,= 0.5/ 7p, AF = Fmax, then P2 = 1, 


Consequently, for the conditions enumerated above in a real receiver of PWM-PM signals it is also possible 
to realize the potential noise stability. 


SUMMARY 
The results of a comparison among various types of multi-channel radiotelemetering systems with time division 
of the channels which was made with respect to real and potential noise stability show that for certain conditions 
the noise stability of almost all systems approaches the potential, 


One of the necessary conditions for the real noise stability of systems with PAM and PWM to approach the 
potential is the use in demodulating of a pulse series of a low-frequency filter having a frequency cutoff equal to the 
maximum frequency of the parameter being measured (AF = F,,,,). The second necessary condition for achieving 


the potential noise stability in the group of systems is the optimum selection of the bandwidth passed by the video 
channel, 


The real noise stability of a PPM-AM type system, even when the conditions given above are observed, is 
worse than the potential by approximately \/2 times. Since real filters cannot have the infinitely great sharpness of 
cutoff, i.e., the width of the noise band (AF) is always greater than F,,,,, then the real noise stability of the systems 
will always be somewhat worse than the potential, 


It should also be noted that when in a radiotelemetering system with time division of the channels a secon- 
dary frequency modulation (PAM-FM, PWM-FM) is used, formula (1) gives a somewhat oversized value of 6 p This 
can be explained by the fact that according to [1] for all pulse modulation systems the noise intensity at the output 


of an ideal receiver has a uniform spectrum, However, a detailed discussion of this question is beyond the scope of 
the article. 


The author expresses deep appreciation to A. P. Manovtsev for valuable advice and interest in the work. 


LITERATURE CITED 
V. A. Kotel’nikov, Theory of Potential Noise Stability (Gosenergoizdat, 1956), 
M. C, Nichols and L, L. Rauch, Radiotelemetry [Russian translation] (Izd, inostr. liter,, 1958), 
V. A. Kashirin, “Noise stability of frequency modulation systems? Avtomatika i telemekhanika, 18, No, 6 (1957). 
V. A. Kashirin and G, A, Shastova, “Noise stability in transmission of telemetering signals over a channel with 
fluctuating noise,” Avtomatika i telemekhanika, 19, No, 8 (1958), 
5. VV. A, Kashirin, “Optimal parameters for noise stability in telemetering systems,” Avtomatika i telemekhanika, 
20, No, 2 (1959), 
6. &, V. Borodich, Multi-Channel Radio Relay Communication Lines {in Russian] (Svyaz‘izdat, 1953), 
7. V. D. Smirnov, Principles of Ultrashort- Wave Radio Communication [in Russian] (Izd, lit, po voprosam svyazi 
i radio, 1957). 





on 
* «© 


> 











DETERMINATION OF DISTURBANCE CHARACTERISTICS 
IN AIRCRAFT ~ AUTOPILOT SYSTEMS 
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A method for determining the mathematical expectation, the autocorrelation function, and 
the disturbance spectrum in aircraft—autopilot systems is presented. The described method 
is also valid for arbitrary linear controlled systems and for automatic control systems. 


the problem of the statistical analysis of disturbance characteristics in a linear aircraft—autopilot system can 
be stated in the following manner. 


The dynamic characteristics of the aircraft and the autopilot and certain given realizations of steady -state 
random processes at their inputs and outputs are assigned. It is necessary to find the equivalent disturbance charac- 
teristics—the mathematical expectation, the autocorrelation function, and the spectrum. 


By the term equivalent disturbance, we shall designate a disturbance representing the over-all action of dis- 
turbances that affect the controlled quantity and are reduced to a certain given point of the system, which is general- 
ly arbitrarily selected. 


Let us consider the flight of an aircraft with an autopilot on a certain course. In the first approximation, the 
linearized equation describing the dynamics of the aircraft — autopilot system can be written in the following form[1]: 
(Jyp? + Myyp + Mi ~=—Q+N, Q= Mi(itkp)¥, (1) 


where ly is the moment of inertia of the aircraft with respect to its vertical axis Y (a right-handed coordinate system 
which is fixed to the aircraft with the coordinate origin at its center of gravity is used); M, = —Q + N is the total 
yawing moment, which acts on the aircraft with respect to the Y 




















v 
M(t) Mylt) | , P axis; My = 3p ry the derivative which characterizes the yawing 
yp? +My Mp on? stability; —_ = ie, is the damping -in-yaw derivative; 
a M* = 3 is the derivative of the moment with respect to the 
rs ¢ rudder deflection; » is the sideslip angle; w, = pp; 5= (1+ kp) » 
My (1p) is the rudder deflection angle, k is the transfer constant in front of 

















the derivative; N(t) is the yawing disturbance, which can be caused 
by wind gusts, air vortices, the effect of banking motion on the 
motion along the flight path, random changes in the air temper- 
ature and pressure along the flight path, etc.; and p = d/dt is the 
differentiation operator. 

The block diagram of the system is shown in the figure. 


In this case, we shall consider the disturbance N(t) to be the equivalent disturbance which acts at the system's 
input. We shall assume that the transfer function of the aircraft 


i 
Mey M°v py + M® 
y+ My" p+ My 











an 


ral- 


n {1}: 


wing 


he 


) » 
of 
sed 


Aue 





and the transfer function of the autopilot 
©, = My (1 + kp). 


are known, 


Moreover, it is assumed that the curve of the yawing moment M, acting on the aircraft is known in a suf- 
ficiently large time interval T +r [2-4], We shall find the mathematical expectation EN(t), the autocorrelation 
function Ry(r), and the spectrum S)(w) of the disturbance N(t) that acts at the system's input. 


For the system under consideration 


N(t) = { M, (¢—) k (0) a0, (2) 
0 


where k(t) is the pulse transfer function which corresponds to the transfer function 1 + @4,. 
By using Eq, (2), we obtain [4] 


Ry() = | Ras we — 0) k(O) a, (3) 
a 


Ryn) =f Ry, + 0) k() dd, «) 
0 


where 


i 


T 
Ry (x) = Fy’ (t) N’ (t — at, 
0 


i 


Ry wl) = x \ My (ON (t — tpt, 


Ry, (*) = a \ My (O My (t— Dat, 


eet ee Mt he 


and N"(t) and M' y(t) are centered inputs, which are determined by the equations 


N'() = N()— EN(), My () = My () — EM, (0, 


where E is the symbol denoting the mathematical expectation, 
The mathematical expectation and the disturbance spectrum are determined by using the equations 


co 
EN (1) = BM, (0) \ k(®) d®, (5) 
0 
co 
Sy) = 2 { Ry (¢) coswt dx. (6) 
0 


We shall assume that My(t) is an ergodic steady -state random process, by the statistical processing of which 
we obtain the following expressions: 


EM, (t) = 0, (7) 
Ry, () = Ce~*'*1, (8) 
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We shall write the pulse transfer function in the following form: 








k(t) = 8 (t) + e* (1 cos at + m sin ot), (9) 
where 5(t) is a delta-function, and 
a 2, — MY k 
Sy “y Viasme —(uen? 





o ¥ 2 
a. w= M4 mM,” . 
—_— a Oe 


By taking into account expressions (7) and (8) and by using Eqs. (4), (3), (5), and (6), we successively obtain the 
following statistical characteristics of the disturbance N(t): 




















EN (t) = 0, (10) 
Ry (*) = ARy (t) + Ce®!*! (F cos at — E sin | +), (11) 
Aa Fe (e* + w* + w?) + Ea,'(e* — w* + @?) (12 
Sy (0) =20 | [e* + (@: + @)*) [e* + (1 — @)*) 
where — 
ma — | (a + e) ma, + | (a — 8) 
A= i i , 
[1+"epeptet Wit @—er ser | 
E = ab— dj, F = db+af, 
ew le mm) _ @+e)— mar l(a —e) + ma, 
4e (e? + wit) ” (a + 8)? + o? (a—e+ a? * 
4e* — 2le* +- mew; + w? (4e — 1) 
” 4s (e* + w?) ; 


fis la, +m(a+e  —lo+m(a—se) 
(a +e) + of (a —e)* + of 





Let us consider an example. Assume that Jy = 450 kg-m-sec’, MYY = 8500 kg-m-sec, MY = 120000 kg-m, 
Mf = 24000 kg-m, k= 1/3 sec, EMy(t) = 0, Ry (7) = ce“OSIT I, 


By using Eqs. (5), (11) and (12), we successively obtain 


EN(t) = 0, 
Ry (8) = Cl A4 5 !*1 4. ¢— 9-441 *1 (9.0872 cos 13.3¢ + 0.0731 sin 13.3 | |], 


1.44 0.29 @* — 950 
Sy (@) = cl arya ~ @— 0.18 wo + 71 wm) 





The case considered above corresponds to an actual dynamic test of an aircraft model with autopilot in an 
aerodynamic wind tunnel. The total moment M,(t) acting on the aircraft was in this case realized by means of 
aerodynamic scales, 


The characteristics of the yawing disturbance N(t) acting on the aircraft can also be found by using the results 
of flight tests if the means for realizing either the sideslip angle y or the rudder deflection angle 5 are available 
to the experimenter. The characteristics of rolling and pitching disturbances can be determined in a similar manner. 
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Wide use is being made today of integrating amplifiers with parallel feedback and with magnetic feedback in 
the case of magnetic amplifiers [1]. To assure accuracy and the possibility of integrating relatively low frequencies, 
the time constant of the system must be sufficiently large. This can be achieved either by increasing the capacitance 
in the feedback loop or by increasing the amplification factor of the amplifier. Each of these methods has its draw - 
backs. 


Increasing the capacitance leads to an increase in the volume of the capacitors and a decrease of their leak- 
age resistance. The latter leads to an increase in the fixed negative feedback and a corresponding decrease in the 
amplification factor, which sometimes nullifies the effect of the increase in capacitance. 


An increase in the amplification factor usually requires an additional amplification stage, which leads to an 
increase in the drift voltage at the amplifier output, as well as to poorer stability of the system. 


Cases are encountered in practice in which both possibilities have been exhausted and it is impossible to in- 
crease any further either the capacitance or the amplification factor, but the time constant of the system must be 
made severa! times as large. 


The problem may also be formulated in another way. It is required to decrease the drift voltage at the out- 
put, at the same time keeping the same value for the time constant of the system. If in order to decrease the drift 
we decrease the amplification factor, then at the same time we decrease the time constant, and the integration errors 
will increase. 


We present below a circuit for an integrating amplifier containing an additional auxiliary amplifier and making 
it possible to some degree to satisfy the above-mentioned contradictory requirements [2], Both amplifiers may be 
either push-pull or single-ended amplifiers. 


The figure below shows circuits for connecting an auxiliary amplifier for the case of parallel feedback, and for 
the case of magnetic feedback when the basic amplifier is magnetic. The amplifiers proper we shall first consider 
to be of electronic type. 


The transfer function for parallel feedback (figure, a) is equal to 


K ub 4 (1) 


W() = TP R/RpT + oF 





where T is the time constant of the system, equal to 


TK ub 


T= TERT R(t Fy thea)’ ? 


kyb and ky are the voltage amplification factors of the basic and auxiliary amplifiers respectively; 1 = RC, For 
I/kyp « (1 + kya) we have 


Pa ag (3) 
TF RR, * * Sua): 





TOrs 


for 





For magnetic feedback the transfer function has the form 
Vour kp gn (1 + PO) 
Sere = uf 
W (p) Tin en +? , (4) 


where kp = Upy;/ Fin is magnetomotive force amplification factor of the basic amplifier; 1 = Rj, C; T is the time 
constant of the system, equal to 


T = (Rp + (1 +kyg) tpwp.1C- (5) 


For | pr| « 1 and Ry, «(1+ kya)kpwp, which is usually the case,we obtain 


kpw 
WO)=TEor T = Cw, kp (1 + kya). (6) 


From expressions (3) and (6) it can be seen that the time constant of the system here is (1 + kya) times as big 
as in the known analogous circuits for integrating amplifiers without auxiliary amplifiers [1}, 


Now we shall consider the questions of the influence of drift, which are very important for integrating am~- 
plifiers. We shall find what the drift voltage Ug, at the output of the basic amplifier willbe when a drift voltage 


Uda 4ppears at the output of the auxiliary amplifier. For this case, for example, with magnetic feedback, the transfer 
function is equal to 


U k pw, pC 
a ee 
Wal?) = 7a thr 





(7) 


From (6) ana ‘7) follows the expression, obvious in this case, for the basic amplifier drift seen at the input: 
Yin =U daa Cp. (8) 


From (7) and (8) and from the circuit of figure b it can be seen that in the steady state a drift voltage Ugg = 
= constant of the auxiliary amplifier will not cause any additional drift in the basic amplifier. Only the drift volt- 
age of. the basic amplifier will appear in the steady state, and this will not increase because of the introduction of 
the auxiliary amplifier. 
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If the drift voltage Uda changes comparatively rapidly, and | pT | >> 1, then the expression (7) will ap- 
proximately be 


re 
i+k, ' 
i.e., in this case the drift at the output caused by the drift voltage of the auxiliary amplifier will be smaller by 4 
factor of 1 + ky,. 


In order to obtain a maximum time constant of the system, the amplification factor ky, should be taken as 
large as possible. However, the value of kgag is limited by the linear range of the characteristic of the auxiliary 


W (p) = (9) 


1279 











amplifier. In practice the maximum obtainable voltage (at the end of the linear zone) at the output of the basic 
amplifier lies between 10 and 100 volts, If we consider the maximum voltage at the output of the auxiliary am- 
plifier to be equal to 100 volts, then it is obvious that the value of k,,, may lie between 10 and 1.0, 


The load current fg of the auxiliary amplifier is the charging (discharging) current of the capacitance. For the 
integration times (frequencies) and capacitance values ordinarily encountered, this current and the output power of 
the auxiliary amplifier are comparatively small, For example, let C= 100 yf, the integration time be equal to 
40 seconds, and the maximum output voltage of the auxiliary amplifier be 100 volts. Then if the nominal signal 
is suddenly connected at the input of the system we find approximately 


ip = M8 ~ 0,25 
5= dt ow Ne ma, 


and the maximum output power is P, ~ 25mw. Thus, the auxiliary amplifier can be comparatively low-power and 
can have a comparatively high output resistance. 


Let us now consider the question of the stability of the above-described integrating system. We shall assume 
that the basic and auxiliary amplifiers are both first-order magnetic elements, with time constants r} and r, res- 
pectively. Then, for example, for the case of magnetic feedback, the transfer function of the system will be equal to 


U i ‘ 
W (p) = 7m 2 ( FPTa)”in (10) 


in i mrt 1 
i +p| 28 + em, « +g) | + 2° E2+ cw ye ) 











From (10) it follows that the given system is stable for the assumptions we have made. Furthermore, it will 
be aperiodic if the condition 





T+ Tt 2 4t 
[Ae 8 + con ct + hua) | > pe ( P+ cwy ). (11) 
is satisfied. 
In practice the following relations are usually observed: 
b+t 
Ey << Oy (1 + ug (12) 
*b 
5, <O- (13) 


Starting from (12) and (13), we obtain instead of (10) the following approximate condition for aperiodicity of 
the system: 


4t 
Cw, (1 +k, > . . (14) 


It should be kept in mind that in the circuit under consideration the capacitor in the feedback network is 
operating at a higher voltage than in a circuit without an auxiliary amplifier. As is known, the leakage resistance 
T; Of some types of capacitors, particularly paper capacitors, decreases as the voltage increases. The presence of a 
finite leakage resistance decreases the time constant of the system in the case, for example, of magnetic feedback, 
by a factor of 





kpwp (1 + kya) 
ie: 
As can be seen from figure a, the auxiliary and basic amplifiers do not have a common ground at the input 
and output. If a tube type auxiliary amplifier is used, this may cause the familiar technical difficulties, If mag- 
netic amplifiers are used, there will of course be no such difficulties. However, if we consider that the amplifica- 





1 to 


3) 


of 


4) 





tion factor ky, is generally no greater than 10, it is obvious that even the tube type amplifier will work sufficiently 
well. Furthermore, if the connection of the auxiliary amplifier output is changed (see the dashed line on figure a), 
we will have the aforementioned common ground, while retaining the general principle of the integrating circuit 
just considered. The difference in this case will be merely that we should substitute ky, for 1 + ky, in the corre- 
sponding formulas mentioned above. 


Experimental investigation of the above-described integrating amplifier completely corroborated the assump- 
tions and relations mentioned in the article, both qualitatively and quantitatively. We investigated circuits, both 
with tube amplifiers and with magnetic amplifiers, with parallel and magnetic feedbacks. We give as an example 
some average parameters of one of the circuits in which we used magnetic amplifiers: kp = 500, wh = 5,000, ky, = 
= 5, C= 60 uf, 7} = 0.25 second, r, = 0.05 second. The introduction of the auxiliary amplifier increased the 
time constant practically by a factor of 6 and made it approximately equal to 800 seconds. The inequalities (12)- 
(14) are satisfied to such a degree that the time constants of both amplifiers could be considerably greater. 


Finally, we wish to emphasize once more that the integrating circuit with an auxiliary amplifier is only supple- 
mentary to the known integration circuits, and its use is advisable in the cases indicated at the beginning of this ar- 
ticle. 
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